Sampling large hyperplane-truncated multivariate normal distributions - FAYOL / DEMO : Décision en Entreprise : Modélisation, Optimisation Access content directly
Journal Articles Computational Statistics Year : 2023

Sampling large hyperplane-truncated multivariate normal distributions

Abstract

Generating multivariate normal distributions is widely used in various fields, including engineering, statistics, finance and machine learning. In this paper, simulating large multivariate normal distributions truncated on the intersection of a set of hyperplanes is investigated. Specifically, the proposed methodology focuses on cases where the prior multivariate normal is extracted from a stationary Gaussian process (GP). It is based on combining both Karhunen-Loève expansions (KLE) and Matheron’s update rules (MUR). The KLE requires the computation of the decomposition of the covariance matrix of the random variables, which can become expensive when the random vector is too large. To address this issue, the input domain is split in smallest subdomains where the eigendecomposition can be computed. Due to the stationary property, only the eigendecomposition of the first subdomain is required. Through this strategy, the computational complexity is drastically reduced. The mean-square truncation and block errors have been calculated. The efficiency of the proposed approach has been demonstrated through both synthetic and real data studies.
Fichier principal
Vignette du fichier
highdimMVN.pdf (664.88 Ko) Télécharger le fichier
Origin : Files produced by the author(s)

Dates and versions

hal-03741860 , version 1 (02-08-2022)
hal-03741860 , version 2 (19-09-2023)

Identifiers

Cite

Hassan Maatouk, Didier Rullière, Xavier Bay. Sampling large hyperplane-truncated multivariate normal distributions. Computational Statistics, 2023, ⟨10.1007/s00180-023-01416-7⟩. ⟨hal-03741860v2⟩
139 View
80 Download

Altmetric

Share

Gmail Facebook X LinkedIn More