在R

Question

给出以下多元正态样本：

mu=rep(0,2)
Sigma=matrix(c(1,0,0,1),2,2)
require(MASS)
X=mvrnorm(n=100,mu,Sigma)

我想计算的是参数密度函数。过去这样做是：

require(rgl)
require(mnormt)
zX=dmnorm(as.matrix(X),mean=colMeans(X),varcov=cov(X))
plot3d(X[,1],X[,2],zX) 

但是这会弹出点图（在RStudio之外）：

我想知道如何获得图形化的RStudio绘制窗口，该窗口也可以以此点绘制表面。我可以使用非参数密度来做到这一点：

ngrid = 50
Zbiksm = kde2d(X[,1],X[,2],n=ngrid)
Zbiksm_den = Zbiksm$z # valor de la densidad en cada punto del grid
persp(x=Zbiksm$x,y=Zbiksm$y,z=Zbiksm_den,xlab="x",ylab="y",
      zlab=expression(K(x,y)),theta=-35,axes=TRUE,box=TRUE)

但在这种情况下不知道该怎么做。

或者，我也可以将surface3d和散点图图像放在一起。