efficacyComputation {PoDBAY}

PoDBAY efficacy equation

Description

Function calculates the PoDBAY efficacy based on the PoD curve parameters and titer distribution parameters (mean, sd) for vaccinated and control groups.

Usage

efficacyComputation(PoDParameters, 
                    means = NA, 
                    standardDeviations = NA,
                    adjustTiters = FALSE,
                    adjustFrom = NA,
                    adjustTo = NA)

Arguments

PoDParameters

named data frame ("pmax", "slope", "et50"): PoD curve parameters

means

named list ("vaccinated", "control"): mean values of vaccinated and control subjects titers

standardDeviations

named list ("vaccinated", "control"): standard deviations of vaccinated and control subjects titers

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

Details

Efficacy = 1 - \frac{E[PoD_{vaccinated}]}{{E[PoD_{control}]} }

.

E[PoD] for each group is calculated as integral from -Inf to Inf of (titer density function) * (PoD Curve); for further details see Example2 andExpectedPoD function.

Value

efficacy: numeric value

Examples

## Data preparation
data(vaccinated)
data(control)
data(PoDParams)

## Example 1
means <- list(vaccinated = vaccinated$mean, control = control$mean)

standardDeviations <- list(vaccinated = vaccinated$stdDev, control = control$stdDev)

efficacyComputation(PoDParams, means, standardDeviations)