reitsma function

Fit the bivariate model of Reitsma et al. (2005) and extensions.