情節3D與R中

我試圖與線框3D繪製線框，這裏是我的代碼：情節3D與R中

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 



...........

是否有可能使用線框或與這些數據的另一個函數來表示二元累計分配功能？

來源

2014-09-12 Quintin Copete

是你的問題「找我的包'lattice'」抑或是「如何轉儲這些變量到函數'格：wireframe'」？ – 2014-09-12 11:48:28

我的問題是如何將這些變量轉儲到函數lattice：wireframe。謝謝 – 2014-09-12 11:56:29

@QuintinCopete請編輯您的問題，使其更清晰。 – 2014-09-12 12:02:59

這是一種方法。

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

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

the surface

#plot points to verify overall shape 

scatter3D(xc,yc,bcdfa)

the data

來源

2014-09-12 13:21:55

非常感謝 – 2014-09-15 07:07:30

回答

相關問題