使用多个条件进行匹配以在R中生成值

问题描述 投票:0回答:1

我目前具有以下数据格式:

df = data.frame(c(rep("A", 12), rep("B", 12)), rep(1:12, 2), seq(-12, 11))
colnames(df) = c("station", "month", "mean")
df

df_master = data.frame(c(rep("A", 10), rep("B", 10)), rep(c(27:31, 1:5), 2), rep(c(rep(1, 5), rep(2, 5)), 2), rep(seq(-4,5), 2))
colnames(df_master) = c("station", "day", "month", "value")
df_master

有效地,df是每个站点的月平均值,我想在df_master数据集中计算一个新变量,该变量计算每个每日观测值与月平均值的差。我设法通过包含所有数据的总体平均值来做到这一点,但由于每个站点的平均值各不相同,因此我想特定于新的可变站点。

我已经尝试使用以下代码来匹配每月的价值,但是目前尚不能解决跨站差异:

df_master$mean = df$mean[match(df_master$month, df$month)]
df_master = df_master %>% mutate(diff = value - mean)

我该如何进一步进行以使每个站点获得平均值?

r dplyr match mean difference
1个回答
1
投票

如果将它们转换为data.tables,则可以通过更新联接添加差异列,并在df_masterdf的值相等的情况下将stationmonth联接。

library(data.table)
setDT(df_master)
setDT(df)

df_master[df, on = .(station, month), 
          diff_monthmean := value - i.mean]

df_master
#     station day month value diff_monthmean
#  1:       A  27     1    -4              8
#  2:       A  28     1    -3              9
#  3:       A  29     1    -2             10
#  4:       A  30     1    -1             11
#  5:       A  31     1     0             12
#  6:       A   1     2     1             12
#  7:       A   2     2     2             13
#  8:       A   3     2     3             14
#  9:       A   4     2     4             15
# 10:       A   5     2     5             16
# 11:       B  27     1    -4             -4
# 12:       B  28     1    -3             -3
# 13:       B  29     1    -2             -2
# 14:       B  30     1    -1             -1
# 15:       B  31     1     0              0
# 16:       B   1     2     1              0
# 17:       B   2     2     2              1
# 18:       B   3     2     3              2
# 19:       B   4     2     4              3
# 20:       B   5     2     5              4

1
投票

使用左连接的dplyr

library(dplyr)
left_join(df_master, df, by = c('station', 'month')) %>% 
        mutate(monthdiff  = value - mean) %>%
        select(-mean)
© www.soinside.com 2019 - 2024. All rights reserved.