Bonjour,
J'essaie d'identifier des trajectoires d'évolution par le fonction jointlcmm. les modèles à une trajectoire tournent correctement avec une bonne convergence. Cependant, les modèles à plus d'une trajectoire ne tournent pas . voila le code (modèle à 2 trajectoires) et le message qui s'affiche après quelques minutes de calculs.
jlcmm(fixed = symp ~ ns(time, knots = 1, Boundary.knots = c(0,
1.8218)), mixture = ~ns(time, knots = 1, Boundary.knots = c(0,
1.8218)), random = ~ns(time, knots = 1, Boundary.knots = c(0,
1.8218)), subject = "id_ckdrein", ng = 2, survival = surv(Tevent,
event) ~ cause1(1) + cause2(1), hazard = c("Weibull", "Weibull"),B=mod1,
hazardtype = c("Specific", "Specific"), link = "beta", data = table_mcj,
maxiter = 200)
Statistical Model:
Dataset: table_mcj
Number of subjects: 2799
Number of observations: 9291
Number of latent classes: 2
Number of parameters: 24
Event 1:
Number of events: 1607
Class-specific hazards and
Weibull baseline risk function
Event 2:
Number of events: 702
Class-specific hazards and
Weibull baseline risk function
Link function for symp: Standardised Beta CdF
Iteration process:
The program stopped abnormally. No results can be displayed.
La fonction gridsearch aussi ne marche pas.
Dans la documentation du package, il est juste précisé que ça arrive si le paramètre conv du retournData est égale à 4 (conv==4) mais je ne vois pas la solution pour contourner ce problème. il faut préciser que ça tourne avec la fonction lcmm sans problème.
Merci pour votre aide.