IDEAS home Printed from https://ideas.repec.org/a/jss/jstsof/v059i04.html
   My bibliography  Save this article

General Purpose Convolution Algorithm in S4 Classes by Means of FFT

Author

Listed:
  • Ruckdeschel, Peter
  • Kohl, Matthias

Abstract

Object orientation provides a flexible framework for the implementation of the convolution of arbitrary distributions of real-valued random variables. We discuss an algorithm which is based on the fast Fourier transform. It directly applies to lattice-supported distributions. In the case of continuous distributions an additional discretization to a linear lattice is necessary and the resulting lattice-supported distributions are suitably smoothed after convolution. We compare our algorithm to other approaches aiming at a similar generality as to accuracy and speed. In situations where the exact results are known, several checks confirm a high accuracy of the proposed algorithm which is also illustrated for approximations of non-central χ2 distributions. By means of object orientation this default algorithm is overloaded by more specific algorithms where possible, in particular where explicit convolution formulae are available. Our focus is on R package distr which implements this approach, overloading operator + for convolution; based on this convolution, we define a whole arithmetics of mathematical operations acting on distribution objects, comprising operators +, -, *, /, and ^.

Suggested Citation

  • Ruckdeschel, Peter & Kohl, Matthias, 2014. "General Purpose Convolution Algorithm in S4 Classes by Means of FFT," Journal of Statistical Software, Foundation for Open Access Statistics, vol. 59(i04).
    • Handle: RePEc:jss:jstsof:v:059:i04
    DOI: http://hdl.handle.net/10.18637/jss.v059.i04
    as

    Download full text from publisher

    File URL: https://www.jstatsoft.org/index.php/jss/article/view/v059i04/v59i04.pdf
    Download Restriction: no
    File URL: https://www.jstatsoft.org/index.php/jss/article/downloadSuppFile/v059i04/distr_2.5.3.tar.gz
    Download Restriction: no
    File URL: https://www.jstatsoft.org/index.php/jss/article/downloadSuppFile/v059i04/v59i04.R
    Download Restriction: no
    File URL: https://libkey.io/http://hdl.handle.net/10.18637/jss.v059.i04?utm_source=ideas
    LibKey link: if access is restricted and if your library uses this service, LibKey will redirect you to where you can use your library subscription to access this item
    ---><---

    References listed on IDEAS

    as
    1. Warnes, Gregory R., 2002. "HYDRA: a Java library for Markov Chain Monte Carlo," Journal of Statistical Software, Foundation for Open Access Statistics, vol. 7(i04).
    2. Grübel, Rudolf & Hermesmeier, Renate, 1999. "Computation of Compound Distributions I: Aliasing Errors and Exponential Tilting," ASTIN Bulletin, Cambridge University Press, vol. 29(2), pages 197-214, November.
    Full references (including those not matched with items on IDEAS)

    Citations

    Citations are extracted by the CitEc Project, subscribe to its RSS feed for this item.
    as


    Cited by:

    1. Biscarri, William & Zhao, Sihai Dave & Brunner, Robert J., 2018. "A simple and fast method for computing the Poisson binomial distribution function," Computational Statistics & Data Analysis, Elsevier, vol. 122(C), pages 92-100.
    2. Valeriy A. Naumov & Yuliya V. Gaidamaka & Konstantin E. Samouylov, 2020. "Computing the Stationary Distribution of Queueing Systems with Random Resource Requirements via Fast Fourier Transform," Mathematics, MDPI, vol. 8(5), pages 1-9, May.
    3. Jenq-Tzong Shiau, 2021. "Analytical Water Shortage Probabilities and Distributions of Various Lead Times for a Water Supply Reservoir," Water Resources Management: An International Journal, Published for the European Water Resources Association (EWRA), Springer;European Water Resources Association (EWRA), vol. 35(11), pages 3809-3825, September.
    4. Taro Ohdoko & Satoru Komatsu, 2023. "Integrating a Pareto-Distributed Scale into the Mixed Logit Model: A Mathematical Concept," Mathematics, MDPI, vol. 11(23), pages 1-22, November.
    5. Richard L. Warr & Cason J. Wight, 2020. "Error Bounds for Cumulative Distribution Functions of Convolutions via the Discrete Fourier Transform," Methodology and Computing in Applied Probability, Springer, vol. 22(3), pages 881-904, September.

    Most related items

    These are the items that most often cite the same works as this one and are cited by the same works as this one.
    1. Vernic, Raluca, 2018. "On the evaluation of some multivariate compound distributions with Sarmanov’s counting distribution," Insurance: Mathematics and Economics, Elsevier, vol. 79(C), pages 184-193.
    2. Richard L. Warr & Cason J. Wight, 2020. "Error Bounds for Cumulative Distribution Functions of Convolutions via the Discrete Fourier Transform," Methodology and Computing in Applied Probability, Springer, vol. 22(3), pages 881-904, September.
    3. Dominique Guegan & Bertrand Hassani, 2009. "A modified Panjer algorithm for operational risk capital calculations," PSE-Ecole d'économie de Paris (Postprint) halshs-00443846, HAL.
    4. Dominique Guegan & Bertrand Hassani, 2009. "A modified Panjer algorithm for operational risk capital calculations," Post-Print halshs-00443846, HAL.
    5. Xiaolin Luo & Pavel V. Shevchenko, 2009. "Computing Tails of Compound Distributions Using Direct Numerical Integration," Papers 0904.0830, arXiv.org, revised Feb 2010.
    6. Hu, Xiang & Duan, Baige & Zhang, Lianzeng, 2017. "De Vylder approximation to the optimal retention for a combination of quota-share and excess of loss reinsurance with partial information," Insurance: Mathematics and Economics, Elsevier, vol. 76(C), pages 48-55.
    7. Paul Embrechts & Marco Frei, 2009. "Panjer recursion versus FFT for compound distributions," Mathematical Methods of Operations Research, Springer;Gesellschaft für Operations Research (GOR);Nederlands Genootschap voor Besliskunde (NGB), vol. 69(3), pages 497-508, July.
    8. Dominique Guegan & Bertrand Hassani, 2009. "A new algorithm for the loss distribution function with applications to Operational Risk Management," Documents de travail du Centre d'Economie de la Sorbonne 09023, Université Panthéon-Sorbonne (Paris 1), Centre d'Economie de la Sorbonne, revised Nov 2009.
    9. Sangüesa, C., 2008. "Error bounds in approximations of random sums using gamma-type operators," Insurance: Mathematics and Economics, Elsevier, vol. 42(2), pages 484-491, April.
    10. Valeriy A. Naumov & Yuliya V. Gaidamaka & Konstantin E. Samouylov, 2020. "Computing the Stationary Distribution of Queueing Systems with Random Resource Requirements via Fast Fourier Transform," Mathematics, MDPI, vol. 8(5), pages 1-9, May.
    11. Wen Wei & Guan Chen, 2008. "JavaStatSoft: design patterns and features," Computational Statistics, Springer, vol. 23(2), pages 235-251, April.
    12. Dazhi Jiao & Yuzhen Ye & Haixu Tang, 2013. "Probabilistic Inference of Biochemical Reactions in Microbial Communities from Metagenomic Sequences," PLOS Computational Biology, Public Library of Science, vol. 9(3), pages 1-11, March.
    13. Dominique Guegan & Bertrand Hassani, 2009. "A new algorithm for the loss distribution function with applications to Operational Risk Management," Post-Print halshs-00384398, HAL.
    14. Thomas Siller, 2013. "Measuring marginal risk contributions in credit portfolios," Quantitative Finance, Taylor & Francis Journals, vol. 13(12), pages 1915-1923, December.
    15. Lorenzo Cappello & Stephen G. Walker, 2018. "A Bayesian Motivated Laplace Inversion for Multivariate Probability Distributions," Methodology and Computing in Applied Probability, Springer, vol. 20(2), pages 777-797, June.
    16. Li Qin & Susan M. Pitts, 2012. "Nonparametric Estimation of the Finite-Time Survival Probability with Zero Initial Capital in the Classical Risk Model," Methodology and Computing in Applied Probability, Springer, vol. 14(4), pages 919-936, December.
    17. Pavel V. Shevchenko, 2010. "Calculation of aggregate loss distributions," Papers 1008.1108, arXiv.org.
    18. Raluca Vernic, 2018. "On the Evaluation of the Distribution of a General Multivariate Collective Model: Recursions versus Fast Fourier Transform," Risks, MDPI, vol. 6(3), pages 1-14, August.

    More about this item

    Statistics

    Access and download statistics

    Corrections

    All material on this site has been provided by the respective publishers and authors. You can help correct errors and omissions. When requesting a correction, please mention this item's handle: RePEc:jss:jstsof:v:059:i04. See general information about how to correct material in RePEc.

    If you have authored this item and are not yet registered with RePEc, we encourage you to do it here. This allows to link your profile to this item. It also allows you to accept potential citations to this item that we are uncertain about.

    If CitEc recognized a bibliographic reference but did not link an item in RePEc to it, you can help with this form .

    If you know of missing items citing this one, you can help us creating those links by adding the relevant references in the same way as above, for each refering item. If you are a registered author of this item, you may also want to check the "citations" tab in your RePEc Author Service profile, as there may be some citations waiting for confirmation.

    For technical questions regarding this item, or to correct its authors, title, abstract, bibliographic or download information, contact: Christopher F. Baum (email available below). General contact details of provider: http://www.jstatsoft.org/ .

    Please note that corrections may take a couple of weeks to filter through the various RePEc services.

    IDEAS is a RePEc service. RePEc uses bibliographic data supplied by the respective publishers.