Question

我想在3D中看到一个类似于：

的图形

此图形由Matlab中的代码生成：

clear all  
  % generate synthetic data
  data=[randn(500,2);
      randn(500,1)+3.5, randn(500,1);];
  % call the routine, which has been saved in the current directory 
  [bandwidth,density,X,Y]=kde2d(data);
  % plot the data and the density estimate
  contour3(X,Y,density,50), hold on
  plot(data(:,1),data(:,2),'r.','MarkerSize',5)

现在，我在R中的实现这是我的数据：

xc1 <- matrix(rnorm(200)*0.5+2, 100, 2) 
xc2 <- matrix(rnorm(200)*0.5+4, 100, 2) 
xc <- rbind(xc1,xc2)

将来，我将使用此逻辑来进行数据的通用分发，这就是我想要使用kde的原因。所以：

fhat<-kde(x=xc)

使用情节：

plot(fhat,display="persp")

My plot

但是，我不明白这个＆＃34;情节＆＃34;做了。查看fhat，我看不到xgrid，ygrid，密度轴。

我希望使用persp3d使用特定网格绘制相同的图形：

seqi <- seq(0,6,0.1)
seqj <- seq(0,6,0.1)

persp3d(seqi, seqj, z1,col = 'red')

如何用kde获得z1？

Tnx非常关注

使用核密度，R-cran的3D图

0 个答案: