从R中的data.table中删除带NA的行

Anj*_*jaM 7 r data.table

我想删除其任何列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 * 0NA与使用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)