0

我使用以下代码创建具有自定义带宽的内核密度估计的三元图。

require(ggplot2)
require(ggtern) 

set.seed(1) 

mydata <- data.frame(
  x = runif(100, min = 0.3, max = 0.6),
  y = runif(100, min = 0.2, max = 0.65),
  z = runif(100, min = 0.5, max = 0.8))   


ggtern(data= mydata,aes(x=x,y=y,z=z)) + stat_density_tern(geom="polygon",n=800,bins=5,
                                                                              aes(fill  = ..level..,
                                                                                  alpha = ..level..))+
  theme_bw() +
  theme_hidetitles() +
  scale_fill_gradient(low = "blue",high = "blue3")  +
  guides(color = "none", fill = "none", alpha = "none")

bin = 5 的 KDE 绘图

现在我想创建一个类似的图,但内核密度具有不同的机密级别而不是带宽。它应该类似于这篇科学论文中的这张图。https://www.sciencedirect.com/science/article/pii/S1048984312000811

谁能建议我如何实现这一目标?

在此处输入图像描述

4

0 回答 0