【发布时间】:2020-11-14 03:11:29
【问题描述】:
我想为每个组替换我的 data.table 中的缺失值,并根据组中的所有值是否缺失,或者只是组中的一些值缺失来填写值。
我可以解决问题,但愿意接受更好的代码(在速度/内存/可读性/灵活性方面)。
我很固执,我更喜欢 data.table 解决方案:)
示例:
起始数据集:
它是一个具有这种结构的data.table:
dt = data.table(
grouping_1 = sort(rep(c('a', 'b', 'c'), 4)),
grouping_2 = c(1,1,2,2,1,1,2,2,1,1,2,2),
value_1 = c(NA, NA, NA, NA, NA, 1, 2, NA, 3, 2,4,NA),
value_2 = c(NA, 2, NA, NA, 2, 5, 2, 7, 10, 5,NA, NA)
)
看起来像这样:
grouping_1 grouping_2 value_1 value_2
1: a 1 NA NA
2: a 1 NA 2
3: a 2 NA NA
4: a 2 NA NA
5: b 1 NA 2
6: b 1 1 5
7: b 2 2 2
8: b 2 NA 7
9: c 1 3 10
10: c 1 2 5
11: c 2 4 NA
12: c 2 NA NA
我想做的是:
我想按列 grouping_1 和 grouping_2 对其进行分组,并替换列 value_1 和 value_2 中的缺失值。
如果给定组没有非缺失值(例如组grrouping_1==a & grouping_2==1),我想用值 9000 替换该组的所有 NA。
如果给定组有一些非缺失值,我想将缺失值替换为 800 if grouping_2==1 并替换为 -800(负 800)if grouping_2==2。如果该值没有丢失,我不想更改它。
我现在该怎么做:
我编写了以下函数,然后将其应用于我要填写缺失值的每一列。该函数通过引用更改原始数据集:
filler_so = function(
data, # the dataset that we will be changing
column, # the column we will be filling in
placeholder_col ='drop_at_the_end', # some temporary column that will disappear in the end
missing_fully = 9000, # value to fill in when all values in group missing
missing_partially_g2_1 = 800, # value to fill when grouping_2 = 1
missing_partially_g2_2 = -800, # value to fill when grouping_2 = 2
g2_col = 'grouping_2', # name of column corresponding to grouping_2 from my example
group_cols = c('grouping_1', 'grouping_2') # names of columns to group by
){
# identify for given column whether all values in group are missing,
# or only some are misisng. The value will be either Infinity (all missig),
# or a real number (none or some missing).
# this info is put in a placeholder column
data[, (placeholder_col) := min(get(column), na.rm = T), by = group_cols]
# if value on a given row is missing, but not all missing in group,
# then fill in the values based on what group is in 2nd grouping column
data[
is.na(get(column)) & (get(placeholder_col) != Inf),
(placeholder_col) := (get(g2_col) == 2) * missing_partially_g2_2 +
(get(g2_col) ==1) * missing_partially_g2_1]
# if all values in group are missing, fill in the "missing_fully" value
data[get(placeholder_col) == Inf, (placeholder_col) := missing_fully]
# put into placeholder column the values that were originally not missing
data[!is.na(get(column)), (placeholder_col) := get(column)]
# drop the original column
data[, (column):=NULL]
# rename the placeholder column to the name of original column
setnames(data, placeholder_col, column)
# if i don't put this here,
# then sometimes the function doesn't return results properly.
# i have no clue why.
data
}
要应用此功能,我需要确定要填充的列,我这样做:
cols_to_fill = colnames(dt)[grep('^value', colnames(dt))]
然后像这样 lapply:
lapply(cols_to_fill, function(x) filler_so(dt, x))
结果:
> dt
grouping_1 grouping_2 value_1 value_2
1: a 1 9000 800
2: a 1 9000 2
3: a 2 9000 9000
4: a 2 9000 9000
5: b 1 800 2
6: b 1 1 5
7: b 2 2 2
8: b 2 -800 7
9: c 1 3 10
10: c 1 2 5
11: c 2 4 9000
12: c 2 -800 9000
我想改进的地方:
- 我的函数可以工作,但很冗长,我希望我可以将代码变成更少的行
- 该函数如果不是很灵活 - 最好传递一个命名向量之类的东西来指定替换逻辑和要填充的值基于
grouping_2 - 我正在寻找速度和内存增益。 (例如,可能有更快的方法来识别缺少所有值的组,然后运行 `min(..., na.rm = TRUE) 然后检查它何时为 Infinity。
- lapply 为我正在填写的每一列打印出更改后的 data.table,这会给控制台带来相当多的垃圾邮件。
- 即使 4 处的问题已解决,我想知道是否有办法代替 lapply 我可以使用
dt[..., (some_column_names) := lapply(.SD, ...), .SDcols = cols_to_fill]之类的东西来做到这一点 - 我愿意接受我没有想到的任何其他改进。
【问题讨论】:
标签: r data.table