# NOT RUN {
data <- adnimerge %>% dplyr::filter(VISCODE == 'bl')
# fit model with one subgroup (DX_bl) and also the entire data
model <- data %>% aba_model() %>%
set_groups(
everyone(),
DX_bl %in% c('MCI', 'AD')
) %>%
set_outcomes(ConvertedToAlzheimers, CSF_ABETA_STATUS_bl) %>%
set_predictors(PLASMA_ABETA_bl, PLASMA_PTAU181_bl, PLASMA_NFL_bl) %>%
set_stats(
stat_glm(std.beta = TRUE)
) %>%
fit()
model_summary <- model %>% aba_summary()
# }
Run the code above in your browser using DataLab