如何在hist3D上打印標籤上的希臘字母？

Question

我正在使用hist3D繪製3D直方圖，並且想要更改xlab，以便它包含帶有下標的希臘字母。

我用

xlab=expression(theta[1]). 

它不工作，我只是得到的標籤字符串 「[1] THETA」。另一方面，它可以正常工作，只需要經常繪圖命令。

那麼如何才能將下標希臘字母引入hist3D圖？

（hist3D來自PLOT3D格式庫）

Answer 1

一種解決方法是使用scatterplot3d：

z <- seq(-10, 10, 0.01) 
x <- cos(z)+1 
y <- sin(z)+1 
scatterplot3d(volcano, highlight.3d=TRUE, col.axis="blue", 
          col.grid="gray", main="scatterplot-greeks", pch=21, 
          xlab=expression(theta[1]), 
          ylab=expression(beta[2]), 
          zlab=expression(alpha[3]),cex.lab=1.5) 

Answer 2

的hist3D功能不接受軸標籤表達式。一種解決方法是對theta，θ使用Unicode字符。

# example data 
library(plot3D) 
x <- seq(-pi, pi, by = 0.2) 
y <- seq(-pi, pi, by = 0.3) 
grid <- mesh(x, y) 
z <- with(grid, cos(x) * sin(y)) 

hist3D(z = z, x = x, y = y, border = "black", xlab = "θ1") 

這將顯示一個符號THETA，但沒有下標1

Answer 3

另一個解決方法是使用text3D函數從同一plot3D包。它接受數學表達式作爲標籤，所以我們可以照常編寫expression(theta[1])。在此解決方法唯一的小問題是要找到

# example data 
library(plot3D) 
x <- seq(-pi, pi, by = 0.2) 
y <- seq(-pi, pi, by = 0.3) 
grid <- mesh(x, y) 
z <- with(grid, cos(x) * sin(y)) 

hist3D(z = z, x = x, y = y, border = "black", ticktype = "detailed", xlab = "", ylab="y", zlab="z") 
text3D(2, -5, 0, labels = expression(theta[1]), add = TRUE, adj = 1) 

我發現2,-5,0是這個例子相當好我們的標籤的合適座標。標籤不像標準標籤那樣旋轉。