如何从循环中保存poLCA全局最大对数似然值数据？

针对你当前循环执行LCA的需求，不需要存储500次的全部结果，只需在循环中实时跟踪并更新全局最大的log-likelihood值即可，修改后的代码如下：

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

# 初始化最大log-likelihood为负无穷（因log-likelihood为负值）
max_llik <- -Inf
# 可选：若需保留对应最优结果的类别概率，初始化变量
best_class_probs <- NULL

## 循环执行LCA，仅保留全局最优结果
for (i in 1:500) {
  gss.lc <- poLCA(f, gss82, nclass=3, maxiter=3000, tol=1e-7)
  
  # 若当前结果的log-likelihood更大，则更新最大值
  if (gss.lc$llik > max_llik) {
    max_llik <- gss.lc$llik
    # 同步更新对应排序后的类别概率（和原逻辑一致）
    o <- order(gss.lc$probs$UNDERSTA[,1], decreasing=T)
    best_class_probs <- gss.lc$P[o]
  }
}

# 输出全局最大的log-likelihood
print(max_llik)
# 若需查看最优结果的类别概率，执行以下代码
# print(best_class_probs)

关键修改说明： #

• 移除了原代码中用于存储500次结果的mlmat矩阵，改用单个变量max_llik跟踪最大值，节省内存占用。
• 初始将max_llik设为负无穷，确保第一次循环的结果能顺利替换初始值。
• 循环内增加条件判断，仅当当前迭代的log-likelihood大于已记录的最大值时，才更新变量。
• 可选保留最优结果对应的类别概率，逻辑和原代码的排序规则一致。

内容的提问来源于stack exchange，提问作者Bruce