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

## ----smoke-data, include=FALSE------------------------------------------------
set.seed(2026)
n_smoke <- 140
p_smoke <- 6
X_smoke <- matrix(rnorm(n_smoke * p_smoke), n_smoke, p_smoke)
colnames(X_smoke) <- paste0("x", seq_len(p_smoke))
beta_smoke <- c(1.0, -0.8, 0.6, 0, 0, 0)
T_smoke <- rexp(n_smoke, rate = 0.05 * exp(drop(X_smoke %*% beta_smoke)))
C_smoke <- rexp(n_smoke, rate = 0.025)
time_smoke <- pmin(T_smoke, C_smoke)
status_smoke <- as.integer(T_smoke <= C_smoke)
library(selSurv)
fit_smoke <- coxlasso_fit(X_smoke, time_smoke, status_smoke, lambda_choice = "1se")

## ----fixed, eval=FALSE--------------------------------------------------------
# fixed <- coxlasso_boot(
#   fit,
#   inference = "conditional",
#   B = 500,
#   conf_level = 0.95,
#   seed = 2026
# )
# summary(fixed)

## ----fixed-smoke, include=FALSE-----------------------------------------------
fixed_smoke <- coxlasso_boot(
  fit_smoke,
  inference = "conditional",
  B = 3L,
  conf_level = 0.95,
  seed = 2026
)
summary(fixed_smoke)