Displaying 1 result from an estimated 1 matches for "dayspt3".
Did you mean:
dayspt
2009 May 18
1
Predicting complicated GAMMs on response scale
...,lab=c(24,4,12),las=1,cex.lab=1.5, cex.axis=1,lwd=2)
lines(p.d$DaysPT,b$fit+b$se.fit*1.96,lty=2,lwd=1.5)
lines(p.d$DaysPT,b$fit-b$se.fit*1.96,lty=2,lwd=1.5)
points(DaysPT,Diff)
However, when I add a correlation structure and/or a variance structure so
that the model may look like:-
gamm(Diff~s(DaysPT3)+AirToC,correlation=corCAR1(form=~DaysPT|Animal),weights=varPower(form=~DaysPT),method="REML")
I get this message at the point of inputting the line
"b<-predict.gam(g.m$gam,p.d,se=TRUE)"
Error in model.frame(formula, rownames, variables, varnames, extras,
extranames, :...