R 中的二维核密度图和透明等高线填充

Question

我尝试使用R中的ks包来绘制带有透明填充的二维核密度图这是我的代码。

library(ks)
B<-matrix(c(3.02,2.71,8.05,1.36,1.18,2.64,1.48,5.51,3.25,2.66,0.15,2.25,3.09,2.49,3.37,4.56,0.15,0.57,2.6,0.24,0.27,0.1,0.34,0.24,0.24,0.24,0.23,0.34,0.36,0.1,12,0.88,0.15,0.2,0.1,11.19,0.1,0.1,0.26,22.3),nrow=20,ncol=2)
f<-kde(B)
png(filename = "KDE_test.png",  units="in", width=5, height=5, res=100)
plot(f,display="filled.contour2",xlim=c(-1,5),ylim=c(-3.5,5),lwd=0,lwd.fc=0.1,drawlabels=FALSE,col=c("transparent",rgb(0,0,255, maxColorValue=255, alpha=50),rgb(0,0,255, maxColorValue=255, alpha=100),rgb(0,0,255, maxColorValue=255, alpha=150)))
dev.off()

但是，我最终得到了一个边缘参差不齐的 KDE 图。抱歉，我无法上传我的图片，但请运行此代码，您将看到我正在谈论的图片。

有谁知道如何解决这个问题还是我做错了什么？

非常感谢。

Answer 1

您可以使用 kde 的 gridsize 选项。

library(ks)
B <- matrix(c(3.02,2.71,8.05,1.36,1.18,2.64,1.48,5.51,3.25,2.66,0.15,2.25,3.09,
  2.49,3.37,4.56,0.15,0.57,2.6,0.24,0.27,0.1,0.34,0.24,0.24,0.24,0.23,
  0.34,0.36,0.1,12,0.88,0.15,0.2,0.1,11.19,0.1,0.1,0.26,22.3),
  nrow=20, ncol=2)

f<-kde(B, gridsize=1024)

plot(f,display="filled.contour2", xlim=c(-1,5), ylim=c(-3.5,5), 
    lwd=0, lwd.fc=0.1, drawlabels=FALSE, 
    col=c("transparent", 
       rgb(0,0,255, maxColorValue=255, alpha=50), 
       rgb(0,0,255, maxColorValue=255, alpha=100), 
       rgb(0,0,255, maxColorValue=255, alpha=150)))

R 中的二维核密度图和透明等高线填充

2D Kernel Density Plot and Transparent Contour Fills in R

r

transparency

contour

fill

kernel-density