We read every piece of feedback, and take your input very seriously.
To see all available qualifiers, see our documentation.
Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.
By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.
Already on GitHub? Sign in to your account
rxode2::.rxWithWd(tempdir(), {# with temporary directory library(nlmixr2) library(rxode2) library(vpc) library(xpose) library(xpose.nlmixr2) library(ggplot2)
pheno <- function() { ini({ tcl <- log(0.008) # typical value of clearance tv <- log(0.6) # typical value of volume ## var(eta.cl) eta.cl + eta.v ~ c(1, 0.01, 1) ## cov(eta.cl, eta.v), var(eta.v) # interindividual variability on clearance and volume prop.err <- 0.1 # residual variability }) model({ cl <- exp(tcl + eta.cl) # individual value of clearance v <- exp(tv + eta.v) # individual value of volume ke <- cl / v # elimination rate constant d/dt(A1) = - ke * A1 # model differential equation cp = A1 / v # concentration in plasma cp ~ prop(prop.err) # define error model }) }
fit_focei <- nlmixr(pheno, pheno_sd, "focei", control=list(print=0), table=list(cwres=TRUE, npde=TRUE))
auto3 <- covarSearchAuto(fit_focei, varsVec = c("v", "cl"), covarsVec = c("WT","APGR"), restart = TRUE, searchType = "backward", pVal=list(fwd=0.01,bck=0.001)) })
.... # it was running for quite a while and everything seems OK, then it comes the error below;
Error in backwardSearch(varsVec, covarsVec, catvarsVec, fitorig = fit, : object 'covInfo' not found
The text was updated successfully, but these errors were encountered:
Thanks for reporting @snailittle
I will look at it when I get a chance.
Sorry, something went wrong.
No branches or pull requests
rxode2::.rxWithWd(tempdir(), {# with temporary directory
library(nlmixr2)
library(rxode2)
library(vpc)
library(xpose)
library(xpose.nlmixr2)
library(ggplot2)
pheno <- function() {
ini({
tcl <- log(0.008) # typical value of clearance
tv <- log(0.6) # typical value of volume
## var(eta.cl)
eta.cl + eta.v ~ c(1,
0.01, 1) ## cov(eta.cl, eta.v), var(eta.v)
# interindividual variability on clearance and volume
prop.err <- 0.1 # residual variability
})
model({
cl <- exp(tcl + eta.cl) # individual value of clearance
v <- exp(tv + eta.v) # individual value of volume
ke <- cl / v # elimination rate constant
d/dt(A1) = - ke * A1 # model differential equation
cp = A1 / v # concentration in plasma
cp ~ prop(prop.err) # define error model
})
}
fit_focei <- nlmixr(pheno, pheno_sd, "focei",
control=list(print=0),
table=list(cwres=TRUE, npde=TRUE))
auto3 <- covarSearchAuto(fit_focei, varsVec = c("v", "cl"), covarsVec = c("WT","APGR"),
restart = TRUE,
searchType = "backward", pVal=list(fwd=0.01,bck=0.001))
})
.... # it was running for quite a while and everything seems OK, then it comes the error below;
Error in backwardSearch(varsVec, covarsVec, catvarsVec, fitorig = fit, :
object 'covInfo' not found
The text was updated successfully, but these errors were encountered: