In this work, we present the implementation of the New Exponentiated Exponential (NEE) distribution within the gamlss package in R, providing users with a flexible and robust tool for statistical modeling. This implementation enables the computation of probabilities, quantiles, and random numbers, as well as parameter estimation and the development of regression models within the GAMLSS framework. To evaluate its performance, we conducted a simulation study to evaluate the estimation process using the mean and mean squared error. The results demonstrate the effectiveness of the NEE distribution in capturing data variability and improving model accuracy. This contribution extends the capabilities of the gamlss package, offering R users an advanced resource for data analysis and regression modeling, especially in contexts where flexibility and precision are crucial.