[R] xyplot: discrete points + continuous curve per panel
Petr Pikal
petr.pikal at precheza.cz
Thu Dec 14 10:17:42 CET 2006
Hi
there is probably better solution but you can try to fiidle with this
idea, which adds stight lines to each panel one after another.
# based on Gabor Grothendieck's code suggestion
# adds straight lines to panels in lattice plots
addLine<- function(...) {
tcL <- trellis.currentLayout()
for(i in 1:nrow(tcL))
for(j in 1:ncol(tcL))
if (tcL[i,j] > 0) {
trellis.focus("panel", j, i, highlight = FALSE)
panel.abline(...)
trellis.unfocus()
}
}
You need to change panel.abline(...) part maybe to panel.curve or
panel.segments or?
HTH
Petr
On 13 Dec 2006 at 23:22, RMan54 wrote:
Date sent: Wed, 13 Dec 2006 23:22:41 -0800 (PST)
From: RMan54 <RMan54 at cox.net>
To: r-help at stat.math.ethz.ch
Subject: [R] xyplot: discrete points + continuous curve per panel
>
> I have a number of x, y observations (Time, Conc) for a number of
> Subjects (with subject number Subj) and Doses. I can plot the
> individual points with xyplot fine:
>
> xyplot(Conc ~ Time | Subj,
> Groups=Dose,
> data=myData,
> panel = function(x,y) {
> panel.xyplot(x, y)
> panel.superpose(???) # Needs more here
> }
> )
>
> I also like to plot on each panel (there is one Subj per panel) a
> continuous curve with predictions that I can calculate from a rather
> complicated function:
>
> myPred <- (time, subj, dose) {
> returns predicted value of Conc for a given time, subj and dose
> }
>
> The predicted curves are different for each panel.
>
> How do I plot the predictions? I have tried to add panel.superinpose
> in the xyplot portion but can't link to the myPred function. I also
> know about panel.curve but couldn't make it work.
>
> My attempt is to calculate the predictions on the fly. Is this
> possible? Or do I need to calculate all predictions first and put the
> results in a data frame. Thanks for any help, Rene -- View this
> message in context:
> http://www.nabble.com/xyplot%3A-discrete-points-%2B-continuous-curve-p
> er-panel-tf2818931.html#a7867892 Sent from the R help mailing list
> archive at Nabble.com.
>
> ______________________________________________
> 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.
Petr Pikal
petr.pikal at precheza.cz
More information about the R-help
mailing list