PoDParamPointEstimation {PoDBAY}

PoD curve point estimate

Description

Function returns PoD curve parameters corresponding to the point estimate of PoD curve.

Usage

PoDParamPointEstimation(resultsPriorReset, 
                        titers = seq(from = 0, to = 20, by = 0.01), 
                        optim_titers = FALSE)

Arguments

resultsPriorReset

named data frame ("pmax", "slope", "et50"): set of estimated PoD curve parameters before resetting the disease status; for further details see PoDParamEstimation function.

titers

numeric vector: a grid of titers for PoD curve point estimate calculation

optim_titers

logical: TRUE for a predefined sequence of titers

Details

For each of estimated PoD curves in resultsPriorReset, the function values (probabilities of disease, PoD) for provided grid of titers are calculated.

Median of function values (PoDs) at each provided titer is calculated.

Subsequently, the PoD curve model is fitted to the median datapoins using fitPoD function, in order to get PoD curve parameters close to this median curve.

Value

paramsPointEstimate: named data frame of PoD curve parameters corresponding to the PoD curve point estimate

Examples

## Data preparation
data(estimatedParameters)

## Example 1
# titers for which we want to optimize the functional values
titers <- seq(from = 0, to = 20, by = 0.01)

# Point estimate of PoD curve
PoDParamPointEstimation(estimatedParameters$resultsPriorReset, titers)