如何在R（瀑布图）中制作3D线图

Question

我想从我的data在R（XYYY）中创建一个瀑布图。

到目前为止，我使用此代码：

load("myData.RData")
ls()
dim(data)

##matrix to xyz coords
library(reshape2)
newData <- melt(data, id="Group.1")
dim(newData)
head(newData)
tail(newData)
newDataO <- newData[c(2,1,3)]
head(newDataO)

##color scale for z axis
myColorRamp <- function(colors, values) {
v <- (values - min(values))/diff(range(values))
x <- colorRamp(colors)(v)
rgb(x[,1], x[,2], x[,3], maxColorValue = 255)
}

cols <- myColorRamp(c("darkblue","yellow","darkorange","red","darkred"),newDataO$value)

##3D scatter
library(rgl)
plot3d(newDataO$variable, newDataO$Group.1, newDataO$value, xlab="", ylab="", zlab="",      type="p", col=cols, box=FALSE, axes=FALSE)

rgl.postscript("persptrial_060514.eps","eps")

得到这个情节：

https://dl.dropboxusercontent.com/u/14906265/persptrial_060514.jpg

我也在2d中使用此选项和多边形，但结果不能正确显示两个图之间的差异效果（左对比）。

我不知道像persp3d这样的东西是否可以完成这项工作，但我对编写代码来实现它并不熟悉。任何帮助将非常感谢。

Answer 1

在我看来，在R中进行瀑布绘图的最简单方法是在循环中手动添加所有行。

library(rgl)

# Function to plot
f <- function(x, y) sin(10 * x * y) * cos(4 * y^3) + x

nx <- 30
ny <- 100
x <- seq(0, 1, length = nx)
y <- seq(0, 1, length = ny)
z <- outer(x, y, FUN = f)

# Plot function and add lines manually
surface3d(x, y, z, alpha = 0.4)
axes3d()
for (i in 1:nx) lines3d(x[i], y, z[i, ], col = 'white', lwd = 2)