ph_fail_example.r

# Load the survival package
library(survival)

# Read in the CSV
data <- read.csv('C:/users/dbatten/google drive/laptop files/personal/time dependent coefficients/tdc_data.csv')

# Fit a cox model of a survival object on the user's gender
fit <- coxph(Surv(time, churn) ~ male, data = data)

# Have a gander at the results
summary(fit)

# Make sure we haven't violated the proportional hazards assumption
cox.zph(fit)