我在data.frame
中有一系列列,其中我想获取最后一个值,不包括任何NA。我用来完成这项工作的功能是
last_value <- function(x) tail(x[!is.na(x)], 1)
我正在使用apply()
在13列中执行此功能,每次观察(按行)。
df$LastVal<-apply(df[,c(116, 561, 1006, 1451, 1896, 2341, 2786, 3231,
3676, 4121, 4566, 5011, 5456)], 1, FUN=last_value)
我的问题是输出是以5336(总观察值)的列表形式出现的,而不是仅按行的最后一个值的向量。答案似乎在那里,但再次,以列表形式。我之前使用过这个功能,它运行良好。当我str()
我的列时,它们都是整数。
如果没有值且只有NA,这个函数会被绊倒吗?
我应该补充一点,当我unlist()
新变量时,我收到的错误是“替换有4649行,数据有5336”,所以我认为这可能与NA有关。
答案 0 :(得分:0)
您可以将您的功能包含在选择中。
<强> 实施例 强>
df <- as.data.frame(matrix(1:12, 3, 4))
> df
V1 V2 V3 V4
1 1 4 7 10
2 2 5 8 11
3 3 6 9 12
last_value <- function(x) tail(x[!is.na(x)], 1)
> df[, last_value(c(3, 4))] # selects last column
[1] 10 11 12
使用NA
进行测试。
df[2, 4] <- NA
> df[, last_value(c(3, 4))]
[1] 10 NA 12
如果你需要apply()
方法,可以使用@Rui Barradas解释得很清楚。如果您依赖速度,请考虑两种解决方案的基准:
Unit: microseconds
expr min lq mean median uq max neval cld
apply(df, 1, last_value) 166.095 172.6005 182.09241 177.449 188.2925 257.179 100 b
df[, last_value(c(3, 4))] 32.147 33.4230 36.12764 34.699 35.5920 131.396 100 a
Apropos-用于列式使用sapply()
。
> sapply(df[, c(3, 4)], FUN=last_value)
V3 V4
9 12
答案 1 :(得分:0)
首先,您需要查看函数last_value
的输出,因为您已使用一行NA
值对其进行了定义。
last_value <- function(x) tail(x[!is.na(x)], 1)
df <- matrix(1:24, 4)
df[2, ] <- NA
df <- as.data.frame(df)
apply(df, 1, last_value)
#[[1]]
#V6
#21
#
#[[2]]
#named integer(0)
#
#[[3]]
#V6
#23
#
#[[4]]
#V6
#24
问题是此列表的第二个成员长度为零。这意味着unlist
无法解决问题
您必须测试长度为零的值。
last_value <- function(x) {
y <- tail(x[!is.na(x)], 1)
if(length(y) == 0) NA else y
}
apply(df, 1, last_value)
#[1] 21 NA 23 24