[R-sig-ME] Logistic regression with spatial autocorrelation structure

[Ben Bolker]

-----BEGIN PGP SIGNED MESSAGE-----
Hash: SHA1

On 11-01-17 09:03 AM, Arnaud Mosnier wrote:
> Dear list,
> 
> Is there a way in R to do a mixed logistic regression with a spatial
> autocorrelation structure ?
> 

  The only straightforward solution I know of is via glmmPQL (MASS
package), although you should be careful (e.g. try fitting to simulated
data where you know the answer first) because penalized quasi-likelihood
is biased for binary data.
   One of the generalized estimating equation packages (geepack, etc.)
may also work.  If you have only spatial autocorrelation (i.e. no random
block effects) then the geoRglm and spatcounts packages may also be of use.

  Ben Bolker
-----BEGIN PGP SIGNATURE-----
Version: GnuPG v1.4.10 (GNU/Linux)
Comment: Using GnuPG with Mozilla - http://enigmail.mozdev.org/

iEYEARECAAYFAk00UzEACgkQc5UpGjwzenMJ4wCeJxVyJLDWXP2hxUZ9MpXxeTsN
pxwAn1t/pb1eduoBg1+WF480SzjvLAqz
=yvA0
-----END PGP SIGNATURE-----