Displaying 1 result from an estimated 1 matches for "yy4".
Did you mean:
yy
2011 Mar 19
1
strange PREDICTIONS from a PIECEWISE LINEAR (mixed) MODEL
...;fictious",
x,
z)
#predict
yy3<-predict(mod3,level=0, newdata=newframe2)
lines(x[order(x)],yy3[order(x)],col="green",lwd=2)
# ADD INTERACTION z:x
mod4<-lme(y~x+x*(x>-1)+
z+
z:x+
z:x*(x>-1)
,random=~x|id,
data=data.frame(x,y,z,id))
#predict
yy4<-predict(mod4,level=0, newdata=newframe2)
lines(x[order(x)],yy4[order(x)],col="violet",lwd=2) #something bizarre
#starts to happen
#in the predicted values
# they begin to jiggle...