R -

Question

それぞれ2面の二変量3Dの法線とその輪郭をプロット私はMASSパッケージで遊んでてきたし、単に例えば画像やPAR（新しい= TRUE）を使用して、通常の2変量をプロットすることができます

# lets first simulate a bivariate normal sample 
library(MASS) 
bivn <- mvrnorm(1000, mu = c(0, 0), Sigma = matrix(c(1, .5, .5, 1), 2)) 
bivn2 <- mvrnorm(1000, mu = c(0, 0), Sigma = matrix(c(1.5, 1.5, 1.5, 1.5), 2)) 

# now we do a kernel density estimate 
bivn.kde <- kde2d(bivn[,1], bivn[,2], n = 50) 
bivn.kde2 <- kde2d(bivn2[,1], bivn[,2], n = 50) 

# fancy perspective 
persp(bivn.kde, phi = 45, theta = 30, shade = .1, border = NA) 
par(new=TRUE) 
persp(bivn.kde2, phi = 45, theta = 30, shade = .1, border = NA) 

これはあまり良く見えませんが、私はちょうど軸と材料で遊んでいなければならないと思います。 しかし、輪郭で同様のアプローチを試してみると、プロットは重ならない。それらは単に置き換えられます：

# fancy contour with image 
image(bivn.kde); contour(bivn.kde, add = T) 
par(new=TRUE) 
image(bivn.kde2); contour(bivn.kde, add = T) 

これは私が欲しいものに最適なアプローチですか、それとも自分自身では難しいですか？どんな提案も大歓迎です。ありがとうございました！

Answer 1

おそらくrglライブラリを使用できます。インタラクティブな3Dプロットを作成することができます。

require(rgl) 

col1 <- rainbow(length(bivn.kde$z))[rank(bivn.kde$z)] 
col2 <- heat.colors(length(bivn.kde2$z))[rank(bivn.kde2$z)] 
persp3d(x=bivn.kde, col = col1) 
with(bivn.kde2, surface3d(x,y,z, color = col2)) 

あなたは二つの面の間の差をプロットしたいなら、あなたは以下のような何かを行うことができます。

res <- list(x = bivn.kde$x, y = bivn.kde$y, z = bivn.kde$z - bivn.kde2$z) 
col3 <- heat.colors(length(res$z))[rank(res$z)] 
persp3d(res, col = col3)