Skip to content

Audit/v200#10

Open
tfm000 wants to merge 101 commits intomainfrom
audit/v200
Open

Audit/v200#10
tfm000 wants to merge 101 commits intomainfrom
audit/v200

Conversation

@tfm000
Copy link
Copy Markdown
Owner

@tfm000 tfm000 commented Mar 25, 2026

  • Implement rigorous tests for Bessel K_v, inverse incomplete gamma functions, and Student's t CDF in test_special.py.
  • Validate against scipy.special and scipy.stats with tight tolerances to catch known errors.
  • Introduce tests for sample skewness and kurtosis in test_stats.py, ensuring consistency with scipy.stats.
  • Create extensive tests for all 11 univariate distributions in test_univariate.py, covering logpdf, CDF, stats, fitting, and edge cases.
  • Include checks for PDF integration, inverse consistency, parameter recovery, and gradient correctness.
  • Address specific findings related to Student-T and IG distributions to ensure accuracy and reliability.

tfm000 added 30 commits March 24, 2026 22:47
- Implement rigorous tests for Bessel K_v, inverse incomplete gamma functions, and Student's t CDF in `test_special.py`.
- Validate against scipy.special and scipy.stats with tight tolerances to catch known errors.
- Introduce tests for sample skewness and kurtosis in `test_stats.py`, ensuring consistency with scipy.stats.
- Create extensive tests for all 11 univariate distributions in `test_univariate.py`, covering logpdf, CDF, stats, fitting, and edge cases.
- Include checks for PDF integration, inverse consistency, parameter recovery, and gradient correctness.
- Address specific findings related to Student-T and IG distributions to ensure accuracy and reliability.
… precision and stability; enhance initial approximation and Newton-Halley refinement process.
…ts for extreme values and refine tolerance levels for better precision.
… handling; add explicit normalization for total PDF integral to ensure accurate results.
…initial estimates; replace kv with log_kv for stability and add sample moments computation.
…g_kv for stability, add expected value computations, and implement negative log-likelihood gradient calculation.
…likelihood computation and update fit method to support EM approach.
…ctor _fit_mle and _fit_ldmle to use sample moments for better convergence.
… _sample_moments method and update _fit_mle to use these estimates for improved fitting.
…re functions for improved parameter estimation.
… enhance fit method to support MOM and MLE options.
tfm000 and others added 30 commits April 16, 2026 18:17
…, enhance clarity, and fix p-value calculations.
…ding references, and enhance clarity in assertions.
- Implemented `_cdf_grid_piecewise` to compute CDF values at a sorted grid of x values using a fixed 16-point Gauss-Legendre quadrature, improving performance for dense grids.
- Updated `_cdf` to utilize the new piecewise CDF computation method.
- Enhanced `_cubic_ppf_solve` to dispatch between numerical and closed-form CDF evaluations, leveraging `_cdf_grid_piecewise` for numerical distributions.
- Introduced `warm_bounds` parameter in PPF functions to optimize bracket discovery during iterative fitting loops, reducing computational overhead.
- Introduced `DataScaler` class for affine rescaling of input data with methods: z-score, min-max, robust, and max-abs.
- Implemented `fit`, `transform`, and `inverse_transform` methods, ensuring compatibility with JAX's JIT, grad, and vmap.
- Added comprehensive tests covering various aspects including round-trip identity, method-specific statistics, new-data stability, zero-variance safety, and integration with `copulax.univariate.Normal`.
- Created documentation for preprocessing module and DataScaler class.
- Introduced `gh_reference_data.py` containing auto-generated reference data for generalized hyperbolic and skewed-t distributions.
- Implemented `test_multivariate_reference.py` to validate logpdf and pdf of `mvt_gh` and `mvt_skewed_t` against the R ghyp package.
- Enhanced `test_copulas_elliptical.py` with manual Sklar verification for Student-t copula.
- Updated `test_univariate.py` to include edge case tests for logpdf and cdf behavior at support boundaries.
- Added tests for skewed-t distribution against scipy's genhyperbolic limit.
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment

Labels

None yet

Projects

None yet

Development

Successfully merging this pull request may close these issues.

1 participant