Je suis nouveau dans l'analyse de survie et j'ai récemment appris qu'il existe différentes façons de le faire en fonction d'un certain objectif. Je suis intéressé par la mise en œuvre réelle et la pertinence de ces méthodes.
On m'a présenté les traditionnels risques proportionnels de Cox , les modèles de temps de défaillance accélérés et les réseaux de neurones (perceptron multicouche) comme méthodes pour obtenir la survie d'un patient compte tenu de son temps, de son statut et d'autres données médicales. L'étude devrait être déterminée sur cinq ans et l'objectif est de donner chaque année des risques de survie pour de nouveaux enregistrements.
J'ai trouvé deux cas où d'autres méthodes ont été choisies par rapport au PH de Cox:
J'ai trouvé " Comment obtenir des prédictions en termes de temps de survie à partir d'un modèle Cox PH " et il a été mentionné que:
Si vous êtes particulièrement intéressé à obtenir des estimations de la probabilité de survie à des moments particuliers, je vous orienterais vers des modèles de survie paramétriques (également appelés modèles de temps de défaillance accélérés) . Ceux-ci sont implémentés dans le package de survie pour R, et vous donneront des distributions de temps de survie paramétriques, dans lesquelles vous pouvez simplement brancher l'heure qui vous intéresse et récupérer une probabilité de survie.
Je suis allé sur le site recommandé et en ai trouvé un dans le
survival
package - la fonctionsurvreg
.Des réseaux de neurones ont été suggérés dans ce commentaire :
... Un avantage des approches neuronales de l'analyse de survie est qu'elles ne reposent pas sur les hypothèses qui sous-tendent l'analyse de Cox ...
Une autre personne avec la question « Modèle de réseau neuronal R avec vecteur cible en sortie contenant des prédictions de survie » a donné un moyen exhaustif de déterminer la survie à la fois dans les réseaux neuronaux et Cox PH.
Le code R pour obtenir la survie serait comme ceci:
mymodel <- neuralnet(T1+T2+T3+T4+T5~covar1+covar2+covar3+..., data=mydata, hidden=1) compute(mymodel,data=mydata)
Je suis allé sur les forums R et j'ai trouvé cette réponse dans la question " prédire.coxph et prédire.survreg " :
En effet, à partir de la
predict()
fonction du,coxph
vous ne pouvez pas obtenir directement des prédictions «temporelles», mais uniquement des scores de risque linéaires et exponentiels. En effet, pour obtenir l'heure, un aléa de base doit être calculé et il n'est pas simple car il est implicite dans le modèle de Cox.
Je me demandais si lequel des trois (ou deux considérant les arguments sur Cox PH) est le meilleur pour obtenir des pourcentages de survie pour les périodes d'intérêt? Je ne sais pas lequel d'entre eux utiliser dans l'analyse de survie.