# NOT RUN {
dat <- theo_md
dat$LOGWT <- log(dat$WT)
dat$TG <- (dat$ID < 6) + 0 # dummy covariate
specs <- list(
fixed = list(lKA = lKA ~ 1, lCL = lCL ~ 1, lV = lV ~ 1),
random = pdDiag(lKA + lCL ~ 1),
start = c(0.5, -3.2, -1)
)
fit0 <- nlme_lin_cmpt(dat, par_model = specs, ncmt = 1)
cv <- list(lCL = c("WT", "TG"), lV = c("WT"))
fit <- frwd_selection(fit0, cv, dat)
print(summary(fit))
# }
Run the code above in your browser using DataLab