Recent studies have highlighted leptin, a key hormone that regulates energy intake and induces satiety, due to the worldwide prevalence of obesity. In this study, we analyzed plasma leptin measurements from 18 women with premenopausal obesity before and after bromocriptine treatment. By using underlying pulses recovered through deconvolution, we modeled the leptin secretory pulses as marked point processes and applied statistical distributions to evaluate the dynamics of leptin, including the interpulse intervals and amplitudes of the secretion. We fit the generalized inverse Gaussian and lognormal distributions to the intervals and the Gaussian, lognormal, and gamma distributions to the amplitudes of pulses. We evaluated the models' goodness of fit using statistical metrics including Akaike's information criterion, Kolmogorov-Smirnov plots, and quantile-quantile plots. Our evaluation results revealed the effectiveness of these statistical distributions in modeling leptin secretion. Although the lognormal and gamma distributions performed the best based on the metrics, we found all distributions capable of accurately modeling the timing of secretory events, leading us to a better understanding of the physiology of leptin secretion and providing a basis for leptin monitoring. In terms of pulse amplitude, the evaluation metrics indicated the gamma distribution as the most accurate statistical representation. We found no statistically significant effect of bromocriptine intake on the model parameters except for one distribution model.