By adding over-relaxation factor to PXEM (Parameter Expanded Expectation Maximization) method, the MOPXEM (Monotonically Overrelaxed Parameter Expanded Expectation Maximization) method is obtained. Compare it with the existing EM (Expectation-Maximization)-like methods. Then, distribute and process five methods and compare them, achieving good performance in convergence speed and result quality.The philosophy of the package is described in Guo G. (2022) <doi:10.1007/s00180-022-01270-z>.

Estimate parameters in Dirichlet-Multinomial and compute log-likelihoods.

Nonparametric kernel density estimation, bandwidth selection,
    and other utilities for analyzing directional data. Implements the estimator
    in Bai, Rao and Zhao (1987) <doi:10.1016/0047-259X(88)90113-3>, the
    cross-validation bandwidth selectors in Hall, Watson and Cabrera (1987)
    <doi:10.1093/biomet/74.4.751> and the plug-in bandwidth selectors in
    García-Portugués (2013) <doi:10.1214/13-ejs821>.

Contains the support functions for the Time Series Analysis book.
	     We present a function to calculate MSE and MAE for inputs
	     of actual and forecast values.  We also have the code for
	     disaggregation as found in Wei and Stram
	     (1990, <doi:10.1111/j.2517-6161.1990.tb01799.x>),
	     and Hodgess and Wei (1996, "Temporal Disaggregation of Time
	     Series").

The twoStepsBenchmark() and threeRuleSmooth() functions allow you to 
    disaggregate a low-frequency time series with higher frequency time series, 
    using the French National Accounts methodology. The aggregated sum of the 
    resulting time series is strictly equal to the low-frequency time series within the 
    benchmarking window. Typically, the low-frequency time series is an annual one, 
    unknown for the last year, and the high frequency one is either quarterly or 
    monthly. See "Methodology of quarterly national accounts", Insee Méthodes 
    N°126, by Insee (2012, ISBN:978-2-11-068613-8, <https://www.insee.fr/en/information/2579410>).