Bivariate Pareto Distribution

Implements the EM algorithm with one-step Gradient Descent method to estimate the parameters of the Block-Basu bivariate Pareto distribution with location and scale. We also found parametric bootstrap and asymptotic confidence intervals based on the observed Fisher information of scale and shape parameters, and exact confidence intervals for location parameters. Details are in Biplab Paul and Arabin Kumar Dey (2023) "An EM algorithm for absolutely continuous Marshall-Olkin bivariate Pareto distribution with location and scale"; E L Lehmann and George Casella (1998) "Theory of Point Estimation"; Bradley Efron and R J Tibshirani (1994) "An Introduction to the Bootstrap"; A P Dempster, N M Laird and D B Rubin (1977) "Maximum Likelihood from Incomplete Data via the EM Algorithm".