similar to: nlme w/no groups and spatially correlated residuals

Dear R users, I am using lme and nlme to account for spatially correlated errors as random effects. My basic question is about being able to correct F, p, R2 and parameters of models that do not take into account the nature of such errors using gls, glm or nlm and replace them for new F, p, R2 and parameters using lme and nlme as random effects. I am studying distribution patterns of 50 tree

Hello, Sorry to be a bit longwinded, but I've struggled quite a bit with the following over the last few days. I've read all entries related to spatial autocorrelation in R help and haven't found what I'm after. If it's okay, I'm going to first describe my general understanding of the process by which a mixed model can account for correlated errors. If possible, please

Dear R users, I want to explain binomial data by a serie of fixed effects. My problem is that my binomial data are spatially correlated. Naively, I thought I could found something similar to gls to analyze such data. After some reading, I decided that lmer is probably to tool I need. The model I want to fit would look like lmer ( cbind(n.success,n.failure) ~ (x1 + x2 + ... + xn)^2 ,

Hi, I've searched online, in a few books, and in the archives, but haven't seen this. I believe that xerror is scaled to rel error on the first split. After fitting an rpart object, is it possible with a little math to determine the percentage of true classifications represented by a xerror value? -seth -- View this message in context:

Hi, I have probability estimates from two predictive models. I have these estimates and also a binary outcome for a validation data set not used in calibrating either model. I would like to calculate the Brier score for both models on this binary outcome and test the hypothesis that the Brier scores are equal from the two models. I have not been able to find an R function to do this, can

Hi, Let's say I have a data frame (called "example") with numeric values stored (columns V1 and V2). I also have a string variable storing this name x1<-"example" Is there a way to use the variable x so that R knows that I want the specified action to occur on the data frame? For example, summary (x) would return a summary of the data frame? I am considering this

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

Hi, I am reading a large space-delimited text file into R (41 columns and many rows) and need to do run each row's values through another R object and then write to another text file. So, far using readLines and writeLines seems to be the best bet. I've gotten the data exchange working except each row is read in as one 'chunk', meaning the row has all values between two quotes

This works fine: fit42<-nls(Vfs~SSlogis(Months,Asym.Int+Asym.Group*Groupdum,xmid,scal), data=df, start=c(Asym.Int=22,Asym.Group=5,xmid=2,scal=6), na.action=na.omit) But this, identical except using gnls, doesn't converge: fit43<-gnls(Vfs~SSlogis(Months,Asym.Int+Asym.Group*Groupdum,xmid,scal), data=df, start=c(Asym.Int=22,Asym.Group=5,xmid=2,scal=6), na.action=na.omit) Error in gnls(Vfs

dear R users, I have troubles performing Moran.I test as suggested on http://www.ats.ucla.edu/stat/r/faq/morans_i.htm my spatial data are longitude and lattitide of communities. The calculation of the inverse distance matrix according to the homepage (using my data) datAL <- read.csv2("C:\\Konvergenz AL.csv", header=T) ALdist <- as.matrix(dist(cbind(datAL$L?nge,

Hi all, I'm running into some computer issues when trying to run a binomial model for spatially correlated data using glmmPQL and was wondering if anyone could help me out. My whole dataset consists of about 300,000 points for which I have a suite of environmental variables (I'm trying to come up with a habitat model for a species of seal, using real (presence) and simulated dives

Hi, I have spatially autocorrelated data (with a binary response variable and continuous predictor variables). I believe I need to do an autologistic model, does anyone know a method for doing this in R? Many thanks C Bell

Dear all, I am trying to gain some experience with the function gnls from the nlme package. I tried to model the Theophyline data by trying to model the presumed dependency of the clearance on the body weight. This is my function call of gnls: gnls(conc~SSfol(Dose,Time,lKe,lKa,lCl),data=Theoph, params=list(lKe~1,lKa~1,lCl~Wt),start=c(-2.4,0.46,-3.22,0.01)) That's been the result: Error

Greetings, May you please suggest a package or function to use for fitting a GLMM (generalized linear mixed model) with spatially correlated random effects? Thank you, Elijah DePalma [[alternative HTML version deleted]]

Hi, I guess this is a final plea, and maybe this should go to R-help but here goes. I am writing a set of functions for calibration and prediction, and to calculate standard errors and intervals I need the variance function to be evaluated at new prediction points. So for instance fit<-gnls(Y~SSlogis(foo,Asym,xmid,scal),weights=varPower())

Hi All, I have question about speficifying a constant in gnls() from package nlme. Here is a testing code: ############# library(nlme) x = exp( rnorm(100)) y = 1/(1+x) + rnorm(100)/10 plot( y ~ x) fm1 = gnls( y ~ 1/(1+(x/v)^w), start=list( v=1, w=1)) a =1; b=1; fm2 = gnls( y ~ a/(b+(x/v)^w), start=list( v=1, w=1)) #This won't work because I don't know to set $a$ and $b$ as

Hi List, does anyone know of an implementation of spatially constrained clustering in R? This is where there is a vector of measurements for points on a plane and only neighbors can be clustered together. I have tried implementint in myself -- but if someone has alkready done it ! I have searched on the obvios terms "spatially constrained clustering" without any luck.

Hi Everyone, I am doing a project based on "Spatially Lagged Predictor Variable Models", I would like to know which package in R would execute this model. Also, I am new to this field of spatial statistics. Any suggestions for a good book on spatial regression analysis would be appreciated. Thanks Again. Cheers Arun -- View this message in context:

>>>>> "VW" == Viechtbauer Wolfgang (STAT) <Wolfgang.Viechtbauer at STAT.unimaas.nl> >>>>> on Fri, 26 Sep 2008 18:00:19 +0200 writes: VW> Hi all, I'm trying to fit a marginal (longitudinal) VW> model with an exponential serial correlation function to VW> the Orange tree data set. However, R crashes frequently VW>

Hi I use gnls to fit non linear models of the form y = alpha * x**beta (alpha and beta being linear functions of a 2nd regressor z i.e. alpha=a1+a2*z and beta=b1+b2*z) with variance function varPower(fitted(.)) which sounds correct for the data set I use. My purpose is to use the fitted models for predictions with other sets of regressors x, z than those used in fitting. I therefore need to