我正在尝试创建类似于其中之一的 3D(4D?)狄利克雷概率密度函数图(来自Wikipedia):
我的数据由 3 列和 100000 行组成,其中每行总和为从rdirichlet
函数获得的 1。我可以使用创建热图版本DirichletReg
,它大致相当于上面维基百科图中右下角的图,如下所示:
library(gtools)
library(DirichletReg)
dat <- rdirichlet(100000, alpha = c(10,1,10))
plot(DR_data(dat), a2d = list(colored = TRUE, c.grid = FALSE, col.scheme = c("entropy")))
或者我可以创建一个类似于左上角 Wiki 图的统一热图版本,如下所示:
dat2 <- rdirichlet(100000, alpha = c(1,1,1))
plot(DR_data(dat2), a2d = list(colored = TRUE, c.grid = FALSE, col.scheme = c("entropy")))
我的问题:是否可以像维基图一样制作三元图 3D,而不是热图?我曾尝试使用plotly
和plot3d
包来做到这一点,但无法接近。任何帮助将非常感激。