r - persp3d图的格式

我有以下3D图:

使用data，我使用以下代码创建了它:

library(rugarch)
library(rgl)
library(fGarch)

fd <- as.data.frame(modelfit, which = 'density')
color <- rgb(85, 141, 85, maxColorValue=255)

x <- seq(-0.2, 0.2, length=100)
y <-c(1:2318)

f <- function(s, t) {
 dged(s,mean=fd[t,'Mu'],sd=fd[t,'Sigma'],nu=fd[t,'Shape'])

}

z <- outer(x, y, f)
persp3d(x, y, z, theta=50, phi=25, expand=0.75, col=color,
      ticktype="detailed", xlab="", ylab="time", zlab="",axes=TRUE)

如何根据z值获得着色？我看了不同的解决方案，例如这个one，但是在这种情况下，我无法根据z值创建颜色。
根据this thread的解决方案如下:

nrz <- nrow(z)
ncz <- ncol(z)
jet.colors <- colorRampPalette( c("#ffcccc", "#cc0000") )
# Generate the desired number of colors from this palette
nbcol <- 100
color <- jet.colors(nbcol)

# Compute the z-value at the facet centres
zfacet <- z[-1, -1] + z[-1, -ncz] + z[-nrz, -1] + z[-nrz, -ncz]
# Recode facet z-values into color indices
facetcol <- cut(zfacet, nbcol)

  persp3d(x, y, z, theta=50, phi=25, expand=0.75, col=color[facetcol],
          ticktype="detailed", xlab="", ylab="time", zlab="",axes=TRUE)

但这并不能带来良好的结果，因为它无法为绘图添加适当的颜色。我想让表面的尖刺成为例如以红色表示，低值例如呈蓝色且平滑过渡良好，但是这种颜色可以使切片着色，因此取决于时间吗？因此，极大的尖峰应在其尖峰处涂成红色，并在底部标出值，例如绿色。我怎么能得到这个？

编辑:我找到了关于轴上日期的上一个问题的解决方案，剩下的唯一问题是取决于z值的适当着色。

最佳答案

试试这个:

nbcol = 100
color = rev(rainbow(nbcol, start = 0/6, end = 4/6))
zcol  = cut(z, nbcol)
persp3d(x, y, z, theta=50, phi=25, expand=0.75, col=color[zcol],
        ticktype="detailed", xlab="", ylab="time", zlab="",axes=TRUE)

如果您希望按时间显示颜色(因此尖峰始终为红色)，则可以为每个时间片设置颜色:

mycut = function(x, breaks) as.numeric(cut(x=x, breaks=breaks)) # to combine different factors
zcol2 = as.numeric(apply(z,2, mycut, breaks=nbcol))
persp3d(x, y, z, theta=50, phi=25, expand=0.75, col=color[zcol2],
        ticktype="detailed", xlab="", ylab="time", zlab="",axes=TRUE)

您已经知道如何正确编辑轴。

关于r - persp3d图的格式，我们在Stack Overflow上找到一个类似的问题：https://stackoverflow.com/questions/17258787/