Optim()使参数在R中相等

M = [ var1    corr    0]
    [ corr    var2    0]
    [  0       0      0]

基本上,矩阵在M[1,2]和M[2,1]位置有两个方差和一个相关性,这两个位置应该相等。为此,我为optim编写了一个函数,如下所示:

M <- matrix(c(NA,NA,0,NA,NA,0,0,0,0),3,3,byrow=TRUE)

objf <- function(pars, model, estimate = TRUE) {

    model$M[is.na(model$M)] <- pars[1:4]

    if (estimate) {
        -logLik(model)
    } else {
         model
    }
}

但是,当我运行代码时,两个相关值的输出是不同的,因此我想知道是否有更好的方法来编写代码,使我的M[1,2]=M[2,1]。我想用电话 model$M[is.na(model$M)]