我有几个数据框架看起来类似于以下数据框架(有更多的列):
id col1 col2 col3 col4 col5
1 4 3 5 4 A
1 3 5 4 9 Z
1 5 8 3 4 H
2 6 9 2 1 B
2 4 9 5 4 K
3 2 1 7 5 J
3 5 8 4 3 B
3 6 4 3 9 C
我想计算按id分组的特定列的标准差(比如col2到col4)。我不知道每个数据帧中的列索引。我只知道要计算标准差的列的名称。
有什么办法我可以很容易做到吗?我的原始数据框架包含大约20列,我只想要10列的标准偏差,这些列有特定的列名,按id分组。
最重要的是,如果我能够根据id将计算出来的标准差直接添加到我的数据框架中作为一个新列,那就太好了,如下所示:
id col1 col2 col3 col4 col5 SD
1 4 3 5 4 A SD1
1 3 5 4 9 Z SD1
1 5 8 3 4 H SD1
2 6 9 2 1 B SD2
2 4 9 5 4 K SD2
3 2 1 7 5 J SD3
3 5 8 4 3 B SD3
3 6 4 3 9 C SD3
发布于 2021-01-25 05:20:09
你可以试试:
library(dplyr)
df %>%
group_by(id) %>%
mutate(SD = sd(unlist(select(cur_data(), col2:col4))))
# id col1 col2 col3 col4 col5 SD
# <int> <int> <int> <int> <int> <chr> <dbl>
#1 1 4 3 5 4 A 2.12
#2 1 3 5 4 9 Z 2.12
#3 1 5 8 3 4 H 2.12
#4 2 6 9 2 1 B 3.41
#5 2 4 9 5 4 K 3.41
#6 3 2 1 7 5 J 2.62
#7 3 5 8 4 3 B 2.62
#8 3 6 4 3 9 C 2.62
发布于 2021-01-25 09:45:05
使用data.table
library(data.table)
setDT(df)[, SD := sd(unlist(.SD)), id, .SDcols = col2:col4]
https://stackoverflow.com/questions/65885374
复制相似问题