2013-03-06 4 views
5

データフレームがあり、値が0より大きい場合は最初の56個の列のすべてのデータポイントを2進化したいと思います1に設定します。それ以外の場合は0に設定されます。これを行う簡単な方法はありますか?ベクトル化ifelseはあなたが行うことができます使用Rのデータフレーム内のデータを2進化する必要があります

+1

[この質問](http://stackoverflow.com/q/14526429/1270695)とさまざまな回答と、速度が懸念される場合のベンチマークも参照できます。この質問との違いは、列のサブセットのみで作業していることです。 – A5C1D2H2I1M1N2O1R2T1

+0

TRUE/FALSEではなく0/1の値が必要な特別な理由はありますか?あなたは、結果をかなりの数の論理的なものとして残してしまうことができます。 – Dason

答えて

8

m[,1:56] <- ifelse(m[,1:56] > 0,1,0) 

たとえば、私たちは、小さな行列でこれをテストすることができます。

m <- matrix(sample(c(-2,2),5*3,rep=T),ncol=5,nrow=3,byrow=T) 
> m 
    [,1] [,2] [,3] [,4] [,5] 
[1,] 2 2 2 2 -2 
[2,] 2 2 -2 2 -2 
[3,] 2 2 2 2 2 
> m[,2:5] <- ifelse(m[,2:5] > 0,1,0) 
> m 
    [,1] [,2] [,3] [,4] [,5] 
[1,] 2 1 1 1 0 
[2,] 2 1 0 1 0 
[3,] 2 1 1 1 1 
+2

私の考えはまさに! (私は別の代替案を投稿しなければならなかったので) – A5C1D2H2I1M1N2O1R2T1

+0

@AnandaMahto、amazing !! –

+0

@RicardoSaporta、バイナリの結果が必要なだけでなく、実際には 'ifelse'が好きです。 – A5C1D2H2I1M1N2O1R2T1

5

あなたがTRUEFALSEがに等しいという事実を利用することができます"1"と "0"を入力してください。

set.seed(1) 
mydf <- data.frame(matrix(rnorm(100), nrow = 10)) 
mydf[, 1:5] <- (mydf[, 1:5] > 0) + 0 
mydf 
# X1 X2 X3 X4 X5   X6   X7   X8   X9  X10 
# 1 0 1 1 1 0 0.3981059 2.40161776 0.475509529 -0.5686687 -0.5425200 
# 2 1 1 1 0 0 -0.6120264 -0.03924000 -0.709946431 -0.1351786 1.2078678 
# 3 0 0 1 1 1 0.3411197 0.68973936 0.610726353 1.1780870 1.1604026 
# 4 1 0 0 0 1 -1.1293631 0.02800216 -0.934097632 -1.5235668 0.7002136 
# 5 1 1 1 0 0 1.4330237 -0.74327321 -1.253633400 0.5939462 1.5868335 
# 6 0 0 0 0 0 1.9803999 0.18879230 0.291446236 0.3329504 0.5584864 
# 7 1 0 0 0 1 -0.3672215 -1.80495863 -0.443291873 1.0630998 -1.2765922 
# 8 1 1 0 0 1 -1.0441346 1.46555486 0.001105352 -0.3041839 -0.5732654 
# 9 1 1 0 1 0 0.5697196 0.15325334 0.074341324 0.3700188 -1.2246126 
# 10 0 1 1 1 1 -0.1350546 2.17261167 -0.589520946 0.2670988 -0.4734006 

+0の考え方は、単に論理値TRUEFALSEを数値に等しくすることです。 すべての行列で作業していて、as.numeric(mydf > 0)を使用している場合、結果のベクトルを行列に再変換する必要があります。しかし、この場合、これは完全にうまくいく(@Dasonによって指摘されているように)。

mydf[, 1:5] <- as.numeric(mydf[, 1:5] > 0) 
+1

ここで 'as.numeric'を使ってみましたか?それは私のための交換でうまく動作します...私たちが交換をしていない場合、違いは重要ですが、このようなもののためにはうまく動作します。 – Dason

+0

@Dason、Good catch! (いいえ、私は試していない - ちょうど私は同じ長さではない行についてのエラーを得たと仮定した)。良い古いリサイクル、私は推測するだろう。私は更新します。ありがとう。 – A5C1D2H2I1M1N2O1R2T1

1

pminおよびpmaxを使用するアプローチ。 (実際にはお勧めしません)

pmin(pmax(m[,2:5], 0),1) 

しかし、それはAnandaMahtoのソリューション@いくつかのベンチマークに

ag <- function() ifelse(m[,2:5] > 0,1,0) 
mn <- function()pmin(pmax(m[,2:5], 0),1) 
am <- function() (m[, 2:5] > 0) + 0 
am2 <- function() as.numeric((m[, 2:5] > 0)) 

library(microbenchmark) 
microbenchmark(ag(),mn(), am(), am2()) 
## Unit: microseconds 
## expr min  lq median  uq  max neval 
## ag() 19.888 20.712 21.9375 22.6430 39.548 100 
## mn() 50.135 51.172 52.2530 53.1055 113.854 100 
## am() 3.076 3.406 4.1755 4.6030 7.912 100 
## am2() 2.623 2.989 3.4640 4.0135 6.995 100 

を追加することができ、さらに高速as.numericアプローチで、明確な勝者です!

+0

最近更新された答えに 'as.numeric'オプションを追加してください(私の宿題をやってくれたDasonに感謝します)。以前の経験によれば、それはさらに速くなるはずです。 – A5C1D2H2I1M1N2O1R2T1

+0

@AnandaMahto - done。はい。 – mnel

関連する問題