R中2d單純形的3d圖

Question

有沒有辦法在R中重現以下圖？

編輯

這就是我可以做persp()的基礎R和plot_ly在plotly。 也有點難看。

x <- seq(0,1,0.01) 
y <- seq(0,1,0.01)
f <- function(x,y){ z <- -x - y + 1 }
z <- outer(x,y,f)
z <- ifelse(z<0,NA,z)
persp(x, y, z, theta = 30, phi = 30, expand = 0.5, col = "lightblue")
plot_ly(x=x,y=y,z=z,type="surface") %>% layout(xaxis=list(range=c(0,1)), yaxis=list(range=c(0,1)), zaxis=list(range=c(0,1)))

順便說一句... matplotlib圖在這里獲得： http : //blog.bogatron.net/blog/2014/02/02/visualizing-dirichlet-distributions/

Answer 1

在基本RI中使用persp可以達到以下目標：

persp(0:1, 0:1, 
      matrix(c(1,0,0,NA), nrow=2), 
      col="green", theta=60, 
      xlab= "theta_1", 
      ylab = "theta_2", 
      zlab="theata_3")

但是我不知道該怎么做，包括軸上的希臘符號 。

我persp它變成一個Wiki，以防persp有專業人士想完成這項工作。

Answer 2

這有點難看/仍不完整，但至少顯示了一種獲取希臘標簽的方法。

pp <- persp(0:1, 0:1, 
      matrix(c(2,0,0,NA), nrow=2), 
      col="green", theta=60, 
      xlab= "",
      ylab ="",
      zlab="",
      ticktype="detailed",
      nticks=1)

text(trans3d(0.5,-0.1,-0.1,pp),labels=expression(theta[1]))