QRMissing for univariate case, with single normal error distribution
remove redundent assign
give more infomation when having error solving roots
QRMissing for bivariate case, with single normal error distribution; elementwise gibbs sampling