1. Grâce à ce formulaire vous pouvez proposer une astuce ou un script sur R.
  2. Votre script doit pouvoir être lancé en l'état. Veuillez penser à :
    • inclure le chargement des "library" nécessaires
    • construire un petit jeu de données si besoin est.
    • commenter les lignes de codes pour en faciliter la compréhension.
  3. Le titre de votre script ou astuce doit être clair et explicite.
  4. Pensez à mettre votre code entre les balises [R] et [/R]; Pour cela, vous pouvez utiliser le bouton

Vous pouvez utiliser vos comptes Facebook, twitter ou google pour vous identifer (google est compatible yahoo, openID...)
L'ideal étant de vous connecter si vous avez un compte utilisateur, ou faire une demande de compte utilisateur si vous n'en avez pas encore.
Créer un compte va vous permettre de pouvoir éditer vos codes et de mettre en avant votre site internet.
Sinon vous pouvez soumettre anonymement en remplissant les champs ci-après.







Choisissez les catégories correspondantes à votre Code:

  • algorithmique
  • Analyse de survie
  • base indispensable
  • bayésien
  • configuration de R
  • exportation de données
  • fonctions utiles
  • graphique
  • importation de données
  • inférence
  • manipulation de données
  • message d'erreur
  • modélisation
  • Non classé
  • optimisation
  • planification
  • programmer avec R
  • regression linéaire
  • Test
  • tidyverse
  • Transformation de données


Analyse de survie

Comment obtenir l’estimateur de Kaplan-Meier pour une analyse de survie avec R ?

4
dans Analyse de survie
- ça ne sert à rien -- c\'est interessant - (score de +7 sur 11 votes)
Loading ... Loading ...
L’estimateur de Kaplan-Meier donne la fonction de survie non paramétrique.
Pour l’obtenir sous R on peut utiliser le package survival.

On se place ici dans un cas très simple où il n’y a ni censure ni troncature.
Pour bien comprendre le code, je vous conseille vivement de regarder la documentation du package en question!!


 
 
 #survival analysis
 
ls()
rm(list=ls())
 
library(survival)
 
#on crée un jeu de données correspondant à des durées (étudiées dans l'analyse de survie)
 
 z<-c(14,14,14,17,17,17,17,17,17,17,17,17,17,17,17,17,17,17,17,17,17,17,20,20,20,20,20,20,20,20,20,20,20,20,20,20,20,20,20,20,20,20,20,20,20,20,20,20,23)
d<-data.frame(delay=z)
 
#on crée une colonne status, ici tous les individus sont "morts" pendant l'expérience 
# mort au sens de l'analyse de survie
 
 d$status<-1 
 
 s<-survfit(Surv(d$delay,d$status)~1)
plot(s,main="survival function")
 
 
 

Proposé par Melen.

Déjà 4 commentaires. Cliquez ici pour réagir.
Formation logiciel R