frwd_selection: Forward covariate selection for nlme-base non-linear mixed...

Description Usage Arguments Value Examples

View source: R/vpc.R

Description

Implements forward covariate selection for nlme-based non-linear mixed effect models

Usage

1
frwd_selection(base, cv, dat, cutoff = 0.05)

Arguments

base

base model

cv

a list of candidate covariate to model parameters

dat

model data

cutoff

significance level

Value

an nlme object of the final model

Examples

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
12
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))

nlmixr documentation built on June 21, 2018, 5:04 p.m.