search for: lmepars

Dear R community, I'm fitting a complex mixed-effects model that requires numerous iterations and function evaluations. I note that nlminb accepts a list of control parameters, including eval.max. Is there a way to change the default eval.max value for nlminb when it is being called from lme? Thanks for any thoughts, Andrew -- Andrew Robinson Department of Mathematics and Statistics

...ied higher settings > gamm3<-gamm(count~offset(offsetter)+s(lon,lat),random=list(code_tripnr=~1),family="poisson", niterPQL=200, control=lmeControl(opt="optim")) Maximum number of PQL iterations: 200 iteration 1 iteration 2 Error in optim(c(coef(lmeSt)), function(lmePars) -logLik(lmeSt, lmePars), : initial value in 'vmmin' is not finite gamm3<-gamm(count~offset(offsetter)+s(lon,lat),random=list(code_tripnr=~1),family="poisson", niterPQL=200,control=lmeControl(minAbsParApV ar=0.0000000000001)) Maximum number of PQL iterations: 200 i...

...ample size is about 2965 and 3 factors: year (5 levels), ssize (4 levels), condition (2 levels). When I issue the following command: > lme(var~year*ssize*condition,random=~ssize+condition|subject,data=smp,method ="ML") I got the following error: Error in logLik.lmeStructInt(lmeSt, lmePars) : Calloc could not allocate (65230 of 8) memory In addition: Warning message: Reached total allocation of 120Mb: see help(memory.size) I'm currently using a Win2000 machine with 128Mb RAM and a 1.2 Gb processor. My version of R is 1.7.1. Thanks in advance, Rodrigo Abt. Department of...

...[1] "site" "time" "shade" "treatp" "aphid" > library(nlme) > aphids<-groupedData(aphid~time|site/shade/treatp,data=aphids,outer=~shade*treatp) > model<-lme(aphid~shade*treatp,random=~time|site) Error in logLik.lmeStructInt(lmeSt, lmePars) : Calloc could not allocate (89359209 of 8) memory In addition: Warning message: Reached total allocation of 254Mb: see help(memory.size) *********************************************************** My response variable is number of aphids/colony on shaded cocoa trees, measured over 15 da...

Dear R-users, I did some more research and I'm still not sure how to set up an ANCOVA with nestedness. Specifically I'm not sure how to express chicks nested within boxes. I will be getting Pinheiro & Bates (Mixed Effects Models in S and S-Plus) but it will not arrive for another two weeks from our interlibrary loan. The goal is to determine if there are urbanization (purban)

Hi All, I've been having a little trouble with creating a loop that will run a a series of t.tests for inspection, Below is the code i've tried, and some checks i've looked at. I've used the get(paste()) idea as i was told previously that the use of the eval should try and be avoided. I've run a single syntax to check that my systax is correct and works without any problems