我想删除其任何列data.table中包含的所有行Inf.到目前为止,我一直在使用这种方法:
DT <- data.table(col1 = c(1,2,3), col2 = c(4,Inf,5))
DT[,drop := apply(.SD, 1, function(x) any(is.infinite(x))), by = 1:nrow(DT)]
DT <- DT[(!drop)][,drop:=NULL]
Run Code Online (Sandbox Code Playgroud)
来自这个Stackoverflow问题.但是,这种方法不能很好地扩展到大量数据.是否有更好的方法来删除行Inf?
sha*_*dow 18
您可以rowSums用来检查行的任何元素是否不是有限的.
DT[is.finite(rowSums(DT))]
Run Code Online (Sandbox Code Playgroud)
或者你可以使用一个事实,即Inf * 0是NA与使用complete.cases
DT[complete.cases(DT*0)]
Run Code Online (Sandbox Code Playgroud)
一些基准测试表明,rowSums对于较小的数据集来说,这是最快的,并且complete.cases是大型数据集的最快解决方案.
require(microbenchmark)
microbenchmark(
DT[is.finite(rowSums(DT))]
,
DT[complete.cases(DT*0)]
,
DT[DT[, Reduce('&', lapply(.SD, is.finite))]]
)
##
## nrow(DT) = 3000
## Unit: microseconds
## expr min lq mean median uq max neval cld
## DT[is.finite(rowSums(DT))] 786.797 839.235 864.0215 852.8465 884.756 1021.988 100 a
## DT[complete.cases(DT * 0)] 1265.658 1326.575 1363.3985 1350.0055 1386.377 1898.040 100 c
## DT[DT[, Reduce("&", lapply(.SD, is.finite))]] 1220.137 1275.030 1319.6226 1308.0555 1348.443 1624.023 100 b
##
## nrow(DT) = 300000
## Unit: milliseconds
## expr min lq mean median uq max neval cld
## DT[is.finite(rowSums(DT))] 21.617935 22.687452 26.698070 25.75765 26.07942 87.56290 100 c
## DT[complete.cases(DT * 0)] 7.209252 7.567393 9.908503 10.17569 10.37473 71.31375 100 a
## DT[DT[, Reduce("&", lapply(.SD, is.finite))]] 11.786773 12.647652 14.128624 14.78512 15.05089 15.39542 100 b
Run Code Online (Sandbox Code Playgroud)
| 归档时间: |
|
| 查看次数: |
10652 次 |
| 最近记录: |