r - R 最尤アルゴリズム

Question

最尤推定量を最初から作成する必要があり、組み込みの mle 関数は許可されていません。私はランダム変数を最大化していますが、平均や分散についてではなく、別のパラメーター (移動平均 (1) プロセスからのベータ) を持っています。ここに私が作成したコードがありますが、「最適化」関数が初期値を見つけることができないことがわかります。私は行列に観測を持っています。行は異なる観測値であり、列は異なるシミュレーションです (そのため、最初に 1 から R=10000 まで実行するループがあります)。大いに感謝します:)

for (i in 1:R) {
llh_regr <- function (β, y, z) {
  kappa <- -((N-1)*log(pi*2)+(N-1)*log(1))/2
  residuals <- (y-c-β*z)
  loglik <- kappa-(1/2)*crossprod(residuals, residuals)
  return(loglik)
                               }

targetfunction <- function(β) {
  return(-llh_regr(β, y = x[,1], z = u_lag[,1]))
         }
  B[1,i] <- optim(runif(1), targetfunction)
}