我有一个 3d 散点图,如下所示:
library(rgl)
N <- 10000
X <- rnorm(N,0,1)
Y <- rnorm(N,0,1)
Z <- X * Y
want <- Z >0 & X>0
palette <- colorRampPalette(c("blue", "green", "yellow", "red"))
col.table <- palette(256)
col.index <- cut(Z, 256)
plot3d(X,Y,Z, col=col.table[col.index])
grid3d(c("x", "y", "z"))
这工作正常。现在我想覆盖另一个情节,所以我尝试了这个:
par(new=F)
plot3d(X[want],Y[want],Z[want], col="black")
然而这失败了——它只是覆盖了旧的情节。有没有办法覆盖新情节?