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
以不同的方式编写它,但它最终会变成更多代码。
我有一个包含大量数据的小标题,但最重要的是,我有一列按数字(例如 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
以不同的方式编写它,但它最终会变成更多代码。