Displaying 6 results from an estimated 6 matches for "lmepar".
Did you mean:
lmepars
2006 Jul 23
1
How to pass eval.max from lme() to nlminb?
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
2009 Jan 22
1
convergence problem gamm / lme
...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...
2003 Nov 10
8
Memory issues..
...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 o...
2005 Dec 14
3
Memory shortage running Repeated Measures (nlme)
...[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 d...
2006 Jan 24
4
nested ANCOVA: still confused
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)
2009 Jan 28
2
t.test in a loop
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