[R-sig-ME] Help to Plot augPred

[Fernanda Balbino]

Dear list,


I need help to plot the graph augPred of a lme adjustment. I'm using the code below for adjust lme for the data "dados":

da <- read.table("dados.txt", header=TRUE, sep="\t", dec=",")
require(nlme)
ajuste <- lapply(split(da, da$Temp), function (a){
mm0 <- lme(modulo~sqrt(Freq)+I(Freq^2), random=~1|bloco,data=a)
mm0
})
lapply(ajuste, summary)
lapply(ajuste, ranef)

But when I try to plot augPred the R program send me an error message, I'm trying this:

plot(augPred(ajuste$Temp[[1]], level=c(0,1))


Can you help me find the right way to do the plot using augPred for this case?

Thank you in advance for any help,
Regards,
Fernanda.
-------------- next part --------------
An embedded and charset-unspecified text was scrubbed...
Name: dados.txt
URL: <https://stat.ethz.ch/pipermail/r-sig-mixed-models/attachments/20110912/934d4a46/attachment-0003.txt>