Sorry, my fault, the model incorporate centered version of the coordinates
so you should read

test_lme <- lme(Dens ~ Speed , random = ~ 1|Year/Date/Survey, data = test,
correlation = corGaus(form = ~ Xprojc + Yprojc))

Arnaud

2010/12/15 Arnaud Mosnier <a.mosnier@gmail.com>

> Thanks for your answer Ben !
>
> So since my last mail, I discovered the possibilities offered by specifying
> the spatial structure in lme using the "correlation" option*.
> *However, when I define the correlation structure with a function like
> corGaus, things start correctly, but after a moment, R simply close without
> any warning on my desktop computer.
> I used the same code on a workstation with larger RAM capacities and Rx64
> and this time it did not close but gave me the following message.
> "Error: evaluation nested too deeply: infinite recursion /
> options(expressions=)?"
>
> Looking at previous mails dealing with that kind of problem, I tried
> options(expressions = 10000), but this time R close without warning as
> previously on the other computer.
> Have you suggestions about what could cause the infinite recursion ?
>
> my model has the following form
>
> test_lme <- lme(Dens ~ Speed , random = ~ 1|Year/Date/Survey, data = test,
> correlation = corGaus(form = ~ Xproj + Yproj))
>
> >str(test)
>
> 'data.frame':   54445 obs. of  7 variables:
>  $ Xprojc: num  -1.91 -1.91 -1.91 -1.91 -1.91 ...
>  $ Yprojc: num  -1.99 -1.96 -1.93 -1.9 -1.87 ...
>  $ Year  : Factor w/ 11 levels "1990","1992",..: 1 1 1 1 1 1 1 1 1 1 ...
>  $ Date  : Factor w/ 34 levels "1990-09-12","1992-09-12",..: 1 1 1 1 1 1 1
> 1 1 1 ...
>  $ Survey: Factor w/ 37 levels "pts_kernP1990-09-12_Real",..: 1 1 1 1 1 1 1
> 1 1 1 ...
>  $ Dens  : num  4.72e-05 8.27e-04 4.65e-03 8.70e-03 5.65e-03 ...
>  $ Speed : num  0.794 0.575 0.572 0.67 0.867 ...
>
> Note: I am working in a windows XP environment
>
> Arnaud
>
>
>
>
>
> Date: Tue, 14 Dec 2010 10:13:41 -0500
> From: Ben Bolker <bbolker@gmail.com>
> To: r-sig-mixed-models@r-project.org
> Subject: Re: [R-sig-ME] Spatial autocorrelation in dependent variable
> Message-ID: <4D0789A5.4040705@gmail.com>
> Content-Type: text/plain; charset=ISO-8859-1
>
> -----BEGIN PGP SIGNED MESSAGE-----
> Hash: SHA1
>
>
> On 10-12-14 10:01 AM, Arnaud Mosnier wrote:
> > Dear list,
> >
> > I am searching a method to include spatial autocorrelation structure of
> my
> > dependent variable into a mixed model.
> > Any suggestions would be much appreciated ?
> >
>
>  see the 'correlation' argument of ?lme in the nlme package.
> Furthermore, see ?corSpatial (and the appropriate chapter of Pinheiro
> and Bates 2000)
>
>  Ben Bolker
> -----BEGIN PGP SIGNATURE-----
> Version: GnuPG v1.4.10 (GNU/Linux)
> Comment: Using GnuPG with Mozilla - http://enigmail.mozdev.org/
>
> iEYEARECAAYFAk0HiaUACgkQc5UpGj
> wzenMB5QCeMFkKEtlm40idYcH+jOPpOyQn
> /fIAnjNq2IeaqYke/xjfWaLtPRYv/9/U
> =okbF
> -----END PGP SIGNATURE-----
>

	[[alternative HTML version deleted]]


