ED.DRMod: Calculate EDp estimator for a DRMod object

Description

This function calculates the effective dose of a specified level for a DRMod object. Here the effective dose of level p is the smallest dose that achieves p percent of the maximum effect observed within the dose range (and not with respect to the asymptotic maximum effect).

Usage

## S3 method for class 'DRMod':
ED(object, p, doseSeq = NULL, lenDose = 101,
   direction = c("increasing", "decreasing"),
   data = getData(object), ...)

Arguments

object

An object of class DRMod

Numeric specifying the level used for the EDp estimate.

doseSeq

Sequence of values allowed for the dose estimator. Should include 0 (ie placebo) as a first value.

lenDose

If doseSeq is not specified an equally spaced grid of lenDose values between placebo and the maximal dose in the study is used.

direction

Character specifying the direction of the dose effect (ie whether the response should increase or decrease with increasing dose).

data

Data frame containing the variables on which the DRMod object was fitted

...

Additional arguments

Value

ED dose estimate

References

Bornkamp B., Pinheiro J. C., and Bretz, F. (2009). MCPMod: An R Package for the Design and Analysis of Dose-Finding Studies, Journal of Statistical Software, 29(7), 1--23

Examples

Run this code

data(biom)
fit <- fitDRModel(resp ~ dose, biom, "logistic")
ED(fit, p = c(0.25,0.5,0.95))