An embedded and charset-unspecified text was scrubbed... Name: not available Url: https://stat.ethz.ch/pipermail/r-help/attachments/20051215/934c6a71/attachment.pl
Lattice graphics with combined plot types
2 messages · Shawn Way, Deepayan Sarkar
On 12/15/05, Shawn Way <Shawn.Way at biogenidec.com> wrote:
The data is of two forms, ie one numeric and another ordinal ie 1.5 and
<.5. The X Value is a factor.
I'm trying to create a lattice conditioned plot with the following
characteristics:
1) The plot is conditioned using the form (Conductivity~Day|Valve)
2) The plot should use a barplot for the ordinal (<.5) and dots for the
numeric (1.5)
3) A line should be created specifying a limit (ie 0.5)
What I would like is the bar on Day 3 to be a dot at 0.5 (a red one would
be even better)
I have the following working:
data <-
data.frame(Conductivity=c(.24,.24,.5),Day=as.factor(c(1,2,3)),Valve=as.factor(c("60-V-234","60-V-234","60-V-234")))
print(barchart(Conductivity~Day|Valve,
data=data,
horizontial=FALSE,
# type="p",
xlab=list(label="Day",cex=3),
ylab=list(label=expression(paste("Endotoxin (EU/mL)")),cex=3),
ylim=c(0,1),
# index.cond=list(c(2,1)),
scales=list(x=list(rot=90,cex=2),y=list(cex=2))
,panel=function (x,y,...) {
panel.barchart(x,y,...)
panel.abline(h=0.25,col="red")
}
# main="Stable Pressure for 70-DS-005 Protocol 02-634-02B"
)
)
You seem to have used a cutoff of 0.25 (rather than 0.5) in your example.
Any thoughts? I know it's something simple, I just cannot see it...
How about
panel=function (x, y, ...) {
below <- y < 0.25
if (any(below))
panel.barchart(x[below], y[below], ...)
if (any(!below))
panel.xyplot(x[!below], y[!below], ...,
col = 'red', cex = 2, pch = 16)
panel.abline(h=0.25,col="red")
}
?
Deepayan
--
http://www.stat.wisc.edu/~deepayan/