Multivariate Subgaussian Stable Distributions in R

We introduce and showcase mvpd (an acronym for multivariate product distributions), a package that uses a product distribution approach to calculating multivariate subgaussian stable distribution functions. The family of multivariate subgaussian stable distributions are elliptically contoured multivariate stable distributions that contain the multivariate Cauchy and the multivariate normal distribution. These distributions can be useful in modeling data and phenomena that have heavier tails than the normal distribution (more frequent occurrence of extreme values). Application areas include log returns for stocks, signal processing for radar and sonar data, astronomy, and hunting patterns of sharks.

Bruce J. Swihart (National Institutes of Health) , John P. Nolan (American University)

Supplementary materials

Supplementary materials are available in addition to this article. It can be downloaded at



Text and figures are licensed under Creative Commons Attribution CC BY 4.0. The figures that have been reused from other sources don't fall under this license and can be recognized by a note in their caption: "Figure from ...".


For attribution, please cite this work as

Swihart & Nolan, "The R Journal: Multivariate Subgaussian Stable Distributions in R", The R Journal, 2022

BibTeX citation

  author = {Swihart, Bruce J. and Nolan, John P.},
  title = {The R Journal: Multivariate Subgaussian Stable Distributions in R},
  journal = {The R Journal},
  year = {2022},
  note = {},
  doi = {10.32614/RJ-2022-056},
  volume = {14},
  issue = {3},
  issn = {2073-4859},
  pages = {293-302}