[R] how to overlay 2d pdf atop scatter plot using ggplot2

Ranjan Maitra maitra at email.com
Sun Oct 8 20:51:43 CEST 2017


Hi,

I am no expert on ggplot2 and I do not know the answer to your question. I looked around a bit but could not find an answer right away. But one possibility could be, if a direct approach is not possible, to draw ellipses corresponding to the confidence regions of the multivariate t density and use geom_polygon to draw this successively? 

I will wait for a couple of days to see if there is a better answer posted and then write some code, unless you get to it first.

Thanks,
Ranjan


On Sun, 8 Oct 2017 09:30:30 -0500 Big Floppy Dog <bigfloppydog at gmail.com> wrote:

> Note: I have posted this on SO also but while the question has been
> upvoted, there has been no answer yet.
> 
> https://stackoverflow.com/questions/46622243/ggplot-plot-2d-probability-density-function-on-top-of-points-on-ggplot
> 
> Apologies for those who have seen it there also but I thought that this
> list of experts may have someone who knows the answer.
> 
> I have the following example code:
> 
> 
> 
> require(mvtnorm)
> require(ggplot2)
> set.seed(1234)
> xx <- data.frame(rmvt(100, df = c(13, 13)))
> ggplot(data = xx,  aes(x = X1, y= X2)) + geom_point() + geom_density2d()
> 
> 
> 
> It yields a scatterplot of X2 against X1 and a KDE contour plot of the
> density (as it should).
> 
> My question is: is it possible to change the contour plot to display
> the contours
> 
> of a two-dimensional density function (say dmvt), using ggplot2?
> 
> The remaining figures in my document are in ggplot2 and therefore I
> am looking for a ggplot2 solution.
> 
> Thanks in advance!
> 
> BFD
> 
> 	[[alternative HTML version deleted]]
> 
> ______________________________________________
> 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.
> 


-- 
Important Notice: This mailbox is ignored: e-mails are set to be deleted on receipt. Please respond to the mailing list if appropriate. For those needing to send personal or professional e-mail, please use appropriate addresses.



More information about the R-help mailing list