# plot growth curves per subject, pooled over conditions # HQ 20120218 sort(unique(TimeC[sel1])) -> xx FE02 <- fixef(kb07.m02) FE02fit <- FE02[1]+FE02[2]*xx+FE02[3]*xx^2 plot( xx+350,(FE02fit), xlab="Time (ms)", ylab="log odds of gaze-on-target", ylim=c(-2.5,0), type="n" ) RE02 <- ranef(kb07.m02)$SubjID for (j in 1:56) { # note different sorting orders in RE and FE parts y <- FE02fit + RE02[j,3]+RE02[j,2]*xx+RE02[j,1]*xx^2 lines( xx+350, y, col="grey" ) } points(xx+350,(FE02fit), pch=16, cex=1.5 )