具有多维参数数组的 R optim/nlm

Question

我正在使用 optim/nlm 进行最大似然估计，并且我的参数位于多维数组中。

似然评估良好，即给定数据 x 和参数 theta 的多维数组，likelihood(theta,x) 给出实数。

然而，使用 optim/nlm，起始值与 theta 的维度相同，评估得很好，我收到以下错误：

Error in theta[1, 1, 1] : incorrect number of dimensions

评估可能性时。事实证明 optim/nlm 将我的多维数组展平为一维数组。无论如何，我可以将 optim/nlm 与多维参数数组一起使用吗？

Answer 1

我不相信 optim 本身可以实现这一点。我的建议是自己恢复形状，例如

optim(
    matrix(1:4, 2, 2),
    function(par) {
        par = matrix(par, 2, 2) # Reshape
        sum((par - matrix(5:8, 2, 2))**2)
    }
)