R ggplot2 stat_density2d网格

Question

在R，ggplot2中，如何将stat_density2d几何对象中使用的网格更改为受三角形区域约束？

我正在绘制不是标准ggplot2类型的三元图，轮廓在轴上溢出。实际三元图的代码非常冗长，但实质上，我需要将等高线约束为等边三角形。

以下是重现问题的示例代码：

library(ggplot2)

triangle <- data.frame(x=c(0,.5,1),y=c(0,1,0))
dummy <- data.frame(x=c(0.25,0.50,0.75),y=c(0.25,0.75,0.25))

p <- ggplot() + geom_polygon(data=triangle,aes(x,y),fill="transparent",color="black",size=1)
p <- p + geom_point(data=dummy,aes(x,y))
p <- p + stat_density2d(data=dummy,aes(x,y)) + theme_bw() +
  theme(panel.grid.major = element_blank(), panel.grid.minor = element_blank(), 
        panel.border = element_blank(), axis.ticks = element_blank(), 
        axis.text.x = element_blank(), axis.text.y = element_blank(),
        axis.title.x = element_blank(), axis.title.y = element_blank())
p

我认为这是由于网格与绘图区域的原始笛卡尔坐标边界相关联，为了演示，请考虑以下代码：

library(ggplot2)

triangle <- data.frame(x=c(0,.5,1),y=c(0,1,0))
dummy <- data.frame(x=c(0.25,0.50,0.75),y=c(0.25,0.75,0.25))

p <- ggplot() + geom_polygon(data=triangle,aes(x,y),fill="transparent",color="black",size=1)
p <- p + geom_point(data=dummy,aes(x,y))
p <- p + stat_density2d(data=dummy,aes(x,y),geom='tile',contour=F,fill="transparent",color="magenta") + theme_bw() +
  theme(panel.grid.major = element_blank(), panel.grid.minor = element_blank(), 
        panel.border = element_blank(), axis.ticks = element_blank(), 
        axis.text.x = element_blank(), axis.text.y = element_blank(),
        axis.title.x = element_blank(), axis.title.y = element_blank())
p

产生以下内容：

那么如何将网格更改为受三角形约束。而不是原始的笛卡尔绘图界限？

Answer 1

自发布ggtern包以来，现在可以使用以下简单代码实现此结果：

#load the library
library(ggtern)

#Create dummy data
df <- data.frame(T=c(.8,.1,.1),
                 L=c(.1,.8,.1),
                 R=c(.1,.1,.8))

#Now plot the diagram
ggtern(data=df,aes(y=T,x=L,z=R)) + geom_point() + geom_density2d() + theme_bw()

Answer 2

这将是一个有点粗糙的选择，只在这里部分证明：

d1 <- data.frame(x = c(0,0.5,0),y = c(0,1,1))
p + geom_polygon(data = d1,aes(x = x,y = y),fill = "white",colour = "black")

您可能会调查其他一些想法here和here，但我认为它们可能涉及丢弃ggplot2。

从理论上讲，你可以使用grid.path做一些奇特的事情，但那里的困难（一如既往）正在干净地在设备和情节坐标系之间移动。