我有这类数据:
df <- data.frame(
A_aoi = c("C*BB*B", "C*B*C*", "B**", "C*B"),
A_dur = c("234,312,222,3456,1112,77", "12,13,14,15,11,1654", "896,45222,55", "5554,322,142"),
B_aoi = c("**ACC", "AC*", "AAA", "C*A*"),
B_dur =c("12,13,15,100,100", "14,55,66", "88,99,100", "1,2,3,4")
)我需要做的是在A_dur和B_dur列中总结那些值,这些值在字符串中的位置与A_aoi和B_aoi列中的运行长度重复值相同。
现在‘正电子’可以被索引捕获,这就是为什么我首先在所有列中拆分字符串的原因:
library(stringr)
df[,c(1,3)] <- lapply(df[,c(1,3)], function(x) str_split(x, ""))
df[,c(2,4)] <- lapply(df[,c(2,4)], function(x) str_split(x, ","))然后,在A_aoi和B_aoi列中得到了运行长度重复的索引。
library(dplyr)
df$index_A_aoi <- sapply(df$A_aoi, function(x) which(x == lead(x)))
df$index_B_aoi <- sapply(df$B_aoi, function(x) which(x == lead(x)))现在,当我试图根据A_dur和index_B_aoi中的索引和sapply来总结index_A_aoi和index_B_aoi中的值时,我陷入了困境:
sapply(as.numeric(df$A_dur), function(x) x[df$index_A_aoi] + lead(x[df$index_A_aoi]))
Error in lapply(X = X, FUN = FUN, ...) :
'list' object cannot be coerced to type 'double'预期结果
df
A_aoi A_dur B_aoi B_dur
1 C*BB*B 234,312,3678,1112,77 **ACC 25,15,200
2 C*B*C* 12,13,14,15,11,1654 AC* 14,55,66
3 B** 896,45277 AAA 287
4 C*B 5554,322,142 C*A* 1,2,3,4我对包括dplyr解决方案在内的任何解决方案持开放态度,并对此表示感谢。
发布于 2021-03-11 09:21:18
下面是一个借助data.table的rleid函数的解决方案
library(data.table)
calculate <- function(p, q) {
mapply(function(x, y) toString(tapply(as.numeric(x), rleid(y), sum)),
strsplit(p, ','), strsplit(q, ''))
}
aoi_cols <- grep('aoi', names(df))
dur_cols <- grep('dur', names(df))
df[dur_cols] <- Map(calculate, df[dur_cols], df[aoi_cols])
df
# A_aoi A_dur B_aoi B_dur
#1 C*BB*B 234, 312, 3678, 1112, 77 **ACC 25, 15, 200
#2 C*B*C* 12, 13, 14, 15, 11, 1654 AC* 14, 55, 66
#3 B** 896, 45277 AAA 287
#4 C*B 5554, 322, 142 C*A* 1, 2, 3, 4calculate函数接受一对列,在逗号上拆分dur,在每个字符上拆分aoi,并使用rleid创建值连续出现的组,并对它们进行sum。我们对数据集中的所有这类对应用calculate函数。
https://stackoverflow.com/questions/66579443
复制相似问题