[英]R data.table Retain the first non-NA value of a group till end of the group
[英]R: find first non-NA observation in data.table column by group
我有一个带有许多缺失值的data.table
,我想要一个变量,它给出了每个组中第一个非data.table
值的1。
说我有这样的数据。表:
library(data.table)
DT <- data.table(iris)[,.(Petal.Width,Species)]
DT[c(1:10,15,45:50,51:70,101:134),Petal.Width:=NA]
现在在开头,结尾和中间都有缺失。 我试过两个版本,一个是:
DT[min(which(!is.na(Petal.Width))),first_available:=1,by=Species]
但它只能找到全局最小值(在这种情况下,setosa得到正确的1),而不是组的最小值。 我认为是这种情况,因为data.table
第一个子集由i,然后按组排序,对吗? 所以它只适用于全局最小值的行which(!is.na(Petal.Width))
,这是第一个非NA值。
j中测试的第二次尝试:
DT[,first_available:= ifelse(min(which(!is.na(Petal.Width))),1,0),by=Species]
它只返回一列1。 在这里,我没有很好的解释为什么它不起作用。
我的目标是:
DT[,first_available:=0]
DT[c(11,71,135),first_available:=1]
但实际上我有数百个团体。 任何帮助,将不胜感激!
编辑: 这个问题确实接近但不是针对NA的,如果我理解正确的话,并没有解决这个问题。 我试过了:
DT <- data.table(DT, key = c('Species'))
DT[unique(DT[,key(DT), with = FALSE]), mult = 'first']
这是一种方式:
DT[!is.na(Petal.Width), first := as.integer(seq_len(.N) == 1L), by = Species]
我们可以尝试
DT[DT[, .I[which.max(!is.na(Petal.Width))] , Species]$V1,
first_available := 1][is.na(first_available), first_available := 0]
或者是一个稍微紧凑的选项
DT[, first_available := as.integer(1:nrow(DT) %in%
DT[, .I[!is.na(Petal.Width)][1L], by = Species]$V1)][]
> DT[!is.na(DT$Petal.Width) & DT$first_available == 1]
# Petal.Width Species first_available
# 1: 0.2 setosa 1
# 2: 1.8 versicolor 1
# 3: 1.4 virginica 1
> rownames(DT)[!is.na(DT$Petal.Width) & DT$first_available == 1]
# [1] "11" "71" "135"
> rownames(DT)[!is.na(DT$Petal.Width) & DT$first_available == 0]
# [1] "12" "13" "14" "16" "17" "18" "19" "20" "21" "22" "23" "24"
# [13] "25" "26" "27" "28" "29" "30" "31" "32" "33" "34" "35" "36"
# [25] "37" "38" "39" "40" "41" "42" "43" "44" "72" "73" "74" "75"
# [37] "76" "77" "78" "79" "80" "81" "82" "83" "84" "85" "86" "87"
# [49] "88" "89" "90" "91" "92" "93" "94" "95" "96" "97" "98" "99"
# [61] "100" "136" "137" "138" "139" "140" "141" "142" "143" "144" "145" "146"
# [73] "147" "148" "149" "150"
声明:本站的技术帖子网页,遵循CC BY-SA 4.0协议,如果您需要转载,请注明本站网址或者原文地址。任何问题请咨询:yoyou2525@163.com.