similar to: nlme and simulation

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

Hi, I would be grateful for help with a problem which is irritating me. I am quite sure that I am doing something stupid, but I can't see what it is. I am running R 1.7 on Windows 2000. The graphics device is the PC screen. The graphics from the nlme demonstration in Bates an Pinheiro's manual work just as advertised. The CO2 data and the Orthodont data dsiplay

Dear all, I tried to create a groupedData object, where the grouping factor is not ordered. Here ist the code: library(nlme) test<-groupedData(conc~Time|Subject,order.groups=F,data=as.data.frame(Theoph)) > getGroups(test) Levels: 6 < 7 < 8 < 11 < 3 < 2 < 4 < 9 < 12 < 10 < 1 < 5 I still get an ordered factor. As always thanks for your help peter

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

Dear Friends, I am seeking for any help on an error message in lme functions. I use mixed model to analyze a data with compound symmetric correlation structure. But I get an error message: "Error in corMatrix.corCompSymm(object) : NA/NaN/Inf in foreign function call (arg 1)". If I change the correlation structure to corAR1, then no error. I have no clue how to solve this problem.

Dear listserv members, In Croissant (2008) “Panel Data Econometrics in R: The plm Package” the authors seem to indicate that the nlme package for R cannot correctly handle unbalanced panel data: “Moreover, economic panel datasets often happen to be unbalanced (i.e., they have a different number of observations between groups), which case needs some adaptation to the methods and is not

Dear R users: I have some difficulties analizing data with mixed effects NLME and the last version of R. More concretely, I have a repeated measures design with a single group and 2 experimental factors (say A and B) and my interest is to compare additive and nonadditive models. suj rv A B 1 s1 4 a1 b1 2 s1 5 a1 b2 3 s1 7 a1 b3 4 s1 1 a2

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 have 16 subjects with 1-4 obs per subject. I need a 3x3 variance-covariance matrix that includes all 3 parameters in order to compute the variance of a specific linear

I asked this before and I am going to try again in more applied terms. I am trying to use R to extract variance components for a two-factor random effects model with both factors crossed. It would also be nice to generate some confidence intervals as well. For example, a data set using SAS Proc Mixed is below followed by the four variance component estimates and the respective confidence

Dear all, I have the following question. I have to fit the hierarchical model for the hypothesis concern the individual-level effects by controlling for the individual -level attributes and national-level contextual effects on individuals by using R. O have to obtain the estimates of the impact of the second-level (national: GDP per capita) effects on individuals ( in this instance the impact

Dear R-experts, I have a problem to formulate non linear mixed effects when more than one explanatory variable is present. I'm using R 1.1.1 under Linux. The version number of the nlme library is 3.1-7. "help(nlme)" says that among several possibilities it is possible to code fixed effects according to fixed = a1 + a2 + ... + an ~ X1 + X2 + ... + Xm where the left hand side

How does one implement a likelihood-ratio test, to test whether the variances of the random effects differ between two groups of subjects? Suppose your data consist of repeated measures on subjects belonging to two groups, say boys and girls, and you are fitting a linear mixed-effects model for the response as a function of time. The within-subject errors (residuals) have the same variance in

Dear all, I need to stack the columns of a matrix with a variable number of columns. I am aware of new<-c(a[,1],a[,2]) if I have a matrix 'a' with e.g. two columns. For a simulation study I need to stack the matrix with a variable number of columns. Is there a way to do this? Many thanks in advance Peter

Dear all, Last autumn there was some discussion on the list of the warning Warning message: fitted probabilities numerically 0 or 1 occurred in: (if (is.empty.model(mt)) glm.fit.null else glm.fit)(x = X, y = Y, when fitting binomial GLMs with many 0 and few 1. Parts of replies: "You should be able to tell which coefficients are infinite -- the coefficients and their standard errors will

I am trying to understand better an analysis mean RT in various conditions in a within subjects design with the overall mean RT / subject as one of the factors. LME seems to be the right way to do this. using something like m<- lme(rt~ a *b *subjectRT, random= ~1|subject) and then anova(m,type = "marginal"). My understanding is that lme is an easy interface for dummy coding

I was recently asked to perform a GLM analysis (the person comes from the JMP world) on a repeated measures design. I have found some things using aov but I cannot find anything with glm. In fact, multiple regressions in general with repeated measures seems to be poorly covered in documentation. I remember SPSS has separate commands to handle them. I have within variables SOA and ec and a

Hi R-users, I am using lme to fit a linear mixed model with the nlme package, does anyone know if it is possible to obtain standard error estimates of the variance components estimators and an adequate method to test the significance of the variance component? Thanks, Berta. [[alternative HTML version deleted]]

Hallo, I have a question that is more statistic related than about nlme and R functioning. I want to fit a complicated nonlinear model with nlme with several different measures of transpiration taken on each of the 220 trees grouped in 8 families. The unknown parameters of the model are three + their variances (and covariances). I want to estimate the variances among families of the

Dear all, I am studying the possibility of using the nlme package in R to analyse field trials of agricultural crops. I have a problem with the syntax for the modelling of variance covariance structures. I can model the within-group covariance structure using the correlation argument and the covariance structure between different random effects of the same grouping level using

Hi, I would like to fit a model for a factorial design that allows for unequal variances in all groups. If I am not mistaken, this can be done in lm by specifying weights. A function intended to specify weights for unequal variance structures is provided in the nlme library with the varIdent function. Is it apropriate to use these weights with lm? If not, is there another possibility to do