Easy to write - here's a quick piece of code (without bells and whistles
though...)
dat1 <- rnorm(1000)
dat2 <- rgamma(3000, shape=3)
biHist <- function(d1, d2) {
col1 <- "#FFCCBB"
col2 <- "#BBCCFF"
hAll <- hist(c(d1, d2))
hD1 <- hist(d1, breaks=hAll$breaks)
hD2 <- hist(d2, breaks=hAll$breaks)
plot(NULL,
xlim = c(-1.05 * max(hD1$counts), 1.05 * max(hD2$counts)),
ylim = c(min(hAll$breaks), max(hAll$breaks)),
xlab = "frequencies",
ylab = "values",
xaxt = "n")
axis(side = 1,
lwd = 0,
lwd.ticks = 1,
at = axTicks(1),
labels = abs(axTicks(1)))
abline(v = 0)
rect(-hD1$counts,
hAll$breaks[-length(hAll$breaks)],
0,
hAll$breaks[-1],
col=col1)
rect(0,
hAll$breaks[-length(hAll$breaks)],
hD2$counts,
hAll$breaks[-1],
col=col2)
}
biHist(dat1, dat2)
Let me know in case this needs modifications that you can't easily make
yourself
:-)
B.
On Jun 2, 2015, at 4:43 PM, David Winsemius <dwinsemius at comcast.net>
wrote:
>
> On Jun 2, 2015, at 11:13 AM, Dan Abner wrote:
>
>> Hi all,
>>
>> Does anyone know how to create a graph of propensity scores like the
>> one on the left in the attachment? I can easily generate the one on
>> the right: How does one force the respective histograms to share the
>> same rotated x-axis? Is it possible to set the horizontal white space
>> between these to graphs to 0 so that there is only the single line
>> running down the middle of the graph?
>>
>
> I think you will find that "pyramid" plots of this sort are in
some packages and that similar requests (to get rid of the middle space) have
been made in the past. Pyramid plots are often used by demographers to compare
and display age distributions of both genders. Why not see if that work can be
found with a search engine. I use this MarkMail link on my browser to search the
Archive
>
> http://markmail.org/search/?q=list%3Aorg.r-project.r-help
>
> ....(but you could also direct the Google Advanced Search function to the
real Archives if you wanted:
>
> https://stat.ethz.ch/pipermail/r-help/
>
>
>> Thanks,
>>
>> Dan
>> <graphs.pdf>______________________________________________
>> R-help at r-project.org mailing list -- To UNSUBSCRIBE and more, see
>> 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.
>
> David Winsemius
> Alameda, CA, USA
>
> ______________________________________________
> R-help at r-project.org mailing list -- To UNSUBSCRIBE and more, see
> 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.