Dear R users,

I'm estimating a mixed effects model in which the spatial correlation is
controlled for by the "corGaus" structure. I'm wondering if there is a
document or paper that explains how the spatial correlation structure (such
as "corExp" or "corGaus") works.

Let me use the example and data posted on UCLA's R FAQ webpage to explain
my problems. The link for the webpage is:
http://www.ats.ucla.edu/stat/r/faq/spatial_regression.htm

install.packages("nlme")
library(nlme)

spdata <- read.table("http://www.ats.ucla.edu/stat/R/faq/thick.csv", header
= T, sep = ",")
dummy <- rep(1, 75)
spdata <- cbind(spdata, dummy)

### estimate the null model ###
soil.model <- lme(fixed = thick ~ soil, data = spdata, random = ~ 1 |
dummy, method = "ML")
summary(soil.model)
plot(Variogram(soil.model,form=~north+east))

### updated the model by the spatial correlation structure ####
soil.gaus <- update(soil.model, correlation=corGaus(1,form=~north+east))
summary(soil.gaus)
plot(Variogram(soil.gaus,form=~north+east))

My questions are:

1) Is there a way that I can tell, to what extent the spatial correlation
has been controlled for by the model, besides the improvements of AIC, BIC,
and Loglik?

2) where can I find the formulas for the corExp or corGaus?

Thanks!

Gary

	[[alternative HTML version deleted]]

