大家晚上好。我正在尝试编写一个直方图函数,以3D形式显示给定的直方图。另外,我只能将persp()函数与网格包一起应用。问题是我们定义了theta和phi,我也想添加额外的gamma参数,以使对象在默认平面上旋转。我想获得一个函数的表格:
hist3d(X, alpha, beta, gamma, v1, v2)
我的函数当前如下所示:
hist3d <- function(X, alpha, beta, gamma, v1, v2)
{
l1 <- length(X$v1)
l2 <- length(X$v2)
x1r <- range(X$v1)
x1seq <- seq(x1r[1], x1r[2], length=l1)
x2r <- range(X$v2)
x2seq <- seq(x2r[1], x2r[2], length=l2)
z <- outer(x1seq, x2seq)
persp(x=x1seq, y=x2seq, z=z, theta=alpha, phi=beta)
}
我应该对theta和phi使用六面体对角线公式吗?
gamma <- ((alpha*sqrt(2))/(beta*sqrt(2))
有人可以为我提供使用技巧的任何提示吗?