从数据帧绘制2D内核密度:设置网格位置、带宽和lims的数量

Question

​

我有两个列的数据帧，分别表示我的"x“和"y”坐标。我想将它们绘制成2D内核密度图。我得了91分。因此，我使用了这个脚本：

x<-read.csv("X_Y.csv", sep=",")
d<-as.matrix(x)
Dgeo<-d[,1]
Dgen<-d[,2] 
dens <- kde2d(Dgeo,Dgen, n=?, lims=c(?,?,?,?))
myPal <- colorRampPalette(c("white","blue","gold", "orange", "red"))
plot(Dgeo, Dgen, pch=20,cex=.5)
image(dens, col=transp(myPal(300),.7), add=TRUE) 
abline(lm(Dgen~Dgeo)) 

我不明白如何估计我已经表示为"?“的值。

你知道吗？

图中显示了我想要获得的内容，右边是我使用n=91和lims=c(-.1，1.5，-.5,4)实际获得的内容。

Answer 1

你可以简单地这样做

library(MASS)
set.seed(0)
dat <- data.frame(x=rnorm(100), y=rnorm(100, 2))
dens <- kde2d(dat$x, dat$y)
image(dens)
contour(dens, add=T)

​

​

不确定颜色的transp函数在哪里