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