Skip to content

Instantly share code, notes, and snippets.

@drewconway

drewconway/chat_roulette_surv3.R

Created April 2, 2010 15:28
Show Gist options
  • Save drewconway/353258 to your computer and use it in GitHub Desktop.
Save drewconway/353258 to your computer and use it in GitHub Desktop.
Download ZIP
Raw
chat_roulette_surv3.R
### Perform analysis and make pretty pictures ###
png("surv_plot.png",height=1000,width=1000,res=100)
par(mfrow=c(2,2))
survplot(survfit(Surv(time,censor)~type,data=subset(surv_data,type=="lonely")),what="survival",conf="bands",xlab="Minutes")
title("Survival function for seeing a lonely guy")
survplot(survfit(Surv(time,censor)~type,data=subset(surv_data,type=="penis")),what="survival",conf="bands",xlab="Minutes",main="Survival function for seeing a penis")
title("Survival function for seeing a penis")
survplot(survfit(Surv(time,censor)~type,data=subset(surv_data,type=="drunk")),what="survival",conf="bands",xlab="Minutes",main="Survival function for seeing two or more drunk people")
title("Survival function for seeing two or more drunk people")
survplot(survfit(Surv(time,censor)~type,data=subset(surv_data,type=="chick")),what="survival",conf="bands",xlab="Minutes",main="Survival function for seeing a woman")
title("Survival function for seeing a woman")
dev.off()
# Summary stats and output data
summary(survfit(Surv(time,censor)~type,data=surv_data))
write.csv(surv_data,"cr_surv.csv")
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment