similar to: segfault with correlation structures in nlme

I'm getting an error from nlme that has me stymied. I have a data set ,'mydata', with variables: AChE, Dose, sex, set, and mrid; 'set' and 'mrid' indicate two levels of nesting, with 'set' nested within 'mrid'. I want to fit the model: mod <- nlme(AChE ~ Cexp(Dose, A, B, m), data=mydata, fixed = A+B+M~sex, random=A+B+m~sex | mrid/set,

Greetings R-helpers, I am attempting to fit an lme() while specifying a correlation structure, but I'm getting into trouble long before I get to that point. I am receiving the error: Error in y[revOrder] - Fitted : non-conformable arrays It doesn't seem to matter how simple or complex the model I specify is, it always gives this same error message. This makes me suspect something is

Dear All -- I am trying to use within a little table producing code an anova comparison of two gls fitted objects, contained in a list of such object, obtained using nlme function gls. The anova procedure fails to locate the second of the objects. The following code, borrowed from the help page of anova.gls, exemplifies: --------------- start example code --------------- library(nlme) ##

Hi, Well after trying to fix the assign commands in version 3.1.7 of your nlme package I still cannot get the plot functions to work. Basically I have no other packages installed (apart from those provided with R) so nlme is the first. my R version info is; platform sparc-sun-solaris2.7 arch sparc os solaris2.7 system sparc, solaris2.7 status

Hi there, I try to compare the mixed model package "lme" by Splus and R. I used the dataset "Ovary" and the following code assuming AR(1) model for the error term: lme(follicles ~ sin(2*pi*Time) + cos(2*pi*Time), data=Ovary, random = pdDiag(~sin(2*pi*Time) ) , correlation=corAR1() ) But I got different results! And then I used a simpler model: lme(follicles ~

Hi, All: Is there a way to get random effects for ARMA parameters? Consider the following example from the 'corARMA' help page: fm1Ovar.lme <- lme(follicles ~ sin(2*pi*Time) + cos(2*pi*Time), data = Ovary, random = pdDiag(~sin(2*pi*Time))) fm5Ovar.lme <- update(fm1Ovar.lme, corr = corARMA(p = 1, q = 1)) fm5Ovar.lme Linear

There has been some recent discussion on this list about the value of using intervals with lme() to check for whether a model is ill-defined. My question is, what else can drive very large confidence intervals for the variance components (or cause the error message "Error in intervals.lme(Object) : Cannot get confidence intervals on var-cov components: Non-positive definite approximate

Hi there fellow R-users, I am trying to follow an example of modelling a serial correlation structure in the textbook "Mixed Effects Model in S and Splus". However, I am getting some very odd results. Here is what I am trying to run: library(nlme) data(Ovary) fm1<-lme(follicles~sin(2*pi*Time)+cos(2*pi*Time),data=Ovary,random=pdDiag(~s in(2*pi*Time))) ### The example is fine up

Dear R developers, first of all let me join the chorus of congratulations for the release of R 1.0.0. Well, done! Unfortunately, I find it necessary to e-mail in a bug report regarding the `nlme' package. On my office machine I experience the following trouble: bossiaea:/opt/R$ R CMD check -c nlme Checking package `nlme' ... Massaging examples into `nlme-Ex.R' ... Running

Trying following example from Pinheiro and Bates in order to fit an ARMA(1,1) model: library(nlme) fm1Ovary.lme<-lme(follicles~sin(2*pi*Time)+cos(*pi*Time),data=Ovary,random=p dDiag(~sin(2*pi*Time))) fm5Ovary.lme<-update(fm1Ovary.lme,corr=corARMA(p=1,q=1)) I get follwing error message: Error in "coef<-.corARMA"(`*tmp*`, value = c(62.3428455941166, 62.3428517930051 :

Hi all ! I'm new in this list and newbie about R I'm trying to use R scripts (as in the attached file) for creating some distributions plots of data retrieved by a workflow(with Rserve, to be precise). I was able to do it (even if not in a beatiful way, I have to improve it especially about labels and coordinates) with number inputs like :

Dear Users, >From previous analysis (semi-variograms using package gstat), I found spatial autocorrelation in my dataset. The best fitted model to this spatial correlation structure is the Gaussian model (Spherical, Exponential, Linear tested and comparison done by Sum of Square errors). So I used corGaus function to define this spatial autocorrelation in my lme model using the option

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 am having trouble fitting correlation structures within nlme. I would like to fit corCAR1, corGaus and corExp correlation structures to my data. I either get the error "step halving reduced below minimum in pnls step" or alternatively R crashes. My dataset is similar to the CO2 example in the nlme package. The one major difference is that in my case the 'conc' steps are

Dear R-helpers, I'd like to use glmmPQL to predict binary responses based on a data.frame data1 containing N entries (N<1000): target covariate1 covariate2 covariate3 ... covariateM cluster 134131 1 -0.30031885 0 0 -2.886870e-07 1 38370 1 -0.04883229 0 1 -1.105720e-07 1 19315 1 -0.11084267

Dear reader, I am trying to simulate 1000 data from nonlinear model in order to be able to do mixed effect analysis. If the program works but give you following warning message, what should I do? Can I still accept the result, which is about the precision of model parameter estimation? FALSE CONVERGENCE. in: ms( ~ - logLik(nlmeSt, nlmePars), start = list(nlmePars = c(coef(nlmeSt))), control

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,

Colleagues, I am using "nls" successfully (2.11.1, OS X) but I am having difficulties retrieving part of the output - residual sum of squares. I have assigned the output to FIT: > > FIT > Nonlinear regression model > model: NEWY ~ PMESOR + PAMPLITUDE * cos(2 * pi * (NEWX - POFFSET)/PERIOD) > data: parent.frame() > PMESOR PAMPLITUDE POFFSET >

Hi, I am trying to plot the original data with the line of the model using the predict function. I want to add SE to the graph, but not sure how to get them out as the predict function for gls does not appear to allow for SE=TRUE argument. Here is my code so far: f1<-formula(MaxNASC40_50~hu3+flcmax+TidalFlag) vf1Exp<-varExp(form=~hu3) B1D<-gls(f1,correlation=corGaus(form=Lat~Lon,

Dear All: I updated my R program as well as associated packages yesterday. Currently my R version is 2.2.1 running under WINXP SP-2. When I tried to list (summary) an nlme object that I developed before, I got the following error message: [ Error in .C("ARMA_constCoef", as.integer(attr(object, "p")), as.integer(attr(object, : C entry point "ARMA_constCoef"