mj1<-Jointlcmm(SF_coxboxpoly(Time)+sex,random=poly(Time),survival=Surv(Tevent,event)sex,hazard="Weibull",subject='ID',data=dfR100,ng=1)
mj2<-Jointlcmm(SF_coxbox ~ poly(Time) + sex,mixture =~ poly(Time),random =~ poly(Time),survival = Surv(Tevent, event) ~ sex,hazard = "Weibull", subject = 'ID', data = dfR100, ng = 2, B = mj1)
mj3<-Jointlcmm(SF_coxbox ~ poly(Time) + sex,mixture =~ poly(Time),random =~ poly(Time),survival = Surv(Tevent, event) ~ sex,hazard = "Weibull", subject = 'ID', data = dfR100, ng = 3, B = mj1)
mj4<-Jointlcmm(SF_coxbox ~ poly(Time) + sex,mixture =~ poly(Time),random =~ poly(Time),survival = Surv(Tevent,event) ~ sex,hazard = "Weibull", subject = 'ID', data = dfR100, ng = 4, B = mj1)
得到
得到 G loglik npm BIC %class1 %class2 %class3 %class4
mj1 1 -4519.358 10 9094.589 100.00000
mj2 2 -4420.839 15 8925.487 72.65918 27.34082
mj3 3 -4420.839 20 8953.423 0.00000 72.65918 27.34082
mj4 4 -4420.845 25 8981.372 0.00000 72.65918 27.34082 0
想改进B值应该怎么自己设定啊,参考什么?哪位老师帮帮忙,万分感谢!