Displaying 1 result from an estimated 1 matches for "varslop".
Did you mean:
varslope
2009 Apr 01
3
How to prevent inclusion of intercept in lme with interaction
...tment is not included (newbone~t-1), but
seems to have no effect in (newone ~ t*treat-1).
What's wrong? Some bracket missing? I tried a few variants.
Dieter Menne
#--------------------------------
library(nlme)
library(lattice)
# Generated data
set.seed(4711)
subject = as.factor(letters[1:5])
varslope = rnorm(length(subject),0,0.02)
cslope = c (0.1,0.15)
grd = expand.grid(t=seq(5,15,by=5),
subject=subject,treat=c("contr","test"))
grd$slope = varslope[grd$subject] + cslope[grd$treat]
grd$newbone = grd$slope*grd$t+rnorm(nrow(grd),0,0.2)
xyplot(newbone~t|treat,groups=sub...