Bonjour chers membres,

J'aimerais avoir votre aide. J'estime un modèle paramétrique de survie à population séparée (split population duration) sur R. et J'ai des soucis, mon code ne fonctionne pas
Si quelqu'un peut m'aider
Code : Sélectionner tout - Visualiser dans une fenêtre à part
1
2
3
4
5
6
7
8
9
10
 
dur.coups <- add_duration(mydata, survie, unitID="rge_id", tID=time,
                          freq="year")
model.coups <- spdur(duree ~ X, survie ~ X, data = dur.coups, distr = "loglog",
                    silent = TRUE)
summary(model.coups)
 
spdur(duree, survie, mydata, last = "64", t.0 = "0",
      fail = "survie", distr = c("weibull", "loglog"), max.iter = 300,
      na.action, silent = FALSE, ...)
j'ai essayer aussi sur stata avec la commande spsurv mais ça ne marche pas.