2

これは私が取り組んでいるスクリプトです:

library(poLCA)

f <- cbind(bq70, bq72_1, bq72_2, bq72_3, bq72_4, bq72_5, 
           bq72_6, bq72_7, bq73a_1, bq73a_2, bq73a_3, bq73a_4) ~ 
     zq88 + zq89 + dm_zq101_2 + dm_zq101_3 + dm_zq101_4 + 
     dm_zq101_5 + dm_zq101_6 + dm_zq101_7 + dm_zq101_8 + dm_zq101_9

for(i in 2:14){
  max_II <- -1000000
  min_bic <- 100000

  for(j in 1:1024){
    res <- poLCA(f, BESDATA, nclass=i, maxiter=1000, 
                 tol=1e-5, na.rm=FALSE, probs.start=NULL, 
                 nrep=1, verbose=TRUE, calc.se=TRUE)
    if(res$bic < min_bic){
      min_bic <- res$bic
      LCA_best_model<-res
    }
  }
}

潜在クラス分析と回帰を実行したいと思います。ただし、上記のコードを実行すると、PC が完了するまでに非常に長い時間がかかります (インテル Core i5 4690k、16GB RAM)。

PoLCA でこれほど時間がかかるのは一般的ですか。

また、グローバルな最大可能性に達したら、各クラスのループを停止するために使用できるコード行はありますか?

N = 2000 前後。

ちなみに、私は R studio を使用しています。

4

1 に答える 1