parameters for analyzing bivariate skew normal data with non monotone missing values is
presented. A simulation study is implemented to investigate the performance of the
presented algorithm. Results of an application are also reported where a Bootstrap
approach is used to find the variances of the parameter estimates.