[R] nlme & SASmixed in 2.0.1

Murray Jorgensen maj at stats.waikato.ac.nz
Wed Apr 6 01:40:43 CEST 2005


I assigned a class the first problem in Pinheiro & Bates, which uses the 
data set PBIB from the SASmixed package. I have recently downloaded 
2.0.1 and its associated packages. On trying

library(SASmixed)
data(PBIB)
library(nlme)
plot(PBIB)

I get a warning message
Warning message:
replacing previous import: coef in: namespaceImportFrom(self, 
asNamespace(ns))

after library(nlme) and a "pairs" type plot of PBIB.

Pressing on I get:

 > lme1 <- lme(response ~ Treatment, data=PBIB, random =~1| Block)
 > summary(lme1)
Error: No slot of name "rep" for this object of class "lme"
Error in deviance(object at rep, REML = REML) :
         Unable to find the argument "object" in selecting a method for 
function "deviance"
 >

Everything works fine under 1.8.1 and plot(PBIB) is of trellis style, 
which is what I think the authors intend.

Cheers,   Murray Jorgensen
-- 
Dr Murray Jorgensen      http://www.stats.waikato.ac.nz/Staff/maj.html
Department of Statistics, University of Waikato, Hamilton, New Zealand
Email: maj at waikato.ac.nz                                Fax 7 838 4155
Phone  +64 7 838 4773 wk    Home +64 7 825 0441    Mobile 021 1395 862




More information about the R-help mailing list