similar to: Specifying spatial correlation Form in nmle

Hi, I have been using the function lme() of the package nlme to model grouped data that is auto-correlated in time and in space (the data was collected on different days via a moving monitor). I am aware that I can use the correlation classes corCAR1 and corExp (among other options) to model the temporal and spatial components of the auto-correlation. However, as far as I can tell, I can only

Let me modify my question about user-defined covariance structures for LME models: Can somebody tell me how I can see the code for the definition of the correlation structures that come with the NLME package. Specifically I like to see the code for the functions coef, corMatrix, and intialize for any of the pre-defined correlation structures, and use this as a template to define a new correlation

Hi: First my apologies for cross-posting. A few days back I posted my queries ar R-sig-geo but did not get any response. Hence this post. I am working on two parcel-level housing dataset to estimate the impact of various variables on home sale prices. I created the spatial weight metrics in ArcGIS 10 using sale year of four nearest houses to assign weights.  Next, I ran LM tests and then ran

Hi, I'm trying to introduce a (spatial) exponential correlation structure (with range=200 and nugget.effet of 0.3) in a lme model of this form: lme(ARBUS~YEAR, random=~1|IDSOUS). The structure of the data is "IDSOUS" "XMIN" "YMAX" "YEAR" "ARBUS" with 2 years of data and 5600 points for each year. I do:

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.

[I think I've seen this reported before but can't locate it any more. I believe this oddity (glitch? feature?) is behind a query that Jean-Baptiste Ferdy asked a year ago <http://finzi.psych.upenn.edu/Rhelp08/2008-October/176449.html>] It appears that glmmPQL looks in the global workspace, not within the data frame specified by the "data" argument, for the variables

Dear listers, I am trying to model the distribution of fox density over years in the Doubs department. Measurements have been taken on 470 plots in March each year and georeferenced. Average density is supposed to be different each year. In a first approach, I would like to use a general model of this type, taking spatial correlation into account:

Hi I am using R 1.7.1 on RH linux 9.0 > sum(unlist(lapply(ls(),function(x)object.size(get(x)))))/1024^2 [1] 2.424263 so I am not using much memory (I have a gig of ram on my machine) now in nlme > gtest<-groupedData(log(X8)~Time|sub,all[,c(names(all)[1:9],"X8")],outer=~A*B) > object.size(gtest)/1024 [1] 59.98438 > plot(gtest,outer=~Dose*chem,key=FALSE,asp=.5) Plotting

Following the Indomethicin example in Pinheiro & Bates, chapter 6, page 277 etc, coming to the following comand: fm2Indom.nlme <- update( fm1Indom.nlme, random = pdDiag(A1 + lrc1 + A2 ~ 1) ) debugging nlme gives the following output: Browse[1]> n debug: modelResid <- ~eval(model, data.frame(data, getParsNlme(plist, fmap, rmapRel, bmap, groups, beta, bvec, b, level,

Hello. I have 16 subjects with 1-4 obs per subject. I am using the package "nlme" to fit a simple random effects (variance components model) with 3 parameters: overall mean (fixed effect), between subject variance (random) and within subject variance (random). I need a 3x3 variance-covariance matrix that includes all 3 parameters in order to compute the variance of a specific

Hi, All: Does 'nlme' require the nonlinear function to be differentiable? I'm fitting structural change models to many related time series, and I'd like to improve the estimates of the change point times through some type of pooling. Unfortunately, I've so far been unable to get 'nlme' to work for me. The following toy example is the closest I've come

A few days ago I posted a question to this discussion group concerning to origin of an error message < subscript out of bounds > while using the nonlinear mixed model (nlme) function in R with a self-starting function. Thanks for those who responded. This posting is to explain what (I think) it causing the error. Pinheiro & Bates (2000, pages 342-347) describe how to construct a

Hi, I am not sure this is a bug but I can repeat it, The functions and data are below. I know this is nasty data, and it is very questionable whether a 4pl model is appropriate, but it is data fed to an automated tool and I would have hoped for an error. Does this repeat for anyone else? My details: > version _ platform i686-pc-linux-gnu

Dear all, I have been using the function gls in the package nlme in R to fit some spatial regressions (as described in Dormann et al.). However, I have been struggling trying to find a way to calculate a measure of effect size from these models, so I wanted to know if any of you had an idea on how to do this. More precisely, I am producing a multiple model with an exponential correlation

Dear all, I am trying to define my own corStruct which is different from the classical one available in nlme. The structure of this correlation is given below. I am wondering to know how to continue with this structure by using specific functions (corMatrix, getCovariate, Initialize,...) in order to get a structure like corAR1, corSymm which will be working for my data. Thanks in advance.

On 13/09/17 13:24, K?tia Emidio wrote: > Dear Rolf, > > Thanks for your help! > > What I need is a spatial window with shape equal to the figure attached. > This figure I made using ArcGis, but it is important to me make it in R. > After having this figure I will make some analysis using spatstat among > others. The points within figure are trees... (1) I believe that

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

Hi all, I am trying to set up a simple path analysis in the SEM package, but I am having some trouble. I keep getting the following error message or something similar with my model, and I'm not sure what I'm doing wrong: Error in solve.default(C) : system is computationally singular: reciprocal condition number = 2.2449e-20 In addition: Warning message: In sem.default(ram = ram, S = S,

Hi! I want to construct a path analysis model that can account for measurement error in totally aggregated parcels, which refer to parcels where all of the items in a scale are summed or averaged. If I am not mistaken, Bollen (1989) advocates the following formula for computing the error variance of each parcel: (1−α(parcel))×variance(parcel), such that α refers to Cronbach's alpha, which

Hi all, I am very knew to this R and Arcgis combination. I have already spent days attempting to find my answer but have failed, therefore could you please if possible explain any answers in simple form. I am trying to access Arcgis through R, following a few steps: 1) Set up working directory 2)connect to Python 3) import arcgis- however, when I get to the import stage it comes up with an error