# NOT RUN {
# use built-in data
data <- adnimerge %>% dplyr::filter(VISCODE == 'bl')
# fit an aba model
model <- data %>% aba_model() %>%
set_groups(everyone()) %>%
set_outcomes(PET_ABETA_STATUS_bl) %>%
set_predictors(
PLASMA_PTAU181_bl,
PLASMA_NFL_bl,
c(PLASMA_PTAU181_bl, PLASMA_NFL_bl)
) %>%
set_covariates(AGE, GENDER, EDUCATION) %>%
set_stats('glm') %>%
fit()
# default aba summary
model_summary <- model %>% aba_summary()
# create an aba control object to customize the summary
my_control <- aba_control(include_covariates = FALSE)
# summarise model with th custom aba control - notice covariates
# wont be included in the tables when you print the summary to console
model_summary2 <- model %>% aba_summary(control = my_control)
# }
Run the code above in your browser using DataLab