R optim/nlm多次元配列のパラメータを持つ

Question

私はoptim/nlmを使って最尤推定を行い、パラメータは多次元配列になっています。

尤度は、データxと与えられたパラメータtheta、likelihood（theta、x）の多次元配列が与えられた場合に、実数を与えます。

しかし、OPTIM/NLMを使用して、シータだけで罰金を評価していたのと同じ寸法を有する開始値で、私は次のエラーを取得しています：

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

ときの評価見込み。 optim/nlmは私の多次元配列を1D配列に平坦化することが判明しました。とにかく、多次元配列のパラメータで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) 
    } 
)