bibliografia
a.a. 2011-2012
Università degli Studi di Trieste
Facoltà di Medicina e Chirurgia a.a. 2011-2012
Scuole di Specialità di Medicina e Chirurgia
Corso di Statistica Medica
docenti:Lucio Torelli, Massimo Borelli, Giulia Barbati, Fabiola Giudici
analisi di sopravvivenza
library(survival)
www <- "http://www.dmi.units.it/~borelli/dataset/breastcancer.txt"
breastcancer <- read.table(www, header = TRUE)
attach(breastcancer)
str(breastcancer)
head(breastcancer)
summary(breastcancer)
kaplan_meier <- survfit(Surv(los, discharge) ~ treatment)
summary(kaplan_meier)
plot(kaplan_meier, lty=c(5,3), col = c("violet", "orange"), xlab="los", lwd = 2)
legend( 15, 0.7, levels(treatment), lty=c(5,3), col = c("violet", "orange") , lwd = 2)
logrank <- survdiff(Surv(los, discharge) ~ treatment)
logrank
modello1 <- survreg(Surv(los, discharge) ~ age + treatment )
summary(modello1)
modello2 <- survreg(Surv(los, discharge) ~ age + drainage + treatment )
summary(modello2)
modello3 <- survreg(Surv(los, discharge) ~ drainage * treatment )
summary(modello3)
modello4 <- survreg(Surv(los, discharge) ~ drainage + treatment )
summary(modello4)
detach(breastcancer)
rm(breastcancer)