hbim-internal function

Internal functions not be to called directly

These are functions that are called by the hbrr and hbpp

functions. They integrate 1, 2 or 3 component models or integrate by simulation. Special functions ending in .rhoeq1 are for when rho equals 1.

hbpp.integrate1(MU, V, A = 1, RP = 0.1, ...)
hbpp.integrate2(MU, V, A = c(1, 1), RP = 0.1, ...)
hbpp.integrate3(MU, V, A = c(1, 1, 1),RP=0.1, ...)
hbpp.simulate(MU, V, A, RP = 0.1, NSIM = 10^5)
hbrr.integrate1(MU, V, A = 1, ...)
hbrr.integrate2(MU, V, A = c(1, 1), ...)
hbrr.integrate2.rhoeq1(MU, V, A = c(1, 1), ...)
hbrr.integrate3(MU, V, A = c(1, 1, 1), ...)
hbrr.integrate3.rhoeq1(MU, V, A = c(1, 1, 1), ...)
hbrr.simulate(MU, V, A, NSIM = 10^4)

Arguments

• MU: mean vector of the log10 dose
• V: variance matrix of the log10 dose
• A: vector of slope parameters in the Hill model, one for each component
• NSIM: number of simulations
• RP: protection bound, an individual is protected if relative risk is greater than RP
• ...: additional parameters to pass to the integrate function

Returns

a numeric value of the expected relative risk or percent protected.

References

Saul A, Fay MP (2007). Human Immunity and the Design of Multi-Component, Single Target Vaccines. PLoS ONE 2(9): e850. doi:10.1371/jounal.pone.0000850

Author(s)

M.P. Fay