[R] [R-sig-ME] Xy plot help

[Jeff Newmiller]

A) Yes, there is something wrong. Code not reproducible. Behavdat not found.

B) Last time I used it, Yahoo Mail could send in plain text. Admittedly that has been awhile, but Google sez read [1] and/or [2].

[1] http://help.yahoo.com/tutorials/mail/mail_persmsg1.html
[2] http://email.about.com/od/yahoomailtips/qt/et/plain_text.htm
---------------------------------------------------------------------------
Jeff Newmiller                        The     .....       .....  Go Live...
DCN:<jdnewmil at dcn.davis.ca.us>        Basics: ##.#.       ##.#.  Live Go...
                                      Live:   OO#.. Dead: OO#..  Playing
Research Engineer (Solar/Batteries            O.O#.       #.O#.  with
/Software/Embedded Controllers)               .OO#.       .OO#.  rocks...1k
--------------------------------------------------------------------------- 
Sent from my phone. Please excuse my brevity.

arun <smartpink111 at yahoo.com> wrote:

>Dear R group,
>
>Tried to turn off the html in yahoo.� So far� not successful.�
>Apologies!
>
>I tried to get xyplot shrink fit for my mixed model.� But, the mixed
>model line is not seen in the graph.
>I would like to know if there is anything wrong with my code.
>Thanks,
>A.K.
>
>
>
>(fm1 <- lmer(Response3 ~1+ Wavelength*Start_Resp*time + (1|resid) +
>(1+time|Subject_BDat), family=binomial, 
>data=Behavdat, REML=0))
>
>df <- coef(lmList(Response ~ time | SubjectBDat, Behavdat)) 
>fclow <- subset(df, '(Intercept)' < -4.1) 
>fchigh <- subset(df, '(Intercept)' > -4.1) 
>cc1 <- as.data.frame(coef(fm1)$Subject) 
>names(cc1) <- c("A", "B","C","D","E","F","G","H","I","K","L","M") 
>df <- cbind(df, cc1) 
>ff <- fixef(fm1)
>SubjectBDat<-as.factor(Behavdat$Subject)
>print(xyplot(Response ~ time | SubjectBDat, Behavdat, 
>������������� type = c("g", "p", "r"), 
>������������ coef.list = df[,3:14], 
>������������ panel = function(..., coef.list) { 
>���������������� panel.xyplot(...) 
>���������������� panel.abline(as.numeric(coef.list[packet.number(),]), 
>����������������������������� col.line =
>trellis.par.get("superpose.line")$col[2], 
>����������������������������� lty =
>trellis.par.get("superpose.line")$lty[2] 
>����������������������������� ) 
>���������������� panel.abline(fixef(fm), 
>����������������������������� col.line =
>trellis.par.get("superpose.line")$col[4], 
>����������������������������� lty =
>trellis.par.get("superpose.line")$lty[4] 
>����������������������������� ) 
>������������ }, 
>������������ index.cond = function(x,y) coef(lm(y ~ x))[1], 
>������������ xlab = "time", 
>������������ ylab = "Response", 
>������������ key = list(space = "top", columns = 3, 
>������������ text = list(c("Within-subject", "Mixed model",
>"Population")), 
>������������ lines = list(col =
>trellis.par.get("superpose.line")$col[c(1:2,4)], 
>������������ lty = trellis.par.get("superpose.line")$lty[c(1:2,4)]))))
>	[[alternative HTML version deleted]]
>
>_______________________________________________
>R-sig-mixed-models at r-project.org mailing list
>https://stat.ethz.ch/mailman/listinfo/r-sig-mixed-models