1

行ごとに異なる列に配置された値のマトリックスがあります。私のコードでやりたいこと:

  1. 行を反復処理 -> 列の値 < しきい値 (例: 1) かどうかを確認します
  2. 行内で、その条件を満たす列が 80% ある場合、行を保持します。それ以外の場合は、行全体を削除します。

これまでのコード:

myfilt <- function(t,x){
             if ((length(which(t[x,] > 1)) / 60) >= 0.8){
               return(1)
             }else{
               return(0)
             }
          }

y=c()
for(i in 1:length(t[,1])){
   y = c(y,myfilt(t,i))
}

しかし、t[v,]すべての行を同じ値で印刷すると:(何が間違っているのかわかりません。また、これを行うためのより短い方法がある場合は、お知らせください。

PS:ここで「t」は私がテストしているマトリックスの名前です

4

1 に答える 1

2

これを行う方法は次のとおりです。

## Parameters
threshold <- 0.8
perc.to.keep <- 0.5
## Example Matrix
set.seed(1337)
m <- matrix(rnorm(25,1,1),nrow=5,ncol=5)

#           [,1]      [,2]        [,3]       [,4]      [,5]
# [1,] 1.7122837 0.8383025 -0.02718379  2.2157099 2.1291008
# [2,] 0.2462742 2.4602621 -0.04117532 -0.6214087 1.4501467
# [3,] 1.0381899 3.0094584  0.12937698  0.9314247 1.0505864
# [4,] 2.1784211 0.9220618  1.85313022  0.9370171 0.8756698
# [5,] 0.8467962 2.3543421  0.37723981  2.0757077 1.9120115

test <- m < threshold
sel <- apply(test,1,function(v) sum(v)/length(v)) < perc
m[sel,]

#           [,1]      [,2]        [,3]      [,4]      [,5]
# [1,] 1.7122837 0.8383025 -0.02718379 2.2157099 2.1291008
# [2,] 1.0381899 3.0094584  0.12937698 0.9314247 1.0505864
# [3,] 2.1784211 0.9220618  1.85313022 0.9370171 0.8756698
# [4,] 0.8467962 2.3543421  0.37723981 2.0757077 1.9120115
于 2013-03-22T16:02:42.553 に答える