Skip to content

Mixed Effect Model - Jackknife error estimate

3 messages · Michael Falkowski, Kenneth Cabrera, Brian Ripley

Michael Falkowski
# 
Greetings,

I?ve fit the following mixed effects model using the NLME package:

hd.impute.lme <- lme(I(log(HEIGHT_M - 1.37)) ~ SPECIES + SPECIES:I(1/(DBH_CM + 2.54)),
        random = ~ I(1/(DBH_CM + 2.54)) | PLOTID,
        data = trees, na.action = na.exclude)

I would now like to extract a jackknife estimate of model error. I tried the following code, however, the estimate produced seems too low.

ss.ok <- 0
for (i in 1:dim(trees)[1])
 lme.msc <- ss.ok + sum(sd(lme(I(log(HEIGHT_M - 1.37)) ~ SPECIES + SPECIES:I(1/(DBH_CM + 2.54)),
        random = ~ I(1/(DBH_CM + 2.54)) | PLOTID,
        data = trees[-i,], na.action = na.exclude)$residuals))

rmse <- lme.msc/(dim(trees)[1]-1)

#output
#rmse = 0.4701232

Any suggestions?

Thanks,
Mike
Kenneth Cabrera
# 
Hello, R users:

I got the following problem when I try to install some packages using
the 2.1.0 patched version on Window$.

This is only an example.
The same happens when I try to install ade4.
(it download the packages and their dependency, but when it finish to 
download
it shows the same error)
What am I doing wrong?

Thank you for your help.


Kenneth
Brian Ripley
# 
Please look in the archives: in particular the *current* R-patched version 
does not do that and this has been discussed within the last week in a 
thread started by Erich Neuwirth.
On Tue, 10 May 2005, Kenneth Cabrera wrote: