The aim of this work was to study the reaction kinetics of neopentyl glycol esterification with different carboxylic acids (propionic acid, isobutyric acid and 2-ethylhexanoic acid). The experiments were carried out in a batch reactor as reactive distillation removing continuously one reaction product (water) by purging nitrogen through the reactor. Commercial heterogeneous gellular (Dowex 50WX2) and macroreticular (Amberlyst 15) type resins and homogeneous para-toluene sulfonic acid were compared as catalysts. The gellular type resin was dried before loading into the reactor. When performing esterification of neopentyl glycol with different acids the observed formation of monoester was clearly much faster than the consecutive esterification of monoester to diester. A disproportionation reaction was also found to take place. The reliability of the estimated Arrhenius parameters was analyzed using Bayesian inference and Markov chain Monte Carlo methods (MCMC). MCMC analysis was also used to study uncertainties in the model predictions and in some additional parameters, such as reaction rate constants calculated from the estimated parameters. This cannot be done by classical regression analysis.