# NOT RUN {
sfit = survfitDistributed(vdra_fit_cox_A)
print(sfit)
plot(sfit)
# From Data Partner 1
sfit = survfitDistributed(vdra_fit_cox_A, ~Exposure, data = vdra_data[, c(3:4, 5:7)])
print(sfit)
plot(sfit)
# From Data Partner 2
sfit = survfitDistributed(vdra_fit_cox_B, ~Race + Sex, data = vdra_data[, 8:11])
print(sfit)
plot(sfit, merge = TRUE)
# }
Run the code above in your browser using DataLab