これは私が取り組んでいるスクリプトです:
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 を使用しています。