‎Maximum likelihood estimation of multivariate distributions needs solving a optimization problem with large dimentions (to the number of unknown parameters) but two‎- ‎stage estimation divides this problem to several simple optimizations‎. ‎It saves significant amount of computational time‎. ‎Two methods are investigated for estimation consistency check‎. ‎We revisit Sankaran and Nair's bivariate Pareto distribution as an example‎. ‎Two data sets (simulated data and real data) have been analyzed for illustrative purposes‎.