fit <- fit_single_optimizer(
formula = FEV1 ~ RACE + SEX + ARMCD * AVISIT + us(AVISIT | USUBJID),
data = fev_data,
weights = rep(1, nrow(fev_data)),
optimizer = "nlminb"
)
best_fit <- refit_multiple_optimizers(fit)
Run the code above in your browser using DataLab