RMST.r

# Load the survival library
 
library(survival)
 
# Read in the data and create a "survival object"
 
netlixx_rmst <- read.csv('C:/Users/dayne/Desktop/NetLixxRMST.csv')
netlixx_rmst$survival <- Surv(netlixx_rmst$time, netlixx_rmst$churned == 1)

# Fit a survival curve for each group

fitted <- survfit(survival ~ credit, data = netlixx_rmst)

# Get the 180-day RMST for each group

print(fitted, print.rmean=getOption("survfit.print.rmean"), rmean = 180)