te_model_ex <- trial_msm(
data = data_subset,
outcome_cov = c("catvarA", "nvarA"),
last_followup = 40,
model_var = "assigned_treatment",
include_followup_time = ~followup_time,
include_trial_period = ~trial_period,
use_sample_weights = FALSE,
quiet = TRUE,
glm_function = "glm"
)