Get the EC50 Estimate from a Model

This function calculates the EC50 (the dose at which 50% of the maximum effect is observed) from a fitted model. It supports generalized linear models (GLMs) and generalized linear mixed models (GLMMs).

Usage

getEC50(mod, approximate = FALSE)

Arguments

mod

A fitted model object. This can be of class "glm" or "glmmPQL", or other models compatible with the ED function.

approximate

logical, passing into backCalcSE

Value

A data frame containing the following columns:

EC50

The estimated EC50 value.

lower

The lower bound of the 95% confidence interval for the EC50.

upper

The upper bound of the 95% confidence interval for the EC50.

se

The standard error of the EC50 estimate.

Examples

# Assuming `fit` is a fitted glm model
fit <- glm(yt ~ log(dose),family= quasibinomial(link="logit"),data=pvi_example)
ec50_result <- getEC50(fit)
print(ec50_result)
#>       EC50    lower    upper       se
#> 1 12.91593 11.03755 15.11396 1.040638