The Megamaser Cosmology Project inferred a value for the Hubble constant given by $H_{0}=73.9 \pm 3.0 $ km/sec/Mpc. This value was obtained using Bayesian inference by marginalizing over six nuisance parameters, corresponding to the velocities of the megamaser galaxy systems. We obtain an independent estimate of the Hubble constant with the same data using frequentist inference. For this purpose, we use profile likelihood to dispense with the aforementioned nuisance parameters. The frequentist estimate of the Hubble constant is given by $H_{0}=73.5^{+3.0}_{-2.9}$ km/sec/Mpc and agrees with the Bayesian estimate to within $0.2\sigma $ , and both approaches also produce consistent confidence/credible intervals. Therefore, this analysis provides a proof-of-principle application of profile likelihood in dealing with nuisance parameters in cosmology, which is complementary to Bayesian analysis.
