R 다차원 배열의 매개 변수가있는 optim/nlm

Question

나는 maximum likelihood estimation을하기 위해 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) 
    } 
)