hbim (version 1.1.2)

hbim-internal: Internal functions not be to called directly

Description

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.

Usage

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)

Value

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

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

Author

M.P. Fay

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