Publica
Hier finden Sie wissenschaftliche Publikationen aus den FraunhoferInstituten. General purpose convolution algorithm in S4 classes by means of FFT
 Journal of statistical software 59 (2014), Nr.4, 25 S. ISSN: 15487660 

 Englisch 
 Zeitschriftenaufsatz, Elektronische Publikation 
 Fraunhofer ITWM () 
Abstract
Object orientation provides a flexible framework for the implementation of the convolution of arbitrary distributions of realvalued random variables. We discuss an algorithm which is based on the fast Fourier transform. It directly applies to latticesupported distributions. In the case of continuous distributions an additional discretization to a linear lattice is necessary and the resulting latticesupported 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 noncentral χ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 .