r - R-3dで2つの2変量法線とそれらの等高線をそれぞれプロットする

Question

私はMASSパッケージで遊んでいて、たとえばimageとpar（new = TRUE）を使用して、2つの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

Perhaps you can use rgl library. It allows you to create interactive 3d plots.

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))

If you want to plot difference between two surfaces then you can do something like below.

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)