R help - Oct 2017 - how to overlay 2d pdf atop scatter plot using ggplot2

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]]

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.

How about geom_contour()?

Am So., 8. Okt. 2017, 20:52 schrieb Ranjan Maitra <maitra at email.com>:
> 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.
>
> ______________________________________________
> 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.
>
	[[alternative HTML version deleted]]