2

私はパッケージscatter3Dから苦労しています。plot3D

コルバーに透明度のしきい値を設定できるようにしたいと思います。たとえば、以下では、すべてのポイントcolvar<100(赤から緑の中央まで) を完全に透明にする必要があります。コルキーの半分（および対応するポイント;））が消えない理由がわかりません。アイデアはありますか？

ここに画像の説明を入力

ここに私のデータがあります: http://s000.tinyupload.com/?file_id=00763309738825363461

そして私のコード：

err_tab2=read.table("Err_tab.txt",header=T)

colfunc<-colorRampPalette(c("red","yellow","green","turquoise2","blue"))


  par(mfrow = c(1, 1))
  panelfirst <- function(pmat) {
    zmin <- min(err_tab2$z)
    XY <- trans3D(err_tab2$x, err_tab2$y,
                  z = rep(zmin, nrow(err_tab2)), pmat = pmat)
    scatter2D(XY$x, XY$y, colvar = err_tab2$colour, pch = ".",
              cex = 1, add = TRUE, colkey = FALSE)
    xmin <- min(err_tab2$x)
    XY <- trans3D(x = rep(xmin, nrow(err_tab2)), y = err_tab2$y,
                  z = err_tab2$z, pmat = pmat)
    scatter2D(XY$x, XY$y, colvar = err_tab2$colour, pch = ".",
              cex = 1, add = TRUE, colkey = FALSE)
  }

  with(err_tab2, scatter3D(x = x, y = y, z = z, colvar= (colour), col=alpha.col(colfunc(length(colour)),
                                                                                alpha =colour<100),
                           pch = 16, cex = 1, xlab = "x", ylab = "y",
                           zlab = "z", clab = c(""),
                           main = paste("Model"), ticktype = "detailed",
                           panel.first = panelfirst, theta = 330,phi=10, y = 2,
                           colkey = list(length = 0.5, width = 0.5, cex.clab = 0.75))
  )

注：実際には、次のような各ポイントの条件が必要です。

alpha=(err_tab2$colour-min(err_tab2$colour))/(max(err_tab2$colour)-min(err_tab2$colour))>0.2

私の考えでは、これにより、「色」がセットの最初の 20% に属するすべてのポイントが非表示になるはずです。F=0=100%transparencyとT=1=100%opacity. _

しかし、これを試してみると、正反対の結果が得られます（前の問題に追加されたと思います）：

ここに画像の説明を入力

何か案が ？:)

4

1 に答える 1

0

これが私が探していた答えです。セットの最初の 20% の値を持つポイントのみをフィルタリングするには:

col_level=sort(as.numeric(levels(as.factor(err_tab2$colour))))

with(err_tab2, scatter3D(x = x, y = y, z = z, colvar= (colour), col=alpha.col(colfunc(length(colour)),
                               alpha =(col_level-min(col_level))/(max(col_level)-min(col_level))<0.20)),
                               pch = 16, cex = 1, xlab = "x", ylab = "y",
                               zlab = "z", clab = c(""),
                               main = paste("Model"), ticktype = "detailed",
                               panel.first = panelfirst, theta = 330,phi=10, y = 2,
                               colkey = list(length = 0.5, width = 0.5, cex.clab = 0.75))
)
于 2016-04-06T08:07:15.513 に答える