[R] Combined plot: Scatter + density plot

[Jannis]

Sorry Ralf,

I was in a hurry when I replied to your plot. The function you are looking for is:

dx=density(x)
plot(dx)

That works with the code you send in your mail. Just adjust the plot limits and change x and y for the density plot of y and it works. A good reference can be found here:

http://onertipaday.blogspot.com/2007/05/rotating-distribution-plot-by-90.html



Hope that helped more than my last post!

Jannis

--- Ralf B <ralf.bierig at gmail.com> schrieb am Di, 21.9.2010:

> Von: Ralf B <ralf.bierig at gmail.com>
> Betreff: [R] Combined plot: Scatter + density plot
> An: "r-help Mailing List" <r-help at r-project.org>
> Datum: Dienstag, 21. September, 2010 05:34 Uhr
> Hi,
> 
> in order to save space for a publication, it would be nice
> to have a
> combined scatter and density plot similar to what is shows
> on
> 
> http://addictedtor.free.fr/graphiques/RGraphGallery.php?graph=78
> 
> I wonder if anybody perhaps has already developed code for
> this and is
> willing to share. This is the reproducible code for the
> histogram
> version obtained from the site:
> 
> def.par <- par(no.readonly = TRUE) # save default, for
> resetting...
> x <- pmin(3, pmax(-3, rnorm(50)))
> y <- pmin(3, pmax(-3, rnorm(50)))
> xhist <- hist(x, breaks=seq(-3,3,0.5), plot=FALSE)
> yhist <- hist(y, breaks=seq(-3,3,0.5), plot=FALSE)
> top <- max(c(xhist$counts, yhist$counts))
> xrange <- c(-3,3)
> yrange <- c(-3,3)
> nf <- layout(matrix(c(2,0,1,3),2,2,byrow=TRUE), c(3,1),
> c(1,3), TRUE)
> par(mar=c(3,3,1,1))
> plot(x, y, xlim=xrange, ylim=yrange, xlab="", ylab="")
> par(mar=c(0,3,1,1))
> barplot(xhist$counts, axes=FALSE, ylim=c(0, top), space=0)
> par(mar=c(3,0,1,1))
> barplot(yhist$counts, axes=FALSE, xlim=c(0, top), space=0,
> horiz=TRUE)
> 
> par(def.par)
> 
> I am basically stuck from line 6 where the bin information
> from the
> histogram is used for determining plotting sizes. Density
> are
> different and don't have (equal) bins and their size would
> need to be
> determined differently. I wonder if somebody here has
> created such a
> diagram already and is willing to share ideas/code/pointers
> to similar
> examples. Your effort is highly appreciated.
> 
> Thanks a lot,
> Ralf
> 
> ______________________________________________
> 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.
>