PmaxEstimation {PoDBAY}

PoD curve paramater, pmax, estimation

Description

Function finds the pmax parameter of the PoD curve using control subjects summary statistics (mean, sd), observed incidence rate and previsouly estimated et50 and slope by PoDEfficacySquaredError function.

Usage

PmaxEstimation(IncidenceRate,
               params, 
               control,
               adjustTiters = FALSE, 
               adjustFrom = NA, 
               adjustTo = NA)

Arguments

IncidenceRate

numeric: observed incidence rate in overall (control) subjects

params

numeric vector: et50 and slope

control

Population-class object: control subjects (mean, sd)

adjustTiters

boolean: set to TRUE if titer values should be adjusted, for details see PoD function

adjustFrom

numeric: value specifying the detection limit, all values below the detection limit will be adjusted to adjustTo value

adjustTo

numeric: value to which titers below the detection limit will be adjusted

Value

PoD curve parameter pmax

Examples

## Example 1
data(vaccinated)
data(control)

# Assigning true efficacy 
TrueEfficacy <- 0.53

# PoD curve parameters (et50, slope) estimation
params <- PoDEfficacySquaredError(TrueEfficacy, vaccinated, control)

# Assigning incidence rate (observed incidence rate)
IncidenceRate <- 0.2

# pmax estimation
pmax <- PmaxEstimation(IncidenceRate, params, control)

# combining PoD curve parameters
PoDParams <- unlist(c(params, pmax))