如何保持 LCA 模型的最大对数似然？

Question

我想为我的 LCA 模型找到更好的对数似然，但我太确定对数似然是尽可能小，还是尽可能重复？ 问题来自 poLCA 论文 (http://www.sscnet.ucla.edu/polisci/faculty/lewis/pdf/poLCA-JSS-final.pdf p14.table1)

我正在复制代码尝试测试

## Example data
library(poLCA)
data(gss82)
f <- cbind(PURPOSE,ACCURACY,UNDERSTA,COOPERAT)~1

mlmat <- matrix(NA,nrow=500,ncol=4)

## Do repreat 500 times LCA for best llik
for (i in 1:500) {
  gss.lc <- poLCA(f,gss82,nclass=3,maxiter=3000,tol=1e-7)
  mlmat[i,1] <- gss.lc$llik
  o <- order(gss.lc$probs$UNDERSTA[,1],decreasing=T)
  mlmat[i,-1] <- gss.lc$P[o]
 }

此外，我想知道如何在我的数据中存储对数似然的不同最大值，如 lc_llik_2754.5、lc_llik_27556、lc_llik_2755.7 ...等

谢谢！