Displaying 3 results from an estimated 3 matches for "fm1rail".
2006 May 08
1
Repeatability and lme
...component (Lessells & Boag, 1987). Usually,
repeatability is used together with one-way ANOVA, but I'd rather stick to
lme. So for the introductionary example to lme in Pinheiro and Bates (2001)
with Rails, I thought it would be appropriate to follow this procedure:
> library(nlme)
> fm1Rail.lme <- lme(travel ~ 1, data = Rail, random = ~ 1 | Rail)
> summary(fm1Rail.lme)
...
Random effects:
Formula: ~1 | Rail
(Intercept) Residual
StdDev: 24.80547 4.020779
...
> as.numeric(VarCorr(fm1Rail.lme)[1]) / (as.numeric(VarCorr(fm1Rail.lme)[1])
+ as.numeric(VarCorr(fm1Rai...
2008 Jan 03
2
confidence interval too small in nlme?
...library('nlme')
M <- 6
n <- 3
beta <- 67
sigma.b <- 25
sigma <- 4
Rail <- rep(1:M, each=n)
set.seed(56820)
B <- 10000
num.wrong <- 0
error.fraction <- Ks <- c()
for (K in 1:B) {
travel <- beta + rep(rnorm(M, sd=sigma.b), each=n) + rnorm(M*n, sd=sigma)
fm1Rail.lme <- lme(travel ~ 1, random = ~ 1 | Rail)
CI <- intervals(fm1Rail.lme, which='fixed')$fixed
if ((CI[1, 'lower'] > beta) || (CI[1, 'upper'] < beta))
num.wrong <- num.wrong + 1
if (K %% 200 == 0) {
error.fraction <- c(error.fraction, num.wrong...
2004 Nov 30
1
lme in R-2.0.0: Problem with lmeControl
...Verbose=TRUE))
but every time i change the defualt values maxIter (e.g. maxIter=1,
niterEM=0) on ones specified by me, the call returns all the iterations
needed until it's converged.
and this is exactly the problem i will to get round. ( e.g. in
example on p.81 of Pinheiro/Bates,2000:
fm1Rail.lme<-lme(...,control=list(maxIter=1,...)) )
so i have tried with option msMaxIter=... and this works.
The other problem is, that i even can not see (in R !!!) the output from
iterations, despite the msVerbose=TRUE specification and setting
options(verbose=TRUE) (The S-plus can do...