An embedded and charset-unspecified text was scrubbed... Name: not available URL: <https://stat.ethz.ch/pipermail/r-help/attachments/20120418/98f81b48/attachment.pl>
how to plot separate lm ablines on the same xyplot by group
3 messages · Louis Plough, Bert Gunter, Richard M. Heiberger
Please read ?panel.superpose again and pay attention to the panel.groups argument, where this can be specified. -- Bert
On Wed, Apr 18, 2012 at 12:34 PM, Louis Plough <lplough at usc.edu> wrote:
Hi,
I am trying to use xyplot to plot the relationship between size and day
(y~x) by a food factor that has two levels, low and high. I have 3 reps per
factor/day. ?I want the plots from each food treatment on the same axiss,
so I used this code:
xyplot(Size ~ Day, groups = Food, data = louis.data.means,col=1,
pch=c(1,17),
panel=function(x,y,groups,...){
panel.superpose(x,y,groups,...)
tmp.lm<-lm(y~x)
panel.abline(tmp.lm)
panel.text(2, 250, label=format(tmp.lm$coefficients[2], digits=4), pos=4)
}
)
This produces a graph of the two treatments (open circles for the low food
?vs triangles for the high food) on the same plot, but only one regression
line (and slope) ?which seems to splits the difference between the two
factors (treats them as the same data set). ?I would like to produce a
separate regression line for the data from each of the two factors (high
food vs low food).
Is there a way to subset the lm by the factor "food"?
Louis
? ? ? ?[[alternative HTML version deleted]]
______________________________________________ R-help at r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code.
Bert Gunter Genentech Nonclinical Biostatistics Internal Contact Info: Phone: 467-7374 Website: http://pharmadevelopment.roche.com/index/pdb/pdb-functional-groups/pdb-biostatistics/pdb-ncb-home.htm
An embedded and charset-unspecified text was scrubbed... Name: not available URL: <https://stat.ethz.ch/pipermail/r-help/attachments/20120418/f31a11af/attachment.pl>