我在R(poLCA)做潜类分析,现在我在作业中有一个问题,我必须对我的潜类模型施加一些限制,例如。
P_ {11} ^ B | X = P_ {11} ^ C | X
或
P_ {11} ^ {A | X} = 1
其中X是潜类,A,B,C是独立的。我试图弄乱起始值(probs.start
),但这没有帮助。
somemodel <- poLCA(cbind(A, B, C, D) ~ 1, df1_expanded, nclass = 2)
有什么想法吗?
更新: 在其他基于脚本的程序(例如Stata)中是否存在具有良好(乳胶)表输出选项的此问题的解决方案?