r - optim() 并解决 R 中的对角矩阵（拟合参数错误）

Question

我正在尝试在以下问题中找到解决方案：

可重复性的任意值：

AZ2k <- seq(1:14) 
noAZ2k <- matrix(seq(1:25), 14, 50)
par <- rep(0.02 , 50)


vw <- function(w)
{

    t(AZ2k - noAZ2k %*% par)%*%w%*%(AZ2k - noAZ2k %*% par ) 

}


vweights <- optim(diag(1,14), vw, gr = NULL,
  method = c("L-BFGS-B"),
  lower = 0, upper = 10000,control = list( factr = 1e4, maxit = 1000, pgtol = .01 ), 
  hessian = FALSE)

当我输入

t(AZ2k - noAZ2k %*% par)%*%diag(1,14)%*%(AZ2k - noAZ2k %*% par )

我得到了一个结果，但是当我尝试运行优化时，它说这些值不合适，这让我感到惊讶。

我可能遗漏了一些完全明显的东西，但我就是找不到哪里出错了，除非 optim 只是使用错误的函数，但我无法找到合适的替代方案。

Answer 1

您的 w 参数将转换为 vector ，因此您需要将其强制转换为具有正确尺寸的矩阵：

vw <- function(w){
    w <- matrix(w,nrow=14,ncol=14,byrow=T)
    t(AZ2k - noAZ2k %*% par)%*%w%*%(AZ2k - noAZ2k %*% par ) 
}