R/plots.R

In barryrowlingson/cransurv: Survival analysis of CRAN packages

library(survival)
library(ggfortify)   # autoplot for survfit

sf <- survfit(Surv(days, event) ~ nc, data=sd2)
autoplot(sf, main = "Survival as a Function of 'Needs Compilation'")

sd2$meanBig <- sd2$meansize > mean(sd2$meansize)
sf <- survfit(Surv(days, event) ~ meanBig, data=sd2)
autoplot(sf, main = "Survival as a Function of 'Larger than Median Size'")

sd2$manyDep <- sd2$Depends + sd2$Imports >= 6
sf <- survfit(Surv(days, event) ~ manyDep, data=sd2)
autoplot(sf, main = "Survival as a Function of 'Many Depends (>= 6)'")

sd2$manySug <- sd2$Suggests >= 4
sf <- survfit(Surv(days, event) ~ manySug, data=sd2)
autoplot(sf, main = "Survival as a Function of 'Many Suggests (>= 4)'")

sd2$years <- sd2$days / 365.25
sd2$manySug <- sd2$Suggests >= 4
sf <- survfit(Surv(years, event) ~ manySug, data=sd2)
autoplot(sf, main = "Survival as a Function of 'Many Suggests (>= 4)'", 
         xlab="Years", ylab="Predicted Probability")