如何在data.table中通过引用更改每个组中的最后一个值

For a data.table DT grouped by site, sorted by time t, I need to change the last value of a variable in each group. I assume it should be possible to do this by reference using :=, but I haven't found a way that works yet.

对于按站点分组的data.table DT,按时间t排序,我需要更改每个组中变量的最后一个值。我认为应该可以通过引用使用:=来做到这一点,但我还没有找到一种方法。

Sample data:

require(data.table)   # using 1.8.11 
DT <- data.table(site=c(rep("A",5), rep("B",4)),t=c(1:5,1:4),a=as.double(c(11:15,21:24)))
setkey(DT, site, t)
DT
#    site t  a
# 1:    A 1 11
# 2:    A 2 12
# 3:    A 3 13
# 4:    A 4 14
# 5:    A 5 15
# 6:    B 1 21
# 7:    B 2 22
# 8:    B 3 23
# 9:    B 4 24

The desired result is to change the last value of a in each group, for example to 999, so the result looks like:

所需的结果是更改每个组中a的最后一个值,例如更改为999,因此结果如下所示:

#    site t   a
# 1:    A 1  11
# 2:    A 2  12
# 3:    A 3  13
# 4:    A 4  14
# 5:    A 5 999
# 6:    B 1  21
# 7:    B 2  22
# 8:    B 3  23
# 9:    B 4 999

It seems like .I and/or .N should be used, but I haven't found a form that works. The use of := in the same statement as .I[.N] gives an error. The following gives me the row numbers where the assignment is to be made:

似乎应该使用.I和/或.N,但我还没有找到一个有效的表单。在与.I [.N]相同的语句中使用:=会产生错误。以下给出了要进行赋值的行号:

DT[, .I[.N], by=site]
#    site V1
# 1:    A  5
# 2:    B  9

but I don't seem to be able to use this with a := assignment. The following give errors:

但我似乎无法使用:=赋值。以下给出了错误:

DT[.N, a:=999, by=site]
# Null data.table (0 rows and 0 cols)

DT[, .I[.N, a:=999], by=site]
# Error in `:=`(a, 999) : 
#   := and `:=`(...) are defined for use in j, once only and in particular ways.
#  See help(":="). Check is.data.table(DT) is TRUE.

DT[.I[.N], a:=999, by=site]
# Null data.table (0 rows and 0 cols)

Is there a way to do this by reference in data.table? Or is this better done another way in R?

有没有办法通过data.table中的引用来做到这一点?或者在R中用另一种方法做得更好?

1 个解决方案

#1

Currently you can use

目前你可以使用

DT[DT[,.I[.N],by=site][['V1']],a:=999]
# or, avoiding the overhead of a second call to `[.data.table`
set(DT, i = DT[,.I[.N],by='site'][['V1']], j = 'a', value = 999L)

alternative approaches:

use replace...

DT[, a := replace(a,seq_len(.N)==.N,999) ,by=site]

or shift the replacement to the RHS, wrapped by {} and return the full vector

或者将替换转移到RHS,由{}包裹并返回完整的向量

DT[,a:={a[.N]<-999L; a},by=site]

or use mult='last' and take advantage of by-without-by. This requires the data.table to keyed by the groups of interest.

或使用mult ='last'并利用by-without-by。这要求data.table由感兴趣的组键入。

 DT[unique(site),a:=999,mult='last']

There is a feature request #2793 that would allow

有一个功能请求#2793允许

DT[, a[.N] := 999]

but this is yet to be implemented

但这还没有实施

#1