## ----include = FALSE----------------------------------------------------------
knitr::opts_chunk$set(
  collapse = TRUE,
  comment = "#>"
)

## ----sim_data-----------------------------------------------------------------
library(SemiParamBernsteinDepCS)

set.seed(2026)
sim_data <- sim_bernstein_depcs(
  n = 100,
  model_type = "PH",
  copula = "gumbel",
  beta_Y = c(0.5, -0.3),
  beta_T = c(-0.2, 0.4)
)

head(sim_data)

## ----fit_model----------------------------------------------------------------
fit_ph <- fit_semiparam_bernstein_depcs(
  formula_Y = delta ~ x1 + x2,
  formula_T = time ~ x1 + x2,
  data = sim_data,
  model_type = "PH",
  copula = "gumbel",
  order_m = 2,
  n_iter = 400,
  n_burn = 100,
  seed = 42
)

print(fit_ph)

## ----summary_model------------------------------------------------------------
summary(fit_ph)

## ----predict_model------------------------------------------------------------
pred_res <- predict(fit_ph, times = seq(0.1, 2.5, length.out = 20))
head(pred_res$predictions_Y)

## ----plot_model, fig.width=6, fig.height=4------------------------------------
plot(fit_ph, type = "survival")

## ----pbc_example--------------------------------------------------------------
data(pbc_depcs)
head(pbc_depcs)