Fitting a fixed effect model while modelling correlation in the residuals in 2 directions
On Tue, 15 Sep 2009, Thomas Bishop wrote:
slate.ar1*ar1 <- gls(Yield~as.factor(Variety),corr = corAR1(form = ~ 1|(FieldRow+FieldColumn),data=Slate))
...data=Slate is one paren too deep. Should be ...+FieldColumn)),data=Slate) :) --Adam