Comparison of Maximum Likelihood Estimation and Bayesian with Generalized Gibbs Sampling for Ordinal Regression Analysis of Ovarian Hyperstimulation Syndrome

Background and Objectives: Analysis of ordinal data outcomes could lead to bias estimates and large variance in sparse one. The objective of this study is to compare parameter estimates of an ordinal regression model under maximum likelihood and Bayesian framework with generalized Gibbs sampling. The models were used to analyze ovarian hyperstimulation syndrome data.   Methods: This study used the data from 138 patients of a clinical trial phase III to compare the efficacy of intravenous Albumin and Cabergoline in prevention of ovarian hyperstimulation syndrome. The original study was done between 2010 to 2011 in Royan institute. We compared maximum likelihood and Bayesian estimation with generalized Gibbs sampling for an ordinal regression model based on confidence intervals and standard errors. The model were fit through R 3.3.2 software version.   Results: Markov Chain Monte Carlo results reduced the standard errors for estimates and consequently, narrower confidence intervals. Autocorrelations for generalized Gibbs sampler reached to zero in compare to standard Gibbs sampler for shorter time.   Conclusion: It seems that confidence intervals of an ordinal regression model are shorter for generalized Gibbs sampler in compare to standard Gibbs and maximum likelihood. It suggests doing more studies to warrant the results.