Survival analysis and log-rank comparison in R

survival.r

# Survival analysis in R

# Load library (install first if needed)
library(survival)
library(survMisc)

# create a sample data set
data <- lung

# fit a survival model
fit <- survfit(Surv(time, status) ~ sex, data = data)

# print the summary of the model
summary(fit)

# plot the survival curve
autoplot(fit, 
            timeTicks = "custom", times = seq(1, 1000, 100), 
            plotTable = TRUE,
            CI = TRUE,
            bands = TRUE)

# compare survival curves
result <- survdiff(Surv(time, status) ~ sex, data = lung)

# print the result of the log-rank test
print(result)