Package: EstemPMM 0.4.0

EstemPMM: Polynomial Maximization Method for Non-Gaussian Regression

Implements the Polynomial Maximization Method ('PMM') for parameter estimation in linear and time series models when error distributions deviate from normality. The 'PMM2' variant achieves lower variance parameter estimates compared to ordinary least squares ('OLS') when errors exhibit significant skewness. The 'PMM3' variant (S=3) targets symmetric platykurtic error distributions, reducing variance when excess kurtosis is negative. Includes automatic method selection ('pmm_dispatch'), linear regression, 'AR'/'MA'/'ARMA'/'ARIMA' models, and bootstrap inference. Methodology described in Zabolotnii, Warsza, and Tkachenko (2018) <doi:10.1007/978-3-319-77179-3_75>, Zabolotnii, Tkachenko, and Warsza (2022) <doi:10.1007/978-3-031-03502-9_37>, and Zabolotnii, Tkachenko, and Warsza (2023) <doi:10.1007/978-3-031-25844-2_21>, and Zabolotnii (2025) <doi:10.48550/arXiv.2511.07059>.

Authors:Serhii Zabolotnii [aut, cre]

EstemPMM_0.4.0.tar.gz
EstemPMM_0.4.0.zip(r-4.7)EstemPMM_0.4.0.zip(r-4.6)EstemPMM_0.4.0.zip(r-4.5)
EstemPMM_0.4.0.tgz(r-4.6-any)EstemPMM_0.4.0.tgz(r-4.5-any)
EstemPMM_0.4.0.tar.gz(r-4.7-any)EstemPMM_0.4.0.tar.gz(r-4.6-any)
EstemPMM_0.4.0.tgz(r-4.6-emscripten)
manual.pdf |manual.html
DESCRIPTION |NEWS
card.svg |card.png
EstemPMM/json (API)

# Install 'EstemPMM' in R:
install.packages('EstemPMM', repos = c('https://szabolotnii.r-universe.dev', 'https://cloud.r-project.org'))

Bug tracker:https://github.com/szabolotnii/estempmm/issues

Datasets:

On CRAN:

Conda:

6.23 score 57 scripts 465 downloads 58 exports 0 dependencies

Last updated from:718925a6dd. Checks:9 OK. Indexed: yes.

TargetResultTimeFilesSyslog
linux-devel-x86_64OK148
source / vignettesOK200
linux-release-x86_64OK148
macos-release-arm64OK177
macos-oldrel-arm64OK153
windows-develOK114
windows-releaseOK128
windows-oldrelOK104
wasm-releaseOK100

Exports:ar_pmm2ar_pmm3arima_pmm2arima_pmm3arma_pmm2arma_pmm3coefcompare_ar_methodscompare_arima_methodscompare_arma_methodscompare_ma_methodscompare_sar_methodscompare_ts_methodscompare_with_olscompute_momentscompute_moments_pmm3confintcreate_sar_matrixcreate_sarma_matrixfittedlm_pmm2lm_pmm3ma_pmm2ma_pmm3nobsplotplot_pmm2_bootstrappmm_arpmm_arimapmm_armapmm_dispatchpmm_gamma6pmm_kurtosispmm_lmpmm_mapmm_sarimapmm_skewnesspmm2_inferencepmm2_monte_carlo_comparepmm2_nonlinear_iterativepmm2_nonlinear_onesteppmm2_variance_factorpmm2_variance_matricespmm3_variance_factorpmm3_variance_matricespredictresidualssar_pmm2sarima_pmm2sarma_pmm2showsma_pmm2summarytest_symmetryts_pmm2ts_pmm2_inferencets_pmm3vcov

Dependencies:

Introduction to PMM2: Polynomial Maximization Method
The Problem: When OLS Falls Short | Practical Implications | The Solution: Polynomial Maximization Method (PMM2) | Theoretical Foundation | Getting Started with EstemPMM | Example 1: Simple Linear Regression with Skewed Errors | Data Generation | Visualize Error Distribution | Model Fitting: PMM2 vs OLS | Bootstrap Inference for Variance Comparison | Example 2: Multiple Regression | Example 3: Polynomial Regression | Diagnostic Plots | Comparing PMM2 with OLS: Monte Carlo Evidence | When to Use PMM2 | Practical Recommendations | Next Steps | Summary | References

Last update: 2026-05-14
Started: 2025-10-23

PMM3: Linear Regression for Symmetric Platykurtic Errors
The Problem: When OLS and PMM2 Are Not Enough | Where Platykurtic Errors Arise | The Solution: PMM3 (Polynomial Maximization Method, S=3) | The Newton-Raphson Algorithm | Getting Started | Example 1: Simple Regression with Uniform Errors | Data Generation | Visualize Error Distribution | Automatic Method Selection | Fit PMM3 | Compare with OLS | Example 2: Three-Way Comparison (OLS vs PMM2 vs PMM3) | Example 3: Real Data Diagnostic — Auto MPG Horsepower | Fit Quadratic Model | Dispatcher Recommendation | Visualize the OLS Fit | Example 4: Real Data — Auto MPG (Linear, PMM2 Case) | Example 5: Multiple Regression | Example 6: No-Intercept Model | Interpreting the PMM3 Summary | Adaptive Mode | Diagnostic Plots | Utility Functions | Testing Symmetry | Computing Moments | Sixth-Order Cumulant | Variance Factor | When to Use PMM3 | Decision Workflow | Efficiency Summary (Monte Carlo Evidence) | Next Steps | References

Last update: 2026-05-14
Started: 2026-03-20

Seasonal Time Series Models with PMM2
Introduction | Load Package | Part 1: Seasonal Autoregressive (SAR) Models | Understanding SAR Models | Example: SAR(1,1)_12 with Monthly Data | Seasonal Patterns Visualization | Fitting SAR Models: PMM2 vs Classical Methods | Understanding PMM2 Efficiency | Residual Diagnostics | Part 2: Seasonal Moving Average (SMA) Models | Understanding SMA Models | Example: SMA(1)_4 with Quarterly Data | Quarterly Seasonal Pattern | Fitting SMA Models | SMA Convergence Behavior | Part 3: Monte Carlo Evidence for Seasonal Models | SAR Monte Carlo Comparison | Part 4: SARMA and SARIMA Models | SARMA(p,P,q,Q)_s Models | SARIMA Models with Differencing | Part 5: Model Comparison and Selection | Comparing Different Methods | Information Criteria for Model Selection | Part 6: Real-World Example - Airline Passengers | Modeling Strategy | Part 7: Practical Guidelines | When to Use Seasonal Models | Choosing Between Seasonal Models | PMM2 vs Classical Methods | Diagnostic Checklist | Part 8: Advanced Topics | Multiple Seasonal Periods | Seasonal Adjustment | Bootstrap Inference for Seasonal Models | Summary | References | Next Steps

Last update: 2026-05-14
Started: 2026-05-14

PMM3 for Time Series: AR, MA, ARMA, and ARIMA Models
Introduction | When PMM3 vs PMM2 vs OLS? | Setup | Part 1: AR Models with PMM3 | AR(1) with Uniform Innovations | Fit AR(1): PMM3 vs MLE | AR(2) Model | Part 2: MA Models with PMM3 | MA(1) Model | MA(2) Model | Part 3: ARMA Models | ARMA(1,1) Model | Diagnostic Plots | Part 4: ARIMA Models with Differencing | ARIMA(1,1,0) Model | Fit ARIMA(1,1,0) | ARIMA(1,1,1) Model | Part 5: Forecasting with PMM3 Models | Forecasting ARMA Models | Part 6: Method Selection with pmm_dispatch() | Part 7: Real-Data Example — DJIA 2002 | Analyze Innovation Properties | Fit AR(1) with All Three Methods | Part 8: Real-Data Example — WTI Crude Oil Prices | Innovation Analysis | Fit and Compare | Part 9: Adaptive Mode | Part 10: Model Comparison via AIC | Practical Guidelines | When to Use PMM3 for Time Series | Diagnostic Workflow | Computational Notes | Efficiency Summary | Available Functions | Next Steps | References

Last update: 2026-03-20
Started: 2026-03-20

PMM2 for Time Series: AR, MA, ARMA, ARIMA, and Seasonal Models
Introduction | Setup | Part 1: Autoregressive (AR) Models | AR(1) Model with Skewed Innovations | Estimate AR(1) Model: PMM2 vs CSS | AR(2) Model | Part 2: Moving Average (MA) Models | MA(1) Model | MA(2) Model | Part 3: ARMA Models | ARMA(1,1) Model | Diagnostic Plots for ARMA Models | Part 4: ARIMA Models with Differencing | ARIMA(1,1,1) Model | Fitting ARIMA Models | Part 5: Bootstrap Inference for Time Series | Part 6: Comparing Methods with Monte Carlo | Part 7: Prediction with Time Series Models | Practical Guidelines | When to Use PMM2 for Time Series | Model Selection Strategy | Computational Considerations | Summary Statistics: Efficiency Gains | Advanced Topics | Custom Innovation Distributions | Part 8: Seasonal Models (SAR and SMA) | Seasonal Autoregressive (SAR) Models | Seasonal Moving Average (SMA) Models | Comparing Seasonal Methods | Practical Applications of Seasonal Models | Conclusion | Next Steps | References

Last update: 2025-11-13
Started: 2025-10-23

Bootstrap Inference for PMM2 Models
Introduction | Bootstrap Basics | The Bootstrap Principle | Setup | Part 1: Bootstrap for Linear Regression | Basic Example: Simple Linear Regression | Understanding Bootstrap Distributions | Part 2: Confidence Interval Methods | Percentile Method | Bias-Corrected (BC) Method | Comparison with OLS | Part 3: Hypothesis Testing | Testing Individual Coefficients | Testing Equality of Coefficients | Part 4: Bootstrap for Time Series Models | AR Model Bootstrap | ARMA Model Bootstrap | Part 5: Choosing Bootstrap Parameters | Number of Bootstrap Samples (B) | Parallel Bootstrap | Part 6: Diagnostic Checks | Convergence of Bootstrap Estimates | Bootstrap Distribution Visualization | Part 7: Practical Guidelines | When to Use Bootstrap for PMM2 | Common Pitfalls and Solutions | Part 8: Advanced Topics | Studentized Bootstrap | Block Bootstrap for Time Series | Summary | Best Practices Checklist | References | Next Steps

Last update: 2025-10-23
Started: 2025-10-23

Readme and manuals

Help Manual

Help pageTopics
Fit an AR model using PMM2 (wrapper)ar_pmm2
Fit an AR model using PMM3ar_pmm3
Fit an ARIMA model using PMM2 (wrapper)arima_pmm2
Fit an ARIMA model using PMM3arima_pmm3
S4 class for storing PMM2 ARIMA model resultsARIMAPMM2-class
S4 class for PMM3 ARIMA model resultsARIMAPMM3-class
Fit an ARMA model using PMM2 (wrapper)arma_pmm2
Fit an ARMA model using PMM3arma_pmm3
S4 class for storing PMM2 ARMA model resultsARMAPMM2-class
S4 class for PMM3 ARMA model resultsARMAPMM3-class
S4 class for storing PMM2 AR model resultsARPMM2-class
S4 class for PMM3 AR model resultsARPMM3-class
Auto MPG Datasetauto_mpg
Virtual S4 class for the PMM2 model familyBasePMM2-class
Virtual S4 class for the PMM3 model familyBasePMM3-class
Extract coefficients from PMM2fit objectcoef,PMM2fit-method
Extract coefficients from PMM3fit objectcoef,PMM3fit-method
Extract coefficients from SARPMM2 objectcoef,SARPMM2-method
Extract coefficients from SMAPMM2 objectcoef,SMAPMM2-method
Extract coefficients from TS2fit objectcoef,TS2fit-method
Extract coefficients from TS3fit objectcoef,TS3fit-method
Compare AR methodscompare_ar_methods
Compare ARIMA methodscompare_arima_methods
Compare ARMA methodscompare_arma_methods
Compare MA methodscompare_ma_methods
Compare SAR model estimation methodscompare_sar_methods
Compare PMM2 with classical time series estimation methodscompare_ts_methods
Compare PMM2 with OLScompare_with_ols
Calculate moments and cumulants of error distributioncompute_moments
Compute central moments for PMM3 from residualscompute_moments_pmm3
Compute PMM2 weights and componentscompute_pmm2_components
Confidence intervals for PMM2fit coefficientsconfint,PMM2fit-method
Confidence intervals for PMM3fit coefficientsconfint,PMM3fit-method
Confidence intervals for TS2fit AR model coefficientsconfint,TS2fit-method
Confidence intervals for TS3fit AR model coefficientsconfint,TS3fit-method
Create design matrix for seasonal AR modelcreate_sar_matrix
Create design matrix for seasonal ARMA modelcreate_sarma_matrix
WTI Crude Oil PricesDCOILWTICO
Dow Jones Industrial Average Daily Data (July-December 2002)djia2002
Extract fitted values from PMM2fit objectfitted,PMM2fit-method
Extract fitted values from PMM3fit objectfitted,PMM3fit-method
Extract fitted values from TS2fit objectfitted,TS2fit-method
Extract fitted values from TS3fit objectfitted,TS3fit-method
Format method for PMMdispatch objectsformat.PMMdispatch
Calculate SARIMAX Jacobian (Numerical)get_sarimax_jacobian
Calculate SARIMAX Residualsget_sarimax_residuals
PMM2: Main function for PMM2 (S=2)lm_pmm2
PMM3: Fit linear model using Polynomial Maximization Method (S=3)lm_pmm3
Extract log-likelihood from PMM2fit objectlogLik.PMM2fit
Extract log-likelihood from PMM3fit objectlogLik.PMM3fit
Extract log-likelihood from TS2fit objectlogLik.TS2fit
Extract log-likelihood from TS3fit objectlogLik.TS3fit
Fit an MA model using PMM2 (wrapper)ma_pmm2
Fit an MA model using PMM3ma_pmm3
S4 class for storing PMM2 MA model resultsMAPMM2-class
S4 class for PMM3 MA model resultsMAPMM3-class
Number of observations in PMM2fit objectnobs,PMM2fit-method
Number of observations in PMM3fit objectnobs,PMM3fit-method
Number of observations in TS2fit objectnobs,TS2fit-method
Number of observations in TS3fit objectnobs,TS3fit-method
Plot bootstrap distributions for PMM2 fitplot_pmm2_bootstrap
Plot diagnostic plots for PMM2fit objectplot,PMM2fit,missing-method
Plot diagnostic plots for PMM3fit objectplot,PMM3fit,missing-method
Build diagnostic plots for TS2fit objectsplot,TS2fit,missing-method
Plot diagnostic plots for TS3fit objectplot,TS3fit,missing-method
Fit an AR model with the polynomial maximization methodpmm_ar
Fit an ARIMA model with the polynomial maximization methodpmm_arima
Fit an ARMA model with the polynomial maximization methodpmm_arma
Automatic PMM method selectionpmm_dispatch
Compute sixth-order cumulant coefficient gamma6pmm_gamma6
Calculate kurtosis from datapmm_kurtosis
Fit a linear model with the polynomial maximization methodpmm_lm
Fit an MA model with the polynomial maximization methodpmm_ma
Fit a seasonal ARIMA model with the polynomial maximization methodpmm_sarima
Calculate skewness from datapmm_skewness
Bootstrap inference for PMM2 fitpmm2_inference
Monte Carlo comparison of PMM2 estimation methodspmm2_monte_carlo_compare
Universal PMM2 estimator (Iterative)pmm2_nonlinear_iterative
Universal PMM2 estimator (One-step Global)pmm2_nonlinear_onestep
Calculate theoretical skewness, kurtosis coefficients and variance reduction factorpmm2_variance_factor
Calculate theoretical variance matrices for OLS and PMM2pmm2_variance_matrices
S4 class for storing PMM2 regression model resultsPMM2fit-class
Calculate PMM3 theoretical variance reduction factorpmm3_variance_factor
Calculate theoretical variance matrices for OLS and PMM3pmm3_variance_matrices
S4 class for PMM3 regression fit resultsPMM3fit-class
Virtual root S4 class for PMM fit objectsPMMfit-class
Virtual S4 class for PMM time-series fit objectsPMMtsfit-class
Prediction method for PMM2fit objectspredict,PMM2fit-method
Predict method for PMM3fit objectspredict,PMM3fit-method
Prediction method for TS2fit objectspredict,TS2fit-method
Predict method for TS3fit objectspredict,TS3fit-method
Print method for PMMdispatch objectsprint.PMMdispatch
Extract residuals from PMM2fit objectresiduals,PMM2fit-method
Extract residuals from PMM3fit objectresiduals,PMM3fit-method
Extract residuals from TS2fit objectresiduals,TS2fit-method
Extract residuals from TS3fit objectresiduals,TS3fit-method
Fit Seasonal AR model using PMM2 methodsar_pmm2
Fit a Seasonal ARIMA model using PMM2 methodsarima_pmm2
S4 class for Seasonal ARIMA model results with PMM2SARIMAPMM2-class
Fit a Seasonal ARMA model using PMM2 methodsarma_pmm2
S4 class for Seasonal ARMA model results with PMM2SARMAPMM2-class
S4 class for Seasonal AR model results with PMM2SARPMM2-class
Show method for PMM2fit objectsshow,PMM2fit-method
Show method for PMM3fit objectsshow,PMM3fit-method
Show method for TS2fit objects (and subclasses)show,TS2fit-method
Show method for TS3fit objects (and subclasses)show,TS3fit-method
Fit a Seasonal MA model using PMM2sma_pmm2
S4 class for Seasonal MA PMM2 resultsSMAPMM2-class
PMM2 step solversolve_pmm2_step
Generic summary method for PMM2fit objectssummary,PMM2fit-method
Summary method for PMM3fit objectssummary,PMM3fit-method
Generic summary method for SARIMAPMM2 objectssummary,SARIMAPMM2-method
Generic summary method for SARMAPMM2 objectssummary,SARMAPMM2-method
Summary method for SARPMM2 objectssummary,SARPMM2-method
Summary method for SMAPMM2 objectssummary,SMAPMM2-method
Generic summary method for TS2fit objectssummary,TS2fit-method
Summary method for TS3fit objectssummary,TS3fit-method
Summary method for PMMdispatch objectssummary.PMMdispatch
Test whether residuals are sufficiently symmetric for PMM3test_symmetry
Fit a time series model using the PMM2 methodts_pmm2
Bootstrap inference for PMM2 time series modelsts_pmm2_inference
Fit a time series model using PMM3ts_pmm3
S4 class for PMM2 time-series fit resultsTS2fit-class
S4 class for PMM3 time-series fit resultsTS3fit-class
Variance-covariance matrix for PMM2fit objectvcov,PMM2fit-method
Variance-covariance matrix for PMM3fit objectsvcov,PMM3fit-method
Variance-covariance matrix for TS2fit AR modelsvcov,TS2fit-method
Variance-covariance matrix for TS3fit AR modelsvcov,TS3fit-method