Displaying 7 results from an estimated 7 matches for "numiter".
Did you mean:
num_iter
2007 Nov 01
1
A question about lme object
I have a question about the lme function in R. My question is: After I got
the object from function lme, why the numIter value of the object is always
NULL? Following is my code:
jjww<-lme(y~x*zz,data=simul,random=~x|group,
control=lmeControl(returnObject=TRUE))
attributes(jjww)
jjww$numIter
the first 20 observation of data simul are:
> simul
y x
zz...
2007 Nov 01
0
Question about lme object
Hi,
I have a question about the lme function in R. My question is: After I got
the object from function lme, why the numIter value of the object is always
NULL? Following is my code:
jjww<-lme(y~x*zz,data=simul,random=~x|group,
control=lmeControl(returnObject=TRUE))
attributes(jjww)
jjww$numIter
the first 20 observation of data simul are:
> simul
y x
zz...
1999 Jul 01
1
lme
...nce (mm)")
points(age[Sex=="0"],distance[Sex=="0"],type="p",pch=0)
points(Age[Sex=="1"],distance[Sex=="1"],type="p",pch=5)
for(i in unique(Subject))
lines(Age[Subject==i],distance[Subject==i],type="l",lty=2)
The function numIter doesn't function
Orthodont.fit1$numIter
produces
NULL
Most distressing is the fact that
summary(Orthodont.fit2)
produces the same fixed effects values
(Intercept) 16.340625
Age 0.784375
Sex 1.032102
Age.Sex -0.0304830
as Pinheiro and Bates 1995 (lme and nlme) but everything...
2004 Jun 17
2
nlme graphics in a loop problem
Hi, I'm fitting mixed effects models using the lme function of the nlme
package. This involves using the various associated plot functions.
However, when I attempt to fit a number of models using an loop, whilst
the models work, the plot functions fail. As a trivial example, the
following works:
library(nlme)
graphics.off()
x<-c(1:10)
y<-c(1:4,7:12)
2006 Apr 23
1
lme: null deviance, deviance due to the random effects, residual deviance
...If yes how to do it ? A lme object provides the following:
> names(glm6)
[1] "modelStruct" "dims" "contrasts" "coefficients"
[5] "varFix" "sigma" "apVar" "logLik"
[9] "numIter" "groups" "call" "method"
[13] "fitted" "residuals" "fixDF" "family"
so no $null.deviance and $deviance elements as in glm objects...
I tried to find out an answer on R-help & Pi...
2004 Jul 26
0
Problem with a while loop embedded in a function.
...nlme.control.parms <- list("maxIter"=maxIter[OK], "returnObject"=returnObject[OK])
fm1 <- try(nlme(nlsList.Object,
fixed = Asym + R0 + lrc ~ 1,
random = Asym ~ 1,
control = nlme.control.parms), silent=TRUE)
ok <- is.numeric(fm1$numIter)
print(OK)
print(nlme.control.parms)
}
return(list(fm1, OK))
}
Test <- M1(maxIter=maxIter, returnObject=returnObject)
---------------------------------
[[alternative HTML version deleted]]
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