test_that("Survival-binary model works for Colorectal data",
          {
            skip_on_cran()
            data("colorectal")

            fitted_model = survbin(data = colorectal, true = surv, trueind = SURVIND,
                                   surrog = responder, trt = TREAT, center = CENTER,
                                   trial = TRIAL, patientid = patientid)
            rsquared_theta = c(round(fitted_model$Trial.R2[1,1], digits = 4), round(fitted_model$Indiv.GlobalOdds[1,1], digits = 4))
            expect_equal(rsquared_theta, c(0.4419,  4.9109), tolerance=1e-3)
          })