[R-sig-ME] UNSOLVED: Re: lme4 or other open source mixed model package code equivalent to asreml-R

Ben Bolker bbolker at gmail.com
Thu Nov 24 23:31:56 CET 2011

David Duffy <David.Duffy at ...> writes:

> On Fri, 18 Nov 2011, Ben Bolker wrote:


> Well, I do not have any experience carrying out spatial modelling, but it 
> did seem to me that the covariances arising from the separable AR1xAR1 
> model must be have a close resemblance to some of the other standard 
> models.  Specifically, page 85 of Haskard's thesis
> http://digital.library.adelaide.edu.au/dspace/bitstream/2440/47972/1/02whole.pdf
> seems to imply it is equivalent to a anisotropic Matern model.  This is 
> offered by the spatialCovariance package, AIUI.
> It took me a while to understand that the example dataset is the Wheat2 
> data from the nmle package, which Pinheiro and Bates (2000, P 263) analyse 
> using spherical and rational quadratic models.  I am unsure if the OP is 
> particularly interested in the exact models he chose for his example, or 
> whether one can generally match ASREML's functionality, or has data that 
> might be comfortably analysed usibg existing lme correlation structures.

  Dave Fournier has posted a solution using AD Model Builder at


(sorry, URL broken to make Gmane happy) that tackles (almost) this problem (it
does a spatial GLMM -- a Poisson model layered on the AR1 x AR1
Gaussian random field).

  Ben Bolker

More information about the R-sig-mixed-models mailing list