powered by
prepare_data(df, df_reference, dependent, relevel = TRUE, drop.nas = c("dependent", "predictors", "all", "none"))
drop_na