pc.asr <- asreml(W1000~1+ + NIRday, random=~sname2+Millday+ idv(ROW):ar1(COL),
,data=data.df)
this model converges without trouble. The model with ar1v(ROW):COL did not converge
The REML estimates are silly ( ie corr of -.76) but anyway this gives the syntax of how to do this in ASReml-R.
Workflow:
Published 

