【问题标题】:Lookup tables in RR中的查找表
【发布时间】:2021-02-11 22:09:31
【问题描述】:

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

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))

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

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

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

【问题讨论】:

标签: r tibble lookup-tables


【解决方案1】:

您也可以使用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

【讨论】:

    【解决方案2】:

    您可以使用 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)
    

    【讨论】:

      【解决方案3】:

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

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

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

      【讨论】:

        猜你喜欢
        • 2012-10-17
        • 1970-01-01
        • 1970-01-01
        • 2013-06-25
        • 2013-02-25
        • 1970-01-01
        • 2014-02-12
        • 1970-01-01
        • 2020-09-04
        相关资源
        最近更新 更多