M dwarfs, though the most abundant star in the galaxy, form only a small subset of stellar hosts with exoplanets with measured radii and masses. In this paper we analyze the Mass-Radius (M-R) relationship of planets around M dwarfs using M-R measurements for 24 exoplanets. In particular, we apply both parametric and nonparametric models and compare the two different fitting methods. We also use these methods to compare the results of the M dwarf M-R relationship with that from the Kepler sample. Using the nonparametric method, we find that the predicted masses for the smallest and largest planets around M dwarfs are smaller than that given by the relation fit to the Kepler data, but that the distribution of masses for 3 Earth Radii planets does not substantially differ between the two datasets. With future additions to the M dwarf M-R relation from TESS and instruments like the Habitable Zone Planet Finder, we will be able to characterize these differences in more detail, which will help illuminate the process of planetary formation and evolution around these stars. We release a publicly available Python code called MRExo which uses the nonparametric algorithm introduced by Ning et al. (2018) to fit the M-R relationship. Such a nonparametric fit does not assume an underlying power law fit to the measurements and hence can be used on samples spanning large mass and radii ranges. Therefore by not assuming a functional form, the fit is less biased. MRExo also offers a tool to predict mass from radius posteriors, and vice versa. This functionality can help inform observational strategies for radial velocity campaigns, such as TESS follow-up studies, as well as predict radii with microlensing planet masses.

A fundamental endeavor in exoplanetary research is to characterize the bulk compositions of planets via measurements of their masses and radii. With future sample sizes of hundreds of planets to come from TESS and PLATO, we develop a statistical method that can flexibly yet robustly characterize these compositions empirically, via the exoplanet M-R relation. Although the M-R relation has been explored in many prior works, they mostly use a power-law model, with assumptions that are not flexible enough to capture important features in current and future M-R diagrams. To address these shortcomings, a nonparametric approach is developed using a sequence of Bernstein polynomials. We demonstrate the benefit of taking the nonparametric approach by benchmarking our findings with previous work and showing that a power-law can only reasonably describe the M-R relation of the smallest planets and that the intrinsic scatter can change non-monotonically with different values of a radius. We then apply this method to a larger dataset, consisting of all the Kepler observations in the NASA Exoplanet Archive. Our nonparametric approach provides a tool to estimate the M-R relation by incorporating heteroskedastic measurement errors into the model. As more observations will be obtained in the near future, this approach can be used with the provided R code to analyze a larger dataset for a better understanding of the M-R relation.