similar to: lme, lmer, gls, and spatial autocorrelation

I'd like to compare two models which were fitted using gls, however I'm having trouble interpreting the results of gls. If any of you could offer me some advice, I'd greatly appreciate it. Short explanation of models: These two models have the same fixed-effects structure (two independent, linear effects), and differ only in that the second model includes a corExp structure for

A couple weeks ago I posted a message on this topic to r-help, the response was that this seemed like odd behavior, and that I ought to post it to one of the developer lists. I posted to r-sig-mixed-models, but didn't get any response. So, with good intentions, I decided to try posting once more, but to this more general list. The goal is (1) FYI, to make you aware of this issue, in case it

Dear all, I am experiencing problems using the glmmPQL function in the MASS package (Venables & Ripley 2002) to model binomial data with spatial autocorrelation. My question - is the presence of birds affected by various hydrological parameters? Presence/absence data were collected from 83 sites and coupled against hydrological data from the same site. The bird survey sampling effort

I have spatial data on a sphere (the Earth) for which I would like to run an gls model assuming that the errors are autcorrelated, i.e. including a corSpatial correlation in the model specification. In this case the distance metric should be calculated on the sphere, therefore metric = "euclidean" in (for example) corSpher would be incorrect. I would be grateful for help on how to

Hi. I'm trying to perform what should be a reasonably basic analysis of some spatial presence/absence data but am somewhat overwhelmed by the options available and could do with a helpful pointer. My researches so far indicate that if my data were normal, I would simply use gls() (in nlme) and one of the various corSpatial functions (eg. corSpher() to be analagous to similar analysis in SAS)

Hi, I have two sets of data for a given set of (non-lattice) locations. I would like to know whether the two are significantly different. This would be simple enough if it wasn't for the fact that the data is spatially autocorrelated. I have come across several possible solutions (including Cliff & Ord which however appears to be for gridded data), or using gls. However, they don't

Hi, I am trying out a generalized least squares method of forecasting that corrects for autocorrelation. I downloaded daily stock data from Yahoo Finance, and am trying to predict Close (n=7903). I have learned to use date functions to extract indicator variables for Monday - Friday (and Friday is missing in the model to prevent it from becoming full rank). When I run the following code...

Dear all, I am currently facing a problem related to the spatial autocorrelation of a sample of stations; these stations supply weekly data for a fixed time-window during the year (namely, 4-6 months per year). For this reason I'm trying to use the R package 'spdep' (specifically Moran's I) in order to get rid of it. Does anyone know how is it possible (if it is...) to

Hi R helpers, Is there a function (package?) in R available which tests "spatial autocorrelation" between points (e.g. vector layer of weather stations)? (e.g. Moran's I...) Via the archives we found out that there is a package 'spdep' which uses grid data for testing spatial autocorrelation. Thanks a lot, Jan

Hi, I have a problem with a regression I try to run. I did an estimation of the market model with daily data. You can see to output below: /> summary(regression_resn) Time series regression with "ts" data: Start = -150, End = -26 Call: dynlm(formula = ror_resn ~ ror_spi_resn) Residuals: Min 1Q Median 3Q Max -0.0255690 -0.0030378 0.0002787

Howdy All, I am looking for some good tutorials (books, websites, whatever) for calculating/testing for Spatial Autocorrelation using R. Specifically, I am wanting to test for autocorrelation of a number of variables measured at a set of discrete locations. Up to this point I have been exploring the "spdep" package and I can get "moran.test" to work, but I am concerned that

Dear list members, I have a large dataset with point data. Each data point has X, Y, number of species present and values for one single explanatory variable. Each point was generated from a polygon grid shapefile and they are at equal distances. I conducted GLS analyses introducing variance structures and correlation structures to improve the model. I plotted variograms before and after adding

Hi, I have collected data on trees from 5 forest plots located within the same landscape. Data within the plots are spatially autocorrelated (calculated using Moran's I). I would like to do a ANCOVA type of analysis combining these five plots, but the assumption that there is no autocorrelation in the residuals is obviously violated. Does anyone have any ideas how to incorporate these spatial

Hello, I have a matrix with the distances among sites. And I have another matrix with the presence and absence of each species in each site. I would like to test the spatial autocorrelation among sites. I have tried to use the function gearymoran of the ade4 package, but error messages keep popping up. Do you know any function for me to test the spatial autocorrelation of my data? Thanks,

Dear all: Does anyone have any suggestions on how to make a spatially explicit landscape with spatial autocorrelation in R? In other words, a landscape where all cells have a spatial reference, and the environment values that are closer in space are more similar (positive spatial autocorrelation). Thank you, Laura

Hello, Has anyone modified the errorsarlm in the R package spdep to allow for more than a single spatial autocorrelation coefficient (i.e. 'lambda')? Or, if not, any initial suggestions on how to make that modification? I have looked at the source code for the function and realize that any attempt to do it on my own would require much dedication, so would like to

Hi, I am new to this list. I have a question regarding including both spatial and temporal random factors in lmer. These two are not nested, and an example of model I try to fit is model1<-lmer(Richness~Y+Canopy+Veg_cm+Treatment+(1|Site/Block/Plot)+(1|Year), family=poisson, REML=FALSE), where richness = integer Y & Treatment = factor Canopy & Veg_cm = numerical, continous

Dear all: I would like to create a landscape of environmental values that follow a uniform frequency distribution and also have spatial autocorrelation in the landscape. I was wondering if there is an algorithm and/or package out there that creates autocorrelation of values that are distributed according to a non-normal frequency distribution. Any suggestions are greatly appreciated. Thank you,

Dear Users: I have been trying to use the geary.test() function in *R*, but am having slight difficulty understanding how I am to apply it in my context. I have 2 matrices: 1) *n x p* matrix of *n* observations with *p* measurements each. It may be noted that this matrix has a spatial dimension to it, as the *n*observations are at different geographical locations on a map. 2) *n x n* spatial

I'm trying to do a simple (?) analysis of a 1D spatial data set, allowing for spatial autocorrelation. (Actually, I'm comparing expected vs. observed for a spatial model of a 1D spatial data set.) I'm using models like gls(obs~exp,correlation=corExp(form=~pos),data=data) or gls(obs~exp,correlation=corLin(form=~pos),data=data) This form is supposed to fit a linear model of