在R中使用线框在3d中进行绘制

Question

I am trying to plot in 3d with wireframe, here is my code:

mecdf = function (nr, nc, x, u) # mecdf is the bivariate cumulative empiric function 
{  k = rep (TRUE, nr)
   for (j in 1:nc) k = k & (x [,j] <= u [j])
   sum (k) / nr
}

xc = round(runif(100), 2)
yc = round(runif(100), 2)
Da = cbind(xc, yc)

bcdfa<-rep(NA,100)
for (i in 1:100) {bcdfa[i]=mecdf(nrow(Da),ncol(Da),Da,Da[i,])}    
bcdfa

x<-cbind(xc,yc,bcdfa) # bcdfa is the value of the bcdfa in every (xc,yc) point



 ...........

Is it possible to use the wireframe or another function with these data to represent the bivariate cumulative distribution function ?

Answer 1

Here's one way.

library(akima)
library(plot3D)
 bcfoo<-interp(xc,yc,bcdfa)

 persp3D(bcfoo$x,bcfoo$y,bcfoo$z)

 #plot points to verify overall shape

 scatter3D(xc,yc,bcdfa)