Displaying 1 result from an estimated 1 matches for "arhnc".
Did you mean:
argc
2009 Feb 02
0
repeated measures with gls
...evels. The
fertilization dates are equal for all subjects (blocks) but not
periodical (-46, 10, 24, 53, 123 days from the event).
The code that I've been using is:
fit.csnC<- gls(dno.C ~ block+Fertil*factor(Samp), data=dDakot,
method="ML", corr=corCompSymm(, form=~1|block))
fit.arhnC<-update(fit.csnC, corr=corAR1(, form=~1|block),
weight=varIdent(form =~1|Samp))
I've plotted the model with :
plot(fit.csnC)
plot(fit.csnC,dno.C ~ fitted(.), abline = c(0,1))
and residuals look fine, same as observed vs. fitted values.
I would like to make sure that the model specification...