0

我通常使用 和 的组合colwise来tapply计算数据框中的分组值。但是，我意外地发现FUNin的参数tapply无法与colwisefrom一起正常工作plyr。示例如下：

数据：

df <- data.frame(a = 1:10, b = rep(1:2, each = 5), c = 2:11)

普通的：

library(plyr)
colwise(tapply)(subset(df, select = c(a, c)), df$b, function(x){sum(x[x > 2])})

以上代码正确，可以正常工作。但是如果我添加FUN，那就错了：

colwise(tapply)(subset(df, select = c(a, c)), df$b, FUN = function(x){sum(x[x > 2])})

错误是：

Error in FUN(X[[1L]], ...) : 

unused arguments (function (X, INDEX, FUN = NULL, ..., simplify = TRUE) 
{
    FUN <- if (!is.null(FUN)) match.fun(FUN)
    if (!is.list(INDEX)) INDEX <- list(INDEX)
    nI <- length(INDEX)
    if (!nI) stop("'INDEX' is of length zero")
    namelist <- vector("list", nI)
    names(namelist) <- names(INDEX)
    extent <- integer(nI)
    nx <- length(X)
    one <- 1
    group <- rep.int(one, nx)
    ngroup <- one
    for (i in seq_along(INDEX)) {
        index <- as.factor(INDEX[[i]])
        if (length(index) != nx) stop("arguments must have same length")
        namelist[[i]] <- levels(index)
        extent[i] <- nlevels(index)
        group <- group + ngroup * (as.integer(index) - one)
        ngroup <- ngroup * nlevels(index)
    }
    if (is.null(FUN)) return(group)
    ans <- lapply(X = split(X, group), FUN = FUN, ...)
    index <- as.integer(names(ans))
    if (simplify && all(unlist(lapply(ans, length)) == 1)) {
        ansmat <- array(dim = extent, dimnames = namelist)

谁能解释原因？先感谢您。

4

1 回答 1

0

好吧，问题是两者lapply都有tapply一个可选FUN参数。请注意，这colwise(tapply)是一个具有以下行的函数：

out <- do.call("lapply", c(list(filtered, .fun, ...), dots))

让我们通过编写调试器转到这一行

ct <- colwise(tapply); trace(ct, quote(browser()), at = 6)

然后运行

ct(subset(df, select = c(a, c)), df$b, FUN = function(x){sum(x[x > 2])})

现在让我们打印c(list(filtered, .fun, ...), dots)。请注意，前三个（未命名的）参数现在是数据框、tapply和db$b，FUN上面的参数排在最后。然而，这个论点被命名为. 由于这是一个do.callon lapply，而不是那个参数成为一个可选参数tapply，它现在成为主要调用 on lapply！所以正在发生的事情是你正在把它变成：

lapply(subset(df, select = c(a, c)), function(x){sum(x[x > 2])}, tapply, df$b)

当然，这毫无意义，如果您手动执行上述操作（仍在调试器中），您将得到完全相同的错误。对于一个简单的解决方法，请尝试：

tapply2 <- function(.FUN, ...) tapply(FUN = .FUN, ...)
colwise(tapply2)(subset(df, select = c(a, c)), df$b, .FUN = function(x){sum(x[x > 2])})

该plyr包应该检查...命名的参数FUN（或任何可能干扰lapply' 工作的东西），但作者似乎没有包含这个。plyr您可以向实现以下任何变通方法的包提交拉取请求：

定义一个本地

.lapply <- function(`*X*`, `*FUN*`, ...) lapply(X = `*X*`, `*FUN*`, ...)

（进一步减少干扰）。

names(list(...))在colwise(tapply)函数中检查Xand FUN（如果作者打算在子调用之前阻止对承诺的评估，可能会引入问题）。

do.call("lapply", ...)使用命名X和显式调用FUN，以便您获得预期的

formal argument "FUN" matched by multiple actual arguments
于 2014-04-05T17:57:09.333 回答