if (FALSE) {
input_data <- pkData
model <-
pkmodel(numCompartments = 2,
data = input_data,
ID = "Subject",
Time = "Act_Time",
A1 = "Amount",
CObs = "Conc",
workingDir = tempdir())
model <-
addCovariate(model,
covariate = "BodyWeight",
direction = "Backward",
center = "Mean",
effect = c("V", "Cl"))
# multicore
multicoreHost <-
hostParams(parallelMethod = "Multicore",
hostName = "multicore",
numCores = 4,
sharedDirectory = tempdir())
# specify scenarios
CovariateEffectNames <- listCovariateEffectNames(model)
combinations <-
combn(c("", CovariateEffectNames),
length(CovariateEffectNames),
simplify = FALSE)
scenarioNames <-
lapply(combinations,
function(x) {paste(x, collapse = " ")})
scenarios <-
lapply(scenarioNames,
function(x, CovariateEffectNames) {
CovariateCombinations <- unlist(strsplit(x, " ", fixed = TRUE))
scenarioIndex <-
paste(which(CovariateEffectNames %in% CovariateCombinations,
arr.ind = TRUE),
collapse = ", ")
NlmeScenario(trimws(x), scenarioIndex)
},
CovariateEffectNames)
res <-
sortfit(model,
hostPlatform = multicoreHost,
params = engineParams(model, numIterations = 5, fastOptimization = TRUE),
sortColumns = SortColumns("Gender"),
scenarios = scenarios)
}
Run the code above in your browser using DataLab