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)