apply() 函数中的错误，但 for 循环有效

Posted 2023-02-14

【中文标题】apply() 函数中的错误，但 for 循环有效

【英文标题】：Error in apply() function, for-loop works though

【发布时间】：2014-04-21 16:23:42

【问题描述】：

我正在尝试使用应用功能。我有一个名为 rowcosty() 的函数：

rowcosty <- function(p,x,y)
    return (-1*y*log(hyp(p,x))-(1-y)*log(1-hyp(p,x)))


hyp <- function(p, input)
    return (sigmoid(p %*% input))


sigmoid <- function(y)
    return (1/(1 + exp(-y)))

在我的示例中，p 是长度为 5000 的向量，x 是维度为 5000 X 400 的矩阵，y 也是长度为 5000 的向量。

现在我想将 rowcosty() 函数应用于 x 中的每一行，然后返回一个长度为 5000 的向量，其中包含 rowcosty() 为每一行返回的值。

我使用 for 循环成功地做到了这一点：

vec <- rep(0,5000)
for(i in 1:5000)
    vec[i] <- rowcosty(p,x[i,],y[i])

但是这个 apply() 函数的调用产生了错误：“Error in -y * log(hyp(p, x)) : non-conformable arrays”

apply(x,1,rowcosty,p = p, y = y)

谁能告诉我这里出了什么问题？

【参考方案1】：

当您执行rowcosty(p,x[i,],y[i]) 时，您同时间接地循环了两个变量：x 和y。 apply、lapply、sapply、vapply 都在单个变量上循环。

mapply 是循环两个或多个变量时的常用方法。你可以这样做：

mapply(rowcosty, split(x, row(x)), y, MoreArgs = list(p = p))

或者您也可以循环单个变量：索引i。

sapply(1:nrow(x), function(i) rowcosty(p,x[i,],y[i]))