In 2012, Bilicki and Seikel (arXiv:1206.5130) showed that $H(z)$ data reconstructed using Gaussian Process Regression from cosmic chronometers and BAO, conclusively rules out the $R_h=ct$ model. These results were disputed by Melia and collaborators in two different works (arXiv:1304.1802 and arXiv:1802.02255), who showed using both an unbinned analysis and Gaussian Process reconstructed $H(z)$ data from chronometers, that $R_h=ct$ is favored over $\Lambda$CDM model. To resolve this imbroglio, we carry out model comparison of $\Lambda$CDM versus $R_h=ct$ by independently reproducing the above claims using latest data. We consider $H(z)$ measurements from both cosmic chronometers as well as a combination of chronometers and BAO, and perform model selection between these two models using frequentist, Bayesian, and three different information theoretical criteria. We find that with only cosmic chronometers, no one model is decisively favored. When we combine with BAO data, the information theory and Bayesian model comparison tests decisively favor the $\Lambda$CDM model.