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

## -----------------------------------------------------------------------------
library(TTM)

data("TTM_data", package = "TTM")

## -----------------------------------------------------------------------------
fit <- TTM(data_long = TTM_data,
                id = "id",
                treatment = "A",
                outcome = "Y",
                covariates = c("X1", "X2"),
                os_time = "OS_time",
                event = "event",
                t_pros = "t_pros",
                jk_block_size = NULL,
                n_inner_knot_list = range(10),
                corstr = "ar1",
                criteria = "RJC")


## -----------------------------------------------------------------------------
fit_linear <- TTM_linear(
  data_long = TTM_data,
  id = "id",
  treatment = "A",
  outcome = "Y",
  covariates = c("X1", "X2"),
  os_time = "OS_time",
  event = "event",
  t_pros = "t_pros",
  knots = 6,
  jk_block_size = NULL,
  corstr = "ar1",
  spline_type = "bs"
)

## -----------------------------------------------------------------------------
get_beta_mu(new_times = 1:10, results = fit, conf.level = 0.95)

## -----------------------------------------------------------------------------
get_beta_A(new_times = 1:10, results = fit, conf.level = 0.95)

## -----------------------------------------------------------------------------
plot_beta_mu(results = fit, conf.level = 0.95)

## -----------------------------------------------------------------------------
plot_beta_A(results = fit, conf.level = 0.95)

## -----------------------------------------------------------------------------
plot_beta_mu(results = fit_linear, conf.level = 0.95)
plot_beta_A(results = fit_linear, conf.level = 0.95)

## -----------------------------------------------------------------------------
test_beta_slope(new_times = c(5, 6, 7), fit = fit_linear, var_name = "X0")

## -----------------------------------------------------------------------------
test_beta_slope(new_times = c(5, 6, 7), fit = fit_linear, var_name = "A")