[R-sig-ME] Fitting a fixed effect model while modelling correlation in the residuals in 2 directions

[Adam D. I. Kramer]

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