[R] Trellis plots with two regression lines

ONKELINX, Thierry Thierry.ONKELINX at inbo.be
Tue Mar 11 13:56:13 CET 2008


Yvonne,

With ggplot2 you could tackle this problem with the code below.

library(ggplot2)
graphs<- read.table("table.txt", header=T)
ggplot(data = graphs, aes(x = bif, y = abund)) + geom_point() +
geom_smooth(method = "lm") + facet_grid(. ~ experiment)

Have a look at the ggplot2 website for more info
(http://had.co.nz/ggplot2/)

HTH,

Thierry

------------------------------------------------------------------------
----
ir. Thierry Onkelinx
Instituut voor natuur- en bosonderzoek / Research Institute for Nature
and Forest
Cel biometrie, methodologie en kwaliteitszorg / Section biometrics,
methodology and quality assurance
Gaverstraat 4
9500 Geraardsbergen
Belgium 
tel. + 32 54/436 185
Thierry.Onkelinx op inbo.be 
www.inbo.be 

Do not put your faith in what statistics say until you have carefully
considered what they do not say.  ~William W. Watt
A statistical analysis, properly conducted, is a delicate dissection of
uncertainties, a surgery of suppositions. ~M.J.Moroney

-----Oorspronkelijk bericht-----
Van: r-help-bounces op r-project.org [mailto:r-help-bounces op r-project.org]
Namens yvo
Verzonden: dinsdag 11 maart 2008 13:18
Aan: r-help op r-project.org
Onderwerp: [R] Trellis plots with two regression lines


http://www.nabble.com/file/p15976467/table.pdf table.pdf 

I want to create 8 trellis (|experiment) scatterplots with regression
lines
for prey and pred in each plot with different colouring. So far i am
creating trellis plots- but the separation in these groups is never
possible.

library(lattice)
graphs<- read.table("table.txt", header=T)
attach(graphs)
names(graphs)
str(graphs)  
xyplot(abund ~ bif | experiment, data=graphs, groups=stage, 
        panel=function(x,y){panel.xyplot(x,y, pch=16)
panel.abline(lm(y~x))})

Thanks a lot for your help!!
Yvonne Fabian
University of Fribourg


-- 
View this message in context:
http://www.nabble.com/Trellis-plots-with-two-regression-lines-tp15976467
p15976467.html
Sent from the R help mailing list archive at Nabble.com.

______________________________________________
R-help op 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.



More information about the R-help mailing list