Bonjour,
J'ai quelques difficultés à bien utiliser la fonction Surv du pack survival. J'essaie de déterminer les paramètres d'une distribution de Weibull grâce à la fonction survreg pour déterminer l'influence de variables sur un taux de casses.
J'ai un tableau de données qui me donne les taux de casses quotidiens.
Problème : la moitié de mes taux sont nuls donc il faut censurer les valeurs nulles à gauche. J'ai donc créé une variable censure (respectivement une variable event) qui affecte le taux initial (resp 1) s'il n'est pas nul ou le minimum des taux (resp 2) s'il est nul.
extrait du tableau de données :
tauxcompletA censure event
1053 7.177240e-07 7.177240e-07 1
1054 1.435448e-06 1.435448e-06 1
1055 1.435448e-06 1.435448e-06 1
1056 1.435448e-06 1.435448e-06 1
1057 1.435448e-06 1.435448e-06 1
1058 7.177240e-07 7.177240e-07 1
1059 7.177240e-07 7.177240e-07 1
1060 7.177240e-07 7.177240e-07 1
1061 0.000000e+00 7.164646e-07 2
1062 0.000000e+00 7.164646e-07 2
1063 0.000000e+00 7.164646e-07 2
1064 0.000000e+00 7.164646e-07 2
1065 7.177240e-07 7.177240e-07 1
1066 7.177240e-07 7.177240e-07 1
1067 7.177240e-07 7.177240e-07 1
1068 7.177240e-07 7.177240e-07 1
Si je lance la commande :
res=survreg(Surv(tauxsurv$tauxcompletA,tauxsurv$censure,tauxsurv$event)~var$eauproduite,dist="weibull")
j'obtiens le message d'erreur ci-dessous :
Erreur dans Surv(tauxsurv$tauxcompletA,tauxsurv$censure,tauxsurv$event) :
Stop time must be > start time
Etnt donné que je ne modélise pas une durée de vie mais un taux de casses comment puis-je adapter l'utilisation de la fonction Surv ?
Merci d'avance pour votre réponse