bayesint

This repository includes the Python code for calculating Bayesian credible interval of ratios of independent beta distributions. It has been used in Methods for calculating credible intervals for ratios of beta distributions with application to relative risks of premature death during the second plague pandemic

Getting Started

Prerequisites

• Python 2.7+

Installation

Clone the repository then run

python setup.py install

Usage

To get the relative risk of a contingency table given by

run

from bayesint import rel_risk
rel_risk(56, 126, 366, 354)
# 236/549

and to obtain the equal-tailed quantile interval for the data given in the contingency table run

from bayesint import eqt_int_frac
eqt_int_frac(56, 126, 366, 354, (0, 0, 0, 0), "risk", 0.05, "estim")
# (236/549, 0.184135819539239, 0.667343920284484)

Authors

Maria Bekker-Nielsen Dunbar and Tom Finnie

Contributing

This project is intended to be a safe, welcoming space for collaboration, and contributors are expected to adhere to the Contributor Covenant code of conduct.

1. Fork it
2. Create your feature branch (git checkout -b my-new-feature)
3. Commit your changes (git commit -am 'Add some feature')
4. Push to the branch (git push origin my-new-feature)
5. Create a new Pull Request

License

Public Health England 2017

This project is licensed under the Open Government License, see LICENSE for details