postM function

Calculate posterior distribution of M and extract statistics (M* and CI)