Several densityplots in single figure
On Wed, Feb 22, 2012 at 8:49 AM, David Winsemius <dwinsemius at comcast.net> wrote:
After going back and constructing a proper dataset, you should be passing 'groups' into the panel function and picking it up inside panel.abline.
Close, but unfortunately things get more complicated when using groups
in densityplot. A straight up panel function wouldn't work for him out
of the box (and you/list be getting more follow ups). He'll need
something like:
densityplot(~value | g1, groups=g2, data= combinedlongformat,
panel=panel.superpose, panel.groups=function(x,...,group.number){
panel.abline(v=mean(x),lwd=(1:0)[group.number]) # seems like he only wants x1
panel.densityplot(x,...)
})
At this stage in the game
trellis.focus('panel',1,1) ; panel.abline(v=calculatedmean1)
trellis.focus('panel',2,1) ; panel.abline(v=calculatedmean2)
...
trellis.unfocus()
Would probably go down easier...
Just a thought
Cheers,
Elai
You
can also recover the 'panel.number()' and then use them both to recover the means from the global environment when inside panel.abline. I probably would have used tapply to create a table with an index (in the enclosing enviroment) rather than creating 4 separate values to choose That way you could use one object name but construct the 4 indices. But for that to work you would have needed to follow my example about binding the 4 datasets into one. One way to find previous worked examples in the rhelp Archives is to search with your favorite engine on strategies like "panel.number groups Sarkar" or "panel.number groups Andrews" Deepayan Sarkar and Felix Andrews are the two persons from whom I have learned the most regarding the fine points of lattice plots. And will you please learn to post in plain text? -- david
? ? ? ? ? ? ? ?}
)
On 22 February 2012 13:48, David Winsemius <dwinsemius at comcast.net> wrote:
On Feb 22, 2012, at 5:28 AM, josh rosen wrote:
Hi,
I have created two separate overlapping density plots- see example code
below.
What I wish now to do is combine them into one figure where they sit side
by side.
Any help would be great!
many thanks in advance, josh.
#####################
thedataA <- data.frame(x1=rnorm(100,1,1),x2=rnorm(100,3,1)) #create data
thedataA.m<-melt(thedataA)
densityplot(~value, thedataA.m, groups=variable,auto.key=list(columns=2),
? ? panel = function(x, y, ...) {
? ? ? ? ? ? panel.densityplot(x, ...)
? ? ? ? ? ? panel.abline(v=0)
? ? }
)
The syntax for grouping (which gives theoverlaid but different colors as
default output) and "separation" is fairly simple. Use the "|" operator for
separated plots and the " .., groups= <var>, .." ?parameter for overlaying
results. It's only going to work easily if they are all in the same dataset.
Try:
?bigset <- cbind( rbind(thedataA.m, thedataB.m), ABgrp=rep(c("datA",
"datB"), each=200) )
?densityplot(~value|ABgrp, data=bigset, groups=variable,
auto.key=list(columns=2),
? ? ?panel = function(x, y, ...) {
? ? ? ? ? ? ?panel.densityplot(x, ...)
? ? ? ? ? ? ?panel.abline(v=0) ? } )
And please work on whatever practice is producing duplicate postings.
--
David.
thedataB <- data.frame(x1=rnorm(100,2,1),x2=rnorm(100,4,1)) #create data
thedataB.m<-melt(thedataA)
I assume that is a copy-paste-fail-to-correct error.
densityplot(~value, thedataB.m, groups=variable,auto.key=list(columns=2),
? ? panel = function(x, y, ...) {
? ? ? ? ? ? panel.densityplot(x, ...)
? ? ? ? ? ? panel.abline(v=0)
? ? }
)
######################
? ? ? [[alternative HTML version deleted]]
David Winsemius, MD West Hartford, CT
______________________________________________ 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.