similar to: error in nlme

Hi out there, I am trying to fit a species accumulation curve (increase in number of species known vs. sampling effort) for multiple regions and several bootstrap samples. The bootstrap samples represent different arrangements of the actual sample sequence. I fitted a series of nlme-models and everything seems OK, but since the observations are correlated I tried to include some correlation

Dear listers, First let me appologize if the same mail arrives multiple times. Recently I had some probelms sending my e-mails to the list. I encountered a problem when running glmmPQL procuedure doing multilevel modeling with a dichotomous outcome. Those are the two error messages I usually get: Error in logLik.reStruct(object, conLin) : NA/NaN/Inf in foreign function call (arg 3)

Hi, lme() users, Can some one tell me how to do this. I model Orthodont with the same G for random variables, but different R{i}'s for boys and girls, so that I can get sigma1_square_hat for boys and sigma2_square_hat for girls. The model is Y{i}=X{i}beta + Z{i}b + e{i} b ~ iid N(0,G) and e{i} ~ iid N(0,R{i}) i=1,2 orth.lme <- lme(distance ~ Sex * age, data=Orthodont, random=~age|Subject,

Hello, Excuse me for posting two questions in one day, but I figured it would be better to ask my questions in separate emails. I will again give the caveat that I'm not a statistician by training, but have a fairly decent understanding of probability and likelihood. As before, I'm trying to fit a nonlinear model to a dataset which has two main factors using nlme. Within the dataset

Hello, I am trying to design new variance structures by using fixed effects variables in combination with the VarPower function. That is, I would like to create and evaluate my own variance function in the data frame and then incorporate it into the model using varPower, with value=.5. As a start, I am trying to recreate the function of VarConstPower by introducing two new variables in the

Dear all, I am estimating a mixed-model in Ubuntu Raring (13.04¸ amd64), with the code: fm0 <- lme(rt ~ run + group * stim * cond, random=list( subj=pdSymm(~ 1 + run), subj=pdSymm(~ 0 + stim)), data=mydat1) When I check the approximate variance-covariance matrix, I get: > fm0$apVar [1] "Non-positive definite

Hi everyone. I'm working with nlme and I have a question regarding nlme fits that fail because of singularity issues. Specifically, there a way to return an nlme object when the estimation process runs into a singular matrix? For example, can the results up to the point of an error such as "Error in solve.default(pdMatrix(a, fact = TRUE)) : system is computationally singular" or

Hi I have found some postings referring to the fact that one can try and understand why a particular model is failing to solve/converge from the verbose output one can generate when fitting a nonlinear mixed model. I am trying to understand this output and have not been able to find out much: **Iteration 1 LME step: Loglik: -237.4517 , nlm iterations: 22 reStruct parameters: subjectno1

Hi All How can I fit a repeated measures analysis using gls? I want to start with a unstructured correlation structure, as if the the measures at the occations are not longitudinal (no AR) but plainly multivariate (corSymm). My data (ignore the prox_pup and gender, occ means occasion): > head(dta,12) teacher occ prox_self prox_pup gender 1 1 0 0.76 0.41 1 2

Greetings, 1) Is there a nice way of extracting the variance estimates from an lme fit? They don't seem to be part of the lme object. 2) In a series of simulations, I am finding that with ML fitting one of my random effect variances is sometimes being estimated as essentially zero with massive CI instead of the finite value it should have, whilst using REML I get the expected value. I guess

Hope one of you could help with the following question/problem: We would like to explain the spatial distribution of juvenile fish. We have 2135 records, from 75 vessels (code_tripnr) and 7 to 39 observations for each vessel, hence the random effect for code_tripnr. The offset (‘offsetter’) accounts for the haul duration and sub sampling factor. There are no extreme outliers in lat/lon. The model

Dear R users, I'm attempting to fit a GLM with random effects using the tweedie family for the error structure. I'm getting the error: iteration 1 Error in logLik.reStruct(object, conLin) : NA/NaN/Inf in foreign function call (arg 3) I'm running V2.1.0 I notice from searching the lists that the same error was reported in May 2004 by Spencer Graves, but no-one was able to

I'm attempting to analyze some survey data comparing multiple docks. I surveyed all of the slips within each dock, but as slips are nested within docks, getting multiple samples per slip, and don't really represent any meaningful gradient, slip is a random effect. There are also an unequal number of slips at each dock. I'm having syntactical issues, however. When I try

Dear Friends, I am writing to seek any help on "lme" error message. I am using lme to do Mixed-model linear regression. I use my simulated data. However, sometimes, I get the following error message, which I do not understand. "Error in solve.default(pdMatrix(a, fact=TRUE)): system is computationally sigular" I would appreciate any help about it. Thanks a lot Jingyuan Fu

Dear All, I was wondering whether any one knows of a matrix bending function in R that can turn non-positive definite matrices into the nearest positive definite matrix. I was hoping there would be something akin to John Henshall's flbend program (http://agbu.une.edu.au/~kmeyer/pdmatrix.html), which allows the standard errors of the estimated matrix elements to be considered in the

When trying to fit a generalized linear mixed model using glmmPQL: > fit0 <- glmmPQL(ifelse(response=="A",1,0)~gender,data=set1, random=~1|subject,family=binomial) iteration 1 Error in solve.default(pdMatrix(a, fact = TRUE)) : Lapack routine dgesv: system is exactly singular Could this be occuring because the paired responses for each subject are always the same? If

Hi, I'm having trouble with glmmPQL from the MASS package. I'm trying to fit a model with a binary response variable, two fixed and two random variables (nested), with a sample of about 200,000 data points. Unfortunately, I'm getting an error message that is difficult to understand without knowing the internals of the glmmPQL function. > model <- glmmPQL(primed ~

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

Hi, I am trying to run a multilevel model with time nested in people and people nested in dyads (3 levels of nesting) by initially running a series of models to test whether the slope/intercept should be fixed or random. The problem that I am experiencing appears to arise between the random intercept, fixed slope equation AND. (syntax: rint<-lme(BDIAFTER~BDI+WEEK+CORUMTO,

The print method for lme *prints out* the fitted correlation matrix for the random effects. Is there any way to get these values as an object in R? I have examined the components of the lme object (called "junk" in the example below) and the components of summary(junk) without finding these numbers. (How I did this: I dumped the entire lme object to a text file and then used egrep to