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)
A vector consist of mean values
A vector consist of variance values
Point estimates and std errors when combining information across MI fits (for approx method)