[R] error in plot.lmList

[Douglas Bates]

Arne.Muller at sanofi-aventis.com wrote:
> Hello,
> 
> in R-2.1.0 I'm trying to prodice trellis plots from an lmList object as described in the help for plot.lmList. I can generate the plots from the help, but on my own data plotting fails with an error message that I cannot interpret (please see below). Any hints are greatly appreciapted.
> 
> 	kind regards,
> 
> 	Arne
> 
> 
>>dim(d)
> 
> [1] 575   4
> 
>>d[1:3,]
> 
>   Level_of_Expression SSPos1 SSPos19 Method
> 1                11.9      G       A   bDNA
> 2                24.7      T       T   bDNA
> 3                 9.8      C       T   bDNA
> 
>>fm <- lmList(Level_of_Expression ~ SSPos1 + SSPos19 | Method, data=d)
>>fm
> 
> Call:
>   Model: Level_of_Expression ~ SSPos1 + SSPos19 | Method 
>    Data: d 
> 
> Coefficients:
>            (Intercept)   SSPos1C    SSPos1G   SSPos1T SSPos19C SSPos19G   SSPos19T
> bDNA          25.75211 -6.379701  -9.193304 10.371056 24.32171 24.06107  9.7357724
> Luciferase    23.79947  4.905679  -7.747861  8.112779 48.95151 48.15064 -0.2646783
> RT-PCR        56.08985 -7.352206 -15.896556 -2.712313 19.91967 24.28425 -2.2317071
> Western       14.03876  2.777038 -14.113157 -7.804959 24.62684 25.50382  8.3864782
> 
> Degrees of freedom: 575 total; 547 residual
> Residual standard error: 25.39981
> 
>>plot(fm, Level_of_Expression ~ fitted(.))
> 
> Error in plot.lmList(fm, Level_of_Expression ~ fitted(.)) : 
>         Object "cF" not found
> 
> what is object cF ...?
> 
> ______________________________________________
> R-help at 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

Could you provide us with a traceback or, perhaps better, a copy of the
data frame d (off-list and disguised, if necessary) so we can reproduce
the problem?

P.S. It may be easier to use xyplot directly as in

 xyplot(Level_of_Expression ~ fitted(fm), d)