PmaxEstimation: PoD curve paramater, pmax, estimation

In PoDBAY: Vaccine Efficacy Estimation Package

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))

PoDBAY documentation built on Sept. 21, 2021, 5:08 p.m.