# test_that("Survival- Survival model works for ovarian data", {
#   skip_on_cran()
#   data("Ovarian")
#
#   Ovarian2 <- Ovarian %>%
#     dplyr::group_by(Center) %>%
#     dplyr::filter(sum(Treat == 1) >= 3 &
#                     sum(Treat == 0) >= 3) %>%
#     dplyr::ungroup()
#
#   fitted_model =  MetaAnalyticSurvSurv(
#     data = Ovarian2,
#     true = Surv,
#     trueind = SurvInd,
#     surrog = Pfs,
#     surrogind = PfsInd,
#     trt = Treat,
#     center = Center,
#     trial = Center,
#     patientid = Patient,
#     copula = "Plackett",
#     adjustment = "weighted"
#   )
#   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.8721, 0.9671), tolerance =
#                  1e-2)
# })