[R-sig-ME] Plot augPred

Fernanda Balbino ferbalbino at yahoo.com.br
Tue Sep 6 19:49:06 CEST 2011


Hi,

I have this code 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 now I need to plot augPred and I'm trying this:

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

But it doesn't work. 
I need help to plot the augPred.

Thanks for your help.
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/20110906/b82467f6/attachment-0003.txt>


More information about the R-sig-mixed-models mailing list