data <- generateSurvData(nP = 5,
nT = 5,
lambdaP = 0.03943723,
lambdaT = 0.5*0.03943723,
endTime = 40,
seed = 2006)
survObj <- survival::Surv(data$time, data$status)
survDataFrame <- as.data.frame(as.matrix(survObj))
y0Index <- which(data$group=="P")
y1Index <- which(data$group=="T")
timeNow <- 4
timeBefore <- 0
computeStatsForLogrank(survDataFrame, y0Index, y1Index, timeNow, timeBefore)
timeNow <- 13
timeBefore <- 4
computeStatsForLogrank(survDataFrame, y0Index, y1Index, timeNow, timeBefore)
Run the code above in your browser using DataLab