次のように作成された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")
ただし、これは失敗します。古いプロットが上書きされるだけです。新しいプロットを重ねる方法はありますか?