0. (important) BiBayes.R single normal Bayesian approach
1. sigmaa/b instead of a,b
2. alpha1/2 instead of c,d
3. beta1, beta2(sp), sigma1 change to 1 dimension
4. sigma21sp and betaysp fixed at 0
5. add corresponding summary, output and plot
6. modified deltabi.f also. which may affect QRMissingBi.R (freq approach)
1. [important] first correct QRMissingBi.R for single beta1 , etc
2. beta1 -> 1 dim, sigma1 -> 1 dim, beta2sp = 0;
   sigma21 -> 1 dim; betaysp = 0, sigma21sp = 0
3. dim(param) = 2q+5 instead of 3q + 5;
   order change to (gamma1(q), gamma2(q), beta1(1), sigma1(1), betay,
   sigma2, logit(p))
4. change correspdonding ll2, residuals function
1. (important) add QRMissingBiMixMle.R for mixture of scaled normal method by MLE
2. re-define ll2Mix
1. discard tuning
2. output mcmc
1. Merge branch 'mixscalemle' into mixture-whole;
2. add QRMissingBiMixMLE.R function
3. correct single Bayes and its deltabi.f
4. correct single MLE for beta1, etc.
5. allow mix MLE feature