У меня есть следующий трехмерный график:
С помощью данных я создал его со следующим кодом:
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? Я рассмотрел различные решения, например. этот один, но я не мог создать раскраску в зависимости от значений z в этом случае. Решение в соответствии с этой нитью было бы следующим:
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.