Hi, thanks to all for the helpful information. In fact our colleague is indeed
keen to check for the difference. More importantly, I was in too much a haste
to drop the message as I now recalled we got around in a similar way in Stata
and SAS for a more extreme case before. JH
> Date: Thu, 7
[EMAIL PROTECTED] wrote:
> Full_Name: Jari Oksanen
> Version: 2.6.2 RC (2008-02-07 r44369)
> OS: Linux
> Submission from: (NULL) (130.231.102.145)
>
>
> Even if function xspline() is called with argument draw=FALSE, it requires a
> graphics device (that it won't use since it was draw=FALSE). I run
[EMAIL PROTECTED] wrote:
> Full_Name: Jari Oksanen
> Version: 2.6.2 RC (2008-02-07 r44369)
> OS: Linux
> Submission from: (NULL) (130.231.102.145)
>
>
> Even if function xspline() is called with argument draw=3DFALSE, it req=
uires a
> graphics device (that it won't use since it was draw=3DFALSE).
Full_Name: Jari Oksanen
Version: 2.6.2 RC (2008-02-07 r44369)
OS: Linux
Submission from: (NULL) (130.231.102.145)
Even if function xspline() is called with argument draw=FALSE, it requires a
graphics device (that it won't use since it was draw=FALSE). I run into this
because I intended to use xs
library(glmpath)
data(heart.data)
# heart.data is a list, $y a vector, $x a matrix
data <- data.frame(x=I(heart.data$x), y = heart.data$y)
> data[1:2,]
x.1 x.2 x.3 x.4 x.5 x.6 x.7 x.8 x.9 y
1 16012 5.73 23.11 149 25.3 97.252 1
2 144 0.01 4.41 28.61 0
Consider the following:
> lm(conc ~ Type, CO2, subset = Plant == "Quebec")
Error in `contrasts<-`(`*tmp*`, value = "contr.treatment") :
contrasts can be applied only to factors with 2 or more levels
Here Type is a factor for which only one level occurs
within the "Quebec" subset. This is some
