From: Sébastien <pomchip_at_free.fr>

Date: Wed, 20 Jun 2007 23:57:13 -0400

etc...

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 and provide commented, minimal, self-contained, reproducible code. Received on Thu 21 Jun 2007 - 04:02:43 GMT

Date: Wed, 20 Jun 2007 23:57:13 -0400

Plot Model Individuals Time Observed Predicted

1 1 A 1 0.05 10 10.2 2 1 A 1 0.10 20 19.5 etc... 10 1 B 1 0.05 10 9.8 11 1 B 1 0.10 20 20.2

etc...

There are p "levels" in mydata$Plot, m in mydata$Model, n in mydata$Individuals and t in mydata$Time (Note that I probably use the word levels improperly as all columns are not factors). Basically, this dataset summarizes the t measurements obtained in n individuals as well as the predicted values from m different modeling approaches (applied to all individuals). Therefore, the observations are repeated m times in the Observed columns, while the predictions appears only once for a given model an a given individual.

What I want to write is a R batch file creating a Trellis graph, where each panel corresponds to one individual and contains the observations (as scatterplot) plus the predicted values for all models (as lines of different colors)... $Plot is just a token: it might be used to not overload graphs in case there are too many tested models. The fun part is that the values of p, m, n and t might vary from one dataset to the other, so everything has to be coded dynamically.

For the plotting part I was thinking about having a loop in my code containing something like that:

for (i in 1:nlevels(mydata$Model)) {

subdata<-subset(mydata,mydata$Model=level(mydata$Model)[i]) xyplot(subset(Observed + Predicted ~ Time | Individuals, data = subdata) #plus additionnal formatting code

}

Unfortunately, this code simply creates a new Trellis plot instead of adding the model one by one on the panels. Any idea or link to a useful command will wellcome.

Sebastien

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 and provide commented, minimal, self-contained, reproducible code. Received on Thu 21 Jun 2007 - 04:02:43 GMT

Archive maintained by Robert King, hosted by
the discipline of
statistics at the
University of Newcastle,
Australia.

Archive generated by hypermail 2.2.0, at Thu 21 Jun 2007 - 13:31:58 GMT.

*
Mailing list information is available at https://stat.ethz.ch/mailman/listinfo/r-help.
Please read the posting
guide before posting to the list.
*