R data.table 根据 j 中的列动态过滤行

Question

我正在处理一个需要为少数预测变量生成 MAPE 的要求。为此，我使用 MAP 函数在 data.table 中生成带有 MAPE 值的汇总表。所以结果输出应该有 4 行 x 6 列，每个城市 1 行和 1 个城市列，后跟 a1、a2、a3、a4、a5，每个单元格中都有 MAPE 值。

以下是我使用的示例数据和代码（注意 - 考虑 a1,a2,a3....作为实际和 p1,p2,p3...作为预测）-

library(data.table)

set.seed(123)
id <- seq(1001,1100,1)
city <- sample(1:4,100,replace = T)
a1 <- sample(1:100,100,replace = T)
a2 <- sample(1:100,100,replace = T)
a3 <- sample(1:100,100,replace = T)
a4 <- sample(1:100,100,replace = T)
a5 <- sample(1:100,100,replace = T)
p1 <- sample(1:100,100,replace = T)
p2 <- sample(1:100,100,replace = T)
p3 <- sample(1:100,100,replace = T)
p4 <- sample(1:100,100,replace = T)
p5 <- sample(1:100,100,replace = T)

df1 <- as.data.table(data.frame(id,city,a1,a2,a3,a4,a5,p1,p2,p3,p4,p5))


sum1 <- df1[, Map(function(x,y) mean(as.numeric(abs(get(x)-get(y))/get(x))*100),
                                                                   paste("a",1:5, sep = ""),
                                                                   paste("p",1:5, sep = "")),by=city]

现在我想为 x==y 然后 x > y 和 x < y..... 的行生成相同的摘要，我认为最简单的方法是将它传递给 i 但是怎么做我没有得到.....当我尝试将它作为function(x,y) get(x)==get(y) 传递时它给出了错误

i 尚未计算为逻辑、整数或双精度

请推荐

Answer 1

更新答案：由于您希望每个 ==、> 和 < 有 3 个单独的数据表，因此我使用 lapply 来一一处理这些操作与您在帖子中发布的代码相同。

运算符（例如==）作为z 传递给Map。由于这是一个“运算符”，因此您不能使用get(z)，因此为了解决此问题，我使用do.call 将fn 函数（即== 运算符）应用于变量列表。

在Map 函数中，我首先根据fn 值（即z）将数据过滤为x1 和y1。然后对这些子集应用 MAPE 公式得到最终结果。

lapply(c("==", ">", "<"), 
       function(z) df1[, 
                       Map(function(x, y, fn = z){
                         x1 <- get(x)[do.call(fn, list(get(x), get(y)))]
                         y1 <- get(y)[do.call(fn, list(get(x), get(y)))]
                         mean(as.numeric(abs(x1 - y1) / x1) * 100)
                         },
                         paste0("a", 1:5),
                         paste0("p", 1:5),
                         z), 
                       by = city])