R中的查找表

Lookup tables in R

我有一个包含大量数据的小标题,但最重要的是,我有一列按数字(例如 1、2、3 等)引用查找 table 中的一行。

df <- tibble(ref = c(1,1,1,2,5)
             data = c(33,34,35,35,32))

lkup <- tibble(CurveID <- c(1,2,3,4,5)
               Slope <- c(-3.8,-3.5,-3.1,-3.3,-3.3)
               Intercept <- c(40,38,40,38,36)
               Min <- c(25,25,21,21,18)
               Max <- c(36,36,38,37,32))

我需要根据查找中引用行中的信息对原始标题中的每一行进行计算 table

df$result <- df$data - lkup$intercept[lkup$CurveID == df$ref]/lkup$slope[lkup$CurveID == df$ref]

想法是根据数据 table 中的数字从正确的查找行 table 中访问斜率或截距等值,并对中的每个数据点执行此操作专栏。但我一直收到错误消息,告诉我我的数据不兼容,并且我的 objects 需要具有相同的长度。

你也可以用 match()

df$result <- df$data - lkup$Intercept[match(df$ref, lkup$CurveID)]/lkup$Slope[match(df$ref, lkup$CurveID)]
df$result
# [1] 43.52632 44.52632 45.52632 45.85714 42.90909

您可以使用 dplyr 包将 tibbles 连接在一起。如果 ref 列和 CurveID 列具有相同的名称,那么 left_join 将通过匹配的行组合两个小标题。

library(dplyr)
df <- tibble(CurveID = c(1,1,1,2,5),
             data = c(33,34,35,35,32))

lkup <- tibble(CurveID = c(1,2,3,4,5),
               Slope = c(-3.8,-3.5,-3.1,-3.3,-3.3),
               Intercept = c(40,38,40,38,36),
               Min = c(25,25,21,21,18),
               Max = c(36,36,38,37,32))

df <- df %>% left_join(lkup, by = "CurveID")

然后对每一行进行计算

df <- df %>% mutate(result = data - (Intercept/Slope)) %>% 
      select(CurveID, data, result)

为了完整起见,这里有一种方法可以真正做到 OP 正在尝试的事情:

library(slider)
df %>%
  mutate(result = slide_dbl(ref, ~ slice(lkup, .x)$Intercept /
                                   slice(lkup, .x)$Slope))

尽管由于 slice 按行号排列,这依赖于 CurveID 等于行号(我们根本没有引用 CurveID)。您可以使用 filter 以不同的方式编写它,但它最终会变成更多代码。