[R] nlme & SASmixed in 2.0.1

From: Murray Jorgensen <maj_at_stats.waikato.ac.nz>
Date: Wed 06 Apr 2005 - 09:40:43 EST


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@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@waikato.ac.nz                                Fax 7 838 4155
Phone  +64 7 838 4773 wk    Home +64 7 825 0441    Mobile 021 1395 862

______________________________________________
R-help@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html
Received on Wed Apr 06 09:45:40 2005

This archive was generated by hypermail 2.1.8 : Fri 03 Mar 2006 - 03:31:02 EST