我正在尝试创建一个类似于其中一个(来自Wikipedia)的3D(4D?)Dirichlet概率密度函数图:
我的数据由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")))
我的问题:除了热图之外,还有可能像维基图一样制作三维图吗?
我尝试使用plotly
和plot3d
套餐来做到这一点,但还未能接近。任何帮助将非常感激。