係数を数値または整数に変換すると、数値としての値ではなく、基礎となるレベル コードが取得されます。
f <- factor(sample(runif(5), 20, replace = TRUE))
## [1] 0.0248644019011408 0.0248644019011408 0.179684827337041
## [4] 0.0284090070053935 0.363644931698218 0.363644931698218
## [7] 0.179684827337041 0.249704354675487 0.249704354675487
## [10] 0.0248644019011408 0.249704354675487 0.0284090070053935
## [13] 0.179684827337041 0.0248644019011408 0.179684827337041
## [16] 0.363644931698218 0.249704354675487 0.363644931698218
## [19] 0.179684827337041 0.0284090070053935
## 5 Levels: 0.0248644019011408 0.0284090070053935 ... 0.363644931698218
as.numeric(f)
## [1] 1 1 3 2 5 5 3 4 4 1 4 2 3 1 3 5 4 5 3 2
as.integer(f)
## [1] 1 1 3 2 5 5 3 4 4 1 4 2 3 1 3 5 4 5 3 2
paste実際の値を取得するには、次の方法に頼る必要があります。
as.numeric(paste(f))
## [1] 0.02486440 0.02486440 0.17968483 0.02840901 0.36364493 0.36364493
## [7] 0.17968483 0.24970435 0.24970435 0.02486440 0.24970435 0.02840901
## [13] 0.17968483 0.02486440 0.17968483 0.36364493 0.24970435 0.36364493
## [19] 0.17968483 0.02840901
係数を数値に変換するより良い方法はありますか?
ベストアンサー1
の警告セクションを参照してください?factor:
特に、
as.numeric因子に を適用しても意味がなく、暗黙的な強制によって発生する可能性があります。因子をf元の数値に近似的に変換するには、as.numeric(levels(f))[f]が推奨され、 よりもわずかに効率的ですas.numeric(as.character(f))。
Rに関するFAQ同様のアドバイスをしている。
as.numeric(levels(f))[f]が よりも効率的なのはなぜですかas.numeric(as.character(f))?
as.numeric(as.character(f))は実質的に なので、値に対してではなく、値as.numeric(levels(f)[f])に対して数値への変換を実行します。速度の違いは、レベル数が少ない長いベクトルの場合に最も顕著になります。値がほとんど一意である場合、速度に大きな違いはありません。変換をどのように行うかに関係なく、この操作がコードのボトルネックになる可能性は低いため、あまり心配する必要はありません。length(x)nlevels(x)
いくつかのタイミング
library(microbenchmark)
microbenchmark(
as.numeric(levels(f))[f],
as.numeric(levels(f)[f]),
as.numeric(as.character(f)),
paste0(x),
paste(x),
times = 1e5
)
## Unit: microseconds
## expr min lq mean median uq max neval
## as.numeric(levels(f))[f] 3.982 5.120 6.088624 5.405 5.974 1981.418 1e+05
## as.numeric(levels(f)[f]) 5.973 7.111 8.352032 7.396 8.250 4256.380 1e+05
## as.numeric(as.character(f)) 6.827 8.249 9.628264 8.534 9.671 1983.694 1e+05
## paste0(x) 7.964 9.387 11.026351 9.956 10.810 2911.257 1e+05
## paste(x) 7.965 9.387 11.127308 9.956 11.093 2419.458 1e+05