Function to implement Rubins 1987 method to estimate point estimates and std errors when combining information across MI fits (for approx method)
Rubin.MI(mean.vec, variance.vec)Point estimates and std errors when combining information across MI fits (for approx method)