R help - May 2005 - Mixed Effect Model - Jackknife error estimate

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

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)
> > utils:::menuInstallPkgs()
> trying URL 
>
'http://cran.at.r-project.org/bin/windows/contrib/2.1/CoCo_0.1.6.5.zip'
> Content type 'application/zip' length 1895592 bytes
> opened URL
> downloaded 1851Kb
>
> package 'CoCo' successfully unpacked and MD5 sums checked
> bundle 'CoCo' successfully unpacked and MD5 sums checked
> Error in sprintf(gettext("unable to move temp installation
'%d' to
> '%s'"),  :
>         use format %s for character objects

What am I doing wrong?

Thank you for your help.


Kenneth