(X'X)^-1 の決定を返す関数を作成しました。X 行列の最適な値を得るために、作成した関数で最適関数を使用する必要があります。最適な動作に問題がありますが、機能は正常に動作しています。
fr <- function(x,nc,nr) {
Xd <- matrix(c(rep(1,nr),x),nrow=nr,ncol=nc) #### Design matrix
det(solve((t(Xd)%*%Xd))) #### det(FIM)
}
x をベクトルと行列の両方として試しました。
x <-
matrix(c(0.1,0.2,0.3,0.2,0.1,0,0.1,0.2,0.1,0.2,0.1,0.2,0.1,0.2,0.1,0.2,0.1,0.2,0.1,0.2),nro w = 10, ncol = 2)
x <- c(0.1,0.2,0.3,0.2,0.1,0,0.1,0.2,0.1,0.2,0.1,0.2,0.1,0.2,0.1,0.2,0.1,0.2,0.1,0.2)
これらは、私が optim で試みたいくつかのバリエーションですが、すべて成功しませんでした。
optim((x, 3, 10), fr)
optim(fn = fr, par = c(x=x, nc=3, nr=10))
optim(c(x,3,10), fr)
optim の x の後にコンマがある場合や、.Internal で nrow が見つからないというエラーが発生する場合があります。