[R] Error Message from Variogram.lme Example

[Prof Brian Ripley]

On Tue, 14 Mar 2006, Prof Brian Ripley wrote:

> So this is most likely a bug in package nlme.  However, we need a
> reproducible example to be able to do anything about it, and without even
> the traceback() we cannot be sure that it is in nlme.
>
> Please follow the bug-reporting procedure.

It does seem to be a bug in Variogram.lme (assuming this is the BodyWeight 
dataset from package nlme).  The lines

         val <- na.omit(val)
         val$n.pairs <- as.vector(table(na.omit(cutDist)))

need to be interchanged since one bin has a zero count and gets omitted.

>
>
> On Mon, 13 Mar 2006, Rick Bilonick wrote:
>
>> When I try to run the example from Variogram with an lme object, I get
>> an error (although summary works):
>>
>> R : Copyright 2005, The R Foundation for Statistical Computing
>> Version 2.2.1  (2005-12-20 r36812)
>> ISBN 3-900051-07-0
>> ...
>>
>>>  fm1 <- lme(weight ~ Time * Diet, BodyWeight, ~ Time | Rat)
>> Error: couldn't find function "lme"
>>>      Variogram(fm1, form = ~ Time | Rat, nint = 10, robust = TRUE)
>> Error: couldn't find function "Variogram"
>>> library(nlme)
>>>  fm1 <- lme(weight ~ Time * Diet, BodyWeight, ~ Time | Rat)
>>>      Variogram(fm1, form = ~ Time | Rat, nint = 10, robust = TRUE)
>> Error in "$<-.data.frame"(`*tmp*`, "n.pairs", value = c(160, 0, 160,
>> 16,  :
>>        replacement has 10 rows, data has 9
>
> [...]
>
>>                Information on package 'nlme'
>>
>> Description:
>>
>> Package:       nlme
>> Version:       3.1-68.1
>> Date:          2006-01-05
>
>
>

-- 
Brian D. Ripley,                  ripley at stats.ox.ac.uk
Professor of Applied Statistics,  http://www.stats.ox.ac.uk/~ripley/
University of Oxford,             Tel:  +44 1865 272861 (self)
1 South Parks Road,                     +44 1865 272866 (PA)
Oxford OX1 3TG, UK                Fax:  +44 1865 272595