Performs comprehensive comparison of dose-response models using multiple criteria including model fit, certainty of protection, and steepness.
Usage
drcCompare(
modRes = NULL,
modList = NULL,
trend = "Decrease",
CI = c("delta", "inv", "bmd-inv"),
...
)Arguments
- modRes
Optional. Result object containing model comparison information. If provided, modList parameter is ignored.
- modList
Optional. List of fitted dose-response models to compare. Required if modRes is NULL.
- trend
Character string specifying the expected trend direction. Must be "Decrease" or "Increase".
- CI
Character string specifying the confidence interval method. One of "delta", "inv", or "bmd-inv". Defaults to "delta".
- ...
Additional arguments passed to internal functions
Value
A data frame containing model comparison results with columns:
Original comparison metrics
Certainty_Protection: Measure of protection certainty
Steepness: Model steepness evaluation
No Effect p-val: P-value for test against no-effect model
Examples
if (FALSE) { # \dontrun{
data("dat_medium")
dat_medium <- dat_medium %>% mutate(Treatment=factor(Dose,levels=unique(Dose)))
dat_medium$Response[dat_medium$Response < 0] <- 0
mod <- drm(Response~Dose,data=dat_medium,fct=LL.3())
fctList <- list(LN.4(),LL.4(),W1.3(),LL2.2())
res <- mselect.plus(mod,fctList = fctList )
modList <- res$modList
# Compare models using existing comparison results
drcCompare(modRes = res, trend = "Decrease")
} # }