[Rd] Problems with nlme (PR#471)

Douglas Bates bates@stat.wisc.edu
07 Mar 2000 08:43:35 -0600


berwin@maths.uwa.edu.au writes:

> bossiaea:/opt/R$ R CMD check -c nlme 
> Checking package `nlme' ...
>  Massaging examples into `nlme-Ex.R' ...
>  Running examples in package `nlme' ...
>  ERROR
> 
> The end of nlme-Ex.Rout is:
> 
> > rm(list = ls(all = TRUE)); .Random.seed <- c(0,rep(7654,3))
> > ###--- >>> `anova.gls' <<<----- Compare Likelihoods of Fitted Objects
> > 
> > 	## alias	 help(anova.gls)
> > 
> > ##___ Examples ___:
> > 
> > library(nlme)
> > data(Ovary)
> > # AR(1) errors within each Mare
> > fm1 <- gls(follicles ~ sin(2*pi*Time) + cos(2*pi*Time), Ovary,
> +            correlation = corAR1(form = ~ 1 | Mare))
> Warning message: 
> NA/Inf replaced by maximum positive value 
> > anova(fm1)
> Denom. DF: 305 
>                    numDF  F-value p-value
> (Intercept)            1 354.7375  <.0001
> sin(2 * pi * Time)     1  18.5035  <.0001
> cos(2 * pi * Time)     1   1.6633  0.1981
> > # variance changes with a power of the absolute fitted values?
> > fm2 <- update(fm1, weights = varPower())
> Error in update.gls(fm1, weights = varPower()) : 
> 	subscript out of bounds
> Execution halted
> 
> 
> (R CMD check MASS is also failing on an example that uses update.gls.)
> 
> On my machine at home I run into troubles even faster when running `R
> CMD check nlme'.  One of the early examples chokes on the line:
> > data(Orthodont)
> 
> I had done a "clean install" on that machine, i.e., I first blew away
> the older installation of R.  It seems that the datasets Orthodont and
> Pixel have been distributed with the nls library in earlier R versions
> but not so in R 1.0.0.  But they are not distributed with the current
> version of nlme either.

Ah, it appears that I need to update the nlme library.  I probably
won't be able to do that until tomorrow at the earliest.
-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-
r-devel mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html
Send "info", "help", or "[un]subscribe"
(in the "body", not the subject !)  To: r-devel-request@stat.math.ethz.ch
_._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._