$$$$$$$$

这是stata板块
你去R版本询问才合适

模型$pprob,例如m2$pprob

代码：m2$pprob

楼主，想咨询一个问题，用lcmm的过程中，ng>1的话不是还要确定初始值嘛，有3种方法，①B=方差协方差的参数，例如m2b <- hlme(normMMSE ~ age65+I(age65^2)+CEP, random =~ age65+I(age65^2), subject = 'ID',data = paquid, ng = 2, mixture =~ age65+I(age65^2), B = c(0, 60, 40, 0, -4, 0, -10, 10, 212.869397, -216.421323,456.229910, 55.713775, -145.715516, 59.351000, 10.072221))；②B=random（m1）;③用grid search ，例如m3g <- gridsearch(hlme(normMMSE ~ age65+I(age65^2)+CEP,  random =~ age65+I(age65^2), subject = 'ID', data=paquid, ng = 3, mixture=~age65+I(age65^2)), rep=100, maxiter=30, minit=m1)。方法③，我跑的时候总出现代码错误，你操作的过程中会这样吗？对于方法①，我一直不知道B值是怎么确定的，想请教一下