I have searched for many years for examples of fitting AR1xAR1 type models
with open-source software other than ASREML.  There are no such examples.
Only ASREML (and SAS PROC MIXED) can fit these models.

Kevin

On Tue, Nov 1, 2011 at 6:27 AM, John Clark <rosbreed.pba@gmail.com> wrote:

> Dear mixed model experts:
>
> I want to fit mixed model using lme4, I need something equivalent to the
> asreml-R codes in lme4 or other free mixed model packages :
>
>    data(nin89)
>    # Model 1: RCB analysis with G and R structure
>     rcb.asr <- asreml(yield ~ Variety, random = ~ idv(Replicate), rcov = ~
> idv(units), data = nin89)
>
>     # Model 2: two-dimensional spatial model with correlation in one
> direction
>      sp.asr <- asreml(yield ~ Variety, rcov = ~ Column:ar1(Row), data =
> nin89)
>
>    # model 3
>     sp.asr <- asreml(yield ~ Variety, rcov = ~ idv(Column):ar1(Row), data
> = nin89, family = asreml.gaussian(dispersion = 1.0))
>
>    # model 4: two-dimensional spatial model
>     sp.asr <- asreml(yield ~ Variety, rcov = ~ ar1(Column):ar1(Row), data
> = nin89)
>       sp.asr <- asreml(yield ~ Variety, random = ~ units, rcov = ~
> ar1(Column):ar1(Row), data = nin89)
>
> Sorry for cross posting, however I did not get any answer in the following
> post and was suggested to email here  to harvest the expertise.
>
>
> http://stackoverflow.com/questions/7961864/r-lme4-code-equivalent-to-asreml-r
>
> Thanks
>
> John
>
>        [[alternative HTML version deleted]]
>
> _______________________________________________
> R-sig-mixed-models@r-project.org mailing list
> https://stat.ethz.ch/mailman/listinfo/r-sig-mixed-models
>



-- 
Kevin Wright

	[[alternative HTML version deleted]]


