4

行列のリストがあり、各行列の各列に関数を適用しようとしています。applyそれを行う自然な方法は、にネストすることだと思いますが、lapplyそうするとエラーが発生します。

library(hydroGOF)  # for the mse function

# Create a vector for comparison
A <- rnorm(10)

# Create a list of matrices to be looped over
B <- rnorm(10 * 5 * 3)
list.element.number <- rep(1:3, 50)
B.list <- split(B, list.element.number)
B.list <- lapply(B.list, matrix, ncol = 5, nrow = 10)  # A 3 element list of 10 x 5 matrices

# Wrapper function for mse
my.mse <- function(sim) {
  mse(sim, A)
}

# I'm trying to loop through each column of B.list and compare it to A
lapply(B.list, apply, MARGIN = 2, FUN = my.mse)
# Error in FUN(X[[1L]], ...) : 
#   unused arguments (function (X, MARGIN, FUN, ...) 
#   {
#     FUN <- match.fun(FUN)
#     dl <- length(dim(X))
#     if (!dl) stop("dim(X) must have a positive length")
#     if (is.object(X)) X <- if (dl == 2) as.matrix(X) else as.array(X)
#     d <- dim(X)
#     dn <- dimnames(X)
#     ds <- seq_len(dl)
#     if (is.character(MARGIN)) {
#       if (is.null(dnn <- names(dn))) stop("'X' must have named dimnames")
#       MARGIN <- match(MARGIN, dnn)
#       if (anyNA(MARGIN)) stop("not all elements of 'MARGIN' are names of dimensions")
#     }
#     s.call <- ds[-MARGIN]
#     s.ans <- ds[MARGIN]
#     d.call <- d[-MARGIN]
#     d.ans <- d[MARGIN]
#     dn.call <- dn[-MARGIN]
#     dn.ans <- dn[MARGIN]
#     d2 <- prod(d.ans)
#     if (d2 == 0) {
#       newX <- array(vector(typeof(X), 1), dim = c(prod(d.call), 1))
#       ans <- FUN(if (length(d.call) < 2) newX[, 1] else array(newX[, 1], d.call, dn.call), ...)
#       return(if (is.null(ans)) ans else if (length(d.ans) < 2) ans[1][-1] else array(ans, d.ans, dn.a

理想的には、これにより、各要素が 5 つの要素ベクトルである 3 つの要素のリストが得られますが、代わりにエラーが発生します。これを修正する方法を知っている人はいますか (洗練されていないと感じるループ以外for)、または何が問題なのですか?

4

1 に答える 1

4

FUN=my.mseへの試みをハイジャックしていますlapply apply。私はあなたが欲しいと思います：

lapply(B.list, function(x) apply(x, 2, my.mse))

FUN両方lapplyとapply同じ呼び出しで指定することはできません(Konrad が提案するように位置一致を使用しない限り)。何が起こっているかというと、あなたの呼び出しは次のように一致しています:

lapply(B.list, FUN=my.mse, ...=list(apply, MARGIN=2))

次に、最初の反復で評価を試みます。

my.mse(B.list[[1]], apply, MARGIN=2)

あなたの予想の代わりにapply(B.list[[1]], 2, my.mse)。1 つの引数しか受け付けないためmy.mse、「未使用の引数」に関するエラーが発生します。よく見ると「未使用の引数」がapply関数本体です。

于 2015-02-24T20:31:33.753 に答える