Re: [R] nlme & SASmixed in 2.0.1

From: Deepayan Sarkar <deepayan_at_stat.wisc.edu>
Date: Wed 06 Apr 2005 - 10:19:22 EST

On Tuesday 05 April 2005 18:40, Murray Jorgensen wrote:
> 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.

SASmixed now depends on lme4, which conflicted (until recently) with nlme. If you had your calls to library(SASmixed) and library(nlme) reversed, you probably would have gotten a warning.

I think the simplest thing you can do is not to load SASmixed at all. Instead, use

data(PBIB, package = "SASmixed")

However, these datasets are (for all practical purposes) vanilla data frames, and you won't get trellis-style plots unless you create nlme groupedData objects yourself. (You should get them if you load the latticeExtra package manually, and then use 'gplot' instead of 'plot' to plot them.)

Deepayan

> 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



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 10:21:14 2005

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