我使用以下代码创建具有自定义带宽的内核密度估计的三元图。
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")
现在我想创建一个类似的图,但内核密度具有不同的机密级别而不是带宽。它应该类似于这篇科学论文中的这张图。https://www.sciencedirect.com/science/article/pii/S1048984312000811
谁能建议我如何实现这一目标?