MarionGuernote / Projet_annuel

0 stars 0 forks source link

script_test1 #1

Open raph06 opened 4 years ago

raph06 commented 4 years ago

Bonjour à tous,

Je vais mettre en ligne un premier script de test. Dans le script vous trouverez des infos, du code à extraire et des schémas de ce qui est attendu.

Vous verrez beaucoup de help(function) (ex: help(rep), help(data) etc), cela vous permet d'accéder à la doc. les infos que vous devez retenir: input, output et arguments ... (@MarionGuernote ! )

@TerryGOUIN @ThibaultMeekel @MaxLrs cela vous aidera à vous familiariser avec le langage.

Le script est simple mais il arbore déjà des difficultés pour lesquelles je propose des solutions. (attach et dico pour répertorier les variables que vous trouvez détaillé plus bas et dans le script)

Lisez bien les commentaires, et posez moi des questions !

milestone au 1 Déc pour lecture et doc.

raph06 commented 4 years ago

Hi,

regarding:

survival model

formula=Surv(lung$time,lung$status)

This computes the necessary information for a survival analysis. simple assignation will do, but 2 parameters for this function.

#   -----                      
# | formula |   <--surv-- ¦time¦  
#      ^-----------surv-- ¦status¦   

often a covariate is used: formula (survival info) ~ (function of) old (age of the patients) this is not only specific to survival analysis but to most regression models The following computes a survival plot in addition to extra information into the variable fit.

fit <- survfit(formula ~ old, data = data.frame(lung)) this computes the stats. fit_table <- coxph(formula ~ old) so

#   -----                      
# | fit |   <--survfit ~ old-- ¦formula¦  
#   -----                      

can do it

or

#   -----                      
# | fit |   <-- ~old-- ¦formula¦  
#   -----