[R] spatial corStruct in lme

From: Lyndon Brooks <lbrooks_at_scu.edu.au>
Date: Wed 07 Jun 2006 - 09:35:02 EST


I'm fitting a relatively simple growth model to some forest plot data. Two species of trees were planted in different mixtures in 10 (nearly-adjacent) plots and measured on four occasions over 10 years. The model is constructed in terms of the diameter increments (per year; DI) in the 3 intervals, in which DI is modelled as a function of mid-interval D and DSQ. The details of the fixed part of the model are not so important here, but four pertinent variables are distance-dependent competition indices: C_spp1fromspp1, C_12, C_21, C_22.

The random structure is: random = ~ 1 | PLOT/TREE.

More complex random structures aren't required, there's no obvious heterogeneity, and the serial correlations (corSymm) are trivial.

I've been trying to fit a spatial correlation structure using the X, Y coordinates of each tree. I've obviously missed a point or so here.

It seems to me that such a structure could be fit for all trees (both species) on their average growth over the 3 intervals, by interval, or by species by interval.

For now, it would be step forward to get any of these working.

Apart from me learning how to do it, the import of this is that, if the CI's are doing what I hope they might, there should be little residual spatial correlation.



R-help@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html Received on Wed Jun 07 10:56:59 2006

Archive maintained by Robert King, hosted by the discipline of statistics at the University of Newcastle, Australia.
Archive generated by hypermail 2.1.8, at Wed 07 Jun 2006 - 12:10:36 EST.

Mailing list information is available at https://stat.ethz.ch/mailman/listinfo/r-help. Please read the posting guide before posting to the list.