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

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