---
title: Mean-Variance Mixture of Normals
url: https://www.emergentmind.com/topics/mean-variance-mixture-of-normals-mvmn
type: topic
---

# Mean-Variance Mixture of Normals

Mean-variance mixture of normals (MVMN), also called the location-scale mixture of normals, is a class of distributions obtained by letting a positive latent scalar simultaneously shift the mean and scale the covariance of a multivariate normal vector. In its standard form, if \(Z\sim N_p(0,I_p)\), \(W>0\) is independent of \(Z\), \(\mu\in\mathbb R^p\) is a location vector, \(\delta\in\mathbb R^p\) is a shape or skewness vector, and \(\Sigma\) is positive definite, then
\[
Y=\mu+W\delta+\sqrt{W}\,\Sigma^{1/2}Z.
\]
Equivalently,
\[
Y\mid W=w\sim N_p(\mu+w\delta,\;w\Sigma),
\]
with marginal density
\[
f(y;\mu,\Sigma,\delta;h(w;\theta))=\int_0^\infty \phi_p(y;\mu+w\delta,w\Sigma)\,h(w;\theta)\,dw.
\]
This construction extends variance mixtures by allowing the same latent variable to govern both mean and covariance, thereby accommodating skewness, heavy tails, and multivariate dependence within a single hierarchical mechanism [2005.06883].

## 1. Formal construction and equivalent representations

The defining feature of MVMN is that the same positive scalar mixing variable controls both conditional location and conditional scale. The notation
\[
Y\sim MVMN_p(\mu,\Sigma,\delta;h(w;\theta))
\]
is used when \(Y\) has the density above. In finance-oriented notation, the same class is often written as
\[
X\overset{d}{=}\mu+\gamma Z+\sqrt{Z}\,A N_d,
\]
where \(A=\Sigma^{1/2}\), \(N_d\sim N_d(0,I_d)\), and \(Z\) is positive and independent of \(N_d\); conditionally on \(Z=z\), one has
\[
X\mid (Z=z)\sim N_d(\mu+z\gamma,\;z\Sigma).
\]
This is the same latent-Gaussian structure with a different symbol set for the mixing variable and skewness vector [2202.02488].

A broader continuous-mixture formulation writes
\[
Y=\xi+R\gamma+SX,
\qquad X\sim N_d(0,\Sigma),\quad (R,S)\perp X,
\]
or, more generally,
\[
Y=\xi+r(U,V)\gamma+s(U,V)X,
\]
with two independent univariate random variables \(U\) and \(V\). The classical mean-variance mixture corresponds to the special choice \(R=V\) and \(S=V^{1/2}\), that is,
\[
Y=\xi+V\gamma+V^{1/2}X.
\]
This general framework places mean mixtures, variance mixtures, and mean-variance mixtures inside a common family of continuous mixtures of multivariate normals [2003.13076].

The same representation underlies several application-specific formulations. In risk and insurance, an \(n\)-dimensional normal mean-variance mixture is written as
\[
\mathbf X\overset{d}{=}\mathbf m(\Theta)+\sqrt{\Theta}\,A\mathbf Z,
\qquad \mathbf m(\Theta)=\boldsymbol\mu+\Theta\boldsymbol\gamma,
\]
so that
\[
\mathbf X\mid \Theta=\theta \sim MVN_n(\boldsymbol\mu+\theta\boldsymbol\gamma,\;\theta\Sigma).
\]
Because linear combinations preserve the class, aggregate losses \(S=X_1+\cdots+X_n\) are again normal mean-variance mixtures [2601.00568].

## 2. Distinction from mean mixtures and variance mixtures

MVMN sits between two simpler normal-mixture constructions. In a variance mixture of normals (VMN),
\[
Y=\mu+\sqrt{W}\,\Sigma^{1/2}Z,
\qquad
Y\mid W=w\sim N_p(\mu,w\Sigma),
\]
so the mean is constant and only the covariance is mixed. In a mean mixture of normals (MMN),
\[
Y=\mu+W\delta+\Sigma^{1/2}Z,
\qquad
Y\mid W=w\sim N_p(\mu+w\delta,\Sigma),
\]
so the covariance is fixed and only the mean is mixed. MVMN combines both mechanisms:
\[
Y\mid W=w\sim N_p(\mu+w\delta,\;w\Sigma).
\]
Accordingly, MVMN reduces to VMN when \(\delta=0\), but MMN is not a special case of MVMN. The skew-normal distribution belongs to the MMN framework rather than to MVMN [2005.06883].

The basic moment structure makes the dual role of the mixing variable explicit. If \(Y\sim MVMN_p(\mu,\Sigma,\delta;h(w;\theta))\), then
\[
E(Y)=\mu+E(W)\delta,
\]
and
\[
\operatorname{cov}(Y)=\operatorname{var}(W)\,\delta\delta^\top+E(W)\Sigma.
\]
The moment generating function is
\[
M_Y(t)=e^{\mu^\top t}\,M_W\!\left(\delta^\top t+\frac12 t^\top\Sigma t\right),
\]
showing that the MVMN mgf is inherited from the mgf of the mixing variable. In the broader \(Y=\xi+R\gamma+SX\) formulation, the same decomposition appears as
\[
E[Y]=\xi+E[R]\gamma,
\qquad
\operatorname{var}(Y)=\operatorname{var}(R)\gamma\gamma^\top+E[S^2]\Sigma.
\]
These formulas isolate mean mixing through \(R\gamma\) and variance mixing through \(S^2\Sigma\) [2005.06883; 2003.13076].

The family retains several normal-like closure properties. Affine transformations preserve the class:
\[
AY+b\sim MVMN_q(A\mu+b,\;A\Sigma A^\top,\;A\delta;\;h(w;\theta)).
\]
Linear combinations with independent normal vectors remain MVMN; marginals are again MVMN; and conditionals are again MVMN with the expected partitioned parameters. This is one reason the family is used when Gaussian tractability is desired but Gaussian symmetry and light tails are inadequate [2005.06883].

Because the same positive latent variable drives both shift and spread, asymmetry arises when \(\delta\neq 0\), while heavy tails arise through random variance scaling. In this sense MVMN is designed to model skewness and heavy tails simultaneously [2005.06883].

## 3. Canonical families, special cases, and extensions

The generalized hyperbolic (GH) distribution is the most prominent special case of MVMN. It is obtained when the mixing variable follows a generalized inverse Gaussian law,
\[
W\sim GIG(\psi,\chi,\lambda).
\]
In the alternative parameterization used in related work,
\[
V\sim GIG(\lambda,\chi,\psi),
\qquad
Y=\xi+V\gamma+V^{1/2}X,
\]
which yields the multivariate GH distribution. GH is the canonical multivariate mean-variance mixture family in this literature [2005.06883; 2003.13076].

Several important distributions arise as GH subfamilies or limits. The literature summarized here identifies the normal inverse Gaussian, variance gamma, and asymmetric Laplace as notable special cases of GH, and also states that GH encompasses Student \(t\), Laplace, and hyperbolic distributions [2005.06883; 1106.2333]. The multivariate skewed variance gamma (MSVG) model is a gamma-mixed instance with
\[
\mathbf y_i\mid \lambda_i\sim \mathcal N_d(\boldsymbol\mu+\boldsymbol\gamma\lambda_i,\;\lambda_i\boldsymbol\Sigma),
\qquad
\lambda_i\sim \mathcal G(\nu,\nu),
\]
which gives
\[
E(\mathbf Y_i)=\boldsymbol\mu+\boldsymbol\gamma,
\qquad
\operatorname{Cov}(\mathbf Y_i)=\boldsymbol\Sigma+\frac1\nu\boldsymbol\gamma\boldsymbol\gamma'.
\]
This places MSVG squarely inside the normal mean-variance mixture framework [1504.01239].

Less familiar examples noted in the survey literature include MVN Birnbaum–Saunders (MVNBS), where \(W\) is Birnbaum–Saunders, and MVN Lindley (MVNL), where \(W\) has a Lindley distribution [2005.06883]. At the opposite extreme, if the mixing law degenerates at \(W=1\), the model collapses to the multivariate normal [2005.06883].

The MVMN architecture has also been generalized by replacing the Gaussian conditional kernel. The multivariate Mixed Tempered Stable model keeps the same mean-variance mixing form,
\[
Y=\mu+\beta V+\sqrt V\,X,
\]
but replaces the normal innovation with a standardized Classical Tempered Stable innovation. When \(\alpha=2\), Mixed Tempered Stable reduces to a Normal Variance Mean Mixture; with Gamma mixing and \(\alpha=2\), one recovers the Variance Gamma model [1609.00926]. This suggests that classical MVMN occupies a central position inside a wider family of latent scale-and-shift mixture constructions.

## 4. Estimation, identifiability, and inferential issues

For parametric MVMN models, the EM algorithm is the standard estimation tool. The conditional Gaussian representation
\[
Y\mid W=w\sim N_p(\mu+w\delta,\;w\Sigma),\qquad W\sim h(w;\theta)
\]
treats \(W\) as missing data and supports likelihood-based estimation by iterating conditional expectation and maximization steps. The survey literature identifies explicit EM implementations for GH distributions, MVNBS, and MVNL, even though it does not derive a full generic MVMN EM scheme in closed form [2005.06883].

Concrete ECM and ECME schemes have been developed for specific MVMN families. For the MSVG model, the complete-data likelihood splits into a normal part and a gamma part; the conditional law of \(\lambda_i\mid \mathbf y_i\) is generalized inverse Gaussian, and closed-form CM updates are available for \(\boldsymbol\mu\), \(\boldsymbol\gamma\), and \(\boldsymbol\Sigma\), while \(\nu\) is updated by Newton–Raphson or by direct maximization of the observed likelihood. The same paper extends the mean to include autoregressive terms, proposes a delta-region bounding device when the density is unbounded for \(\nu\le d/2\), inserts an extra E-step before updating \(\boldsymbol\Sigma\), and computes standard errors by Louis’s method [1504.01239].

A different inferential direction is semiparametric estimation. In the univariate model
\[
X\stackrel{d}{=}\mu\xi+\sqrt{\xi}\,\eta,
\qquad \eta\sim N(0,1),\ \xi\sim G,
\]
the coefficient \(\mu\) can be identified as the unique zero of a monotone functional \(W(\rho)\), estimated by the empirical root
\[
\widehat\mu_n:=\inf\{\rho>0:W_n(\rho)=0\}\wedge M,
\]
and then the unknown mixing density can be recovered nonparametrically by Mellin-transform inversion. The paper proves an \(n^{-1/2}\) convergence rate for \(\widehat\mu_n\) under a moment condition and shows that the plug-in error from using \(\widehat\mu_n\) in the second-stage density estimator is asymptotically negligible in the stated normalization [1705.07578].

Identifiability is more delicate than the latent-Gaussian hierarchy might suggest. In the model
\[
Y_i=X_i+S_i\varepsilon_i,\qquad \varepsilon_i\sim N(0,1),
\]
with i.i.d. latent pairs \((X_i,S_i)\), the mixing distribution \(\Pi\) is generally not identifiable when the latent shift is unbounded. The same work proves identifiability under bounded-shift conditions in the model class it studies, by comparing characteristic functions and using uniqueness of the Laplace transform. It also shows that generalized maximum likelihood can be inconsistent even when the model is identifiable, and that independence of shift and scale does not remove this inconsistency. By contrast, if there is more than one observation per latent realization, the resulting model becomes a mixture of a bounded full-rank exponential family, and the MLE/GMLE exists, is unique, and is consistent [2408.09195].

## 5. Shape theory, geometry, and relation to adjacent models

MVMN densities inherit important shape properties from the mixing distribution. In the univariate model
\[
Y=\mu+\beta X+\sqrt X\,Z,
\]
if the mixing density \(g\) is unimodal, then the mixture density \(f\) is unimodal; if \(g\) is log-concave, then \(f\) is log-concave; and if \(g\) is log-convex on \((0,\infty)\), then \(f\) is log-convex on each of \((-\infty,\mu)\) and \((\mu,\infty)\). If \(g\) is decreasing, or if \(\beta=0\), then the only mode is at \(\mu\). In the multivariate case, the corresponding results use the corrected mixing density
\[
g^*(x)=x^{-(p-1)/2}g(x).
\]
If \(g^*\) is unimodal, the density has only one local maximum, lying on the line \(y=\mu+\beta t\); if \(g^*\) is log-concave, then the multivariate density is log-concave [1106.2333].

The multivariate geometry is correspondingly rigid. For each \(t\), the density has ellipsoidal contours on the hyperplane
\[
H_t=\mu+\beta t+\{z\in\mathbb R^p:\ z^\top(AA^\top)^{-1}\beta=0\},
\]
so the density is spindle-like: it is constant on ellipsoids orthogonal to the direction \(\beta\). In the GH case, these general shape results yield an especially clean classification: all GH densities are unimodal; in the univariate case GH is log-concave iff \(\lambda\ge 1\); and in the multivariate case GH is log-concave iff
\[
\lambda\ge \frac{p+1}{2}.
\]
The same paper presents these results as a short proof of unimodality for all generalized hyperbolic densities [1106.2333].

Several neighboring models are related but not identical to MVMN. Mean-mixtures of multivariate normals randomize only the mean,
\[
\mathbf Y=\boldsymbol\xi+\boldsymbol\omega(\boldsymbol\delta U+\mathbf Z),
\]
with fixed conditional covariance, and are explicitly described as not obtained from the MVMN class because they lack the \(\sqrt U\) variance mixing term [2006.10018]. Conversely, finite discrete mixtures of multivariate normals, such as two-regime portfolio models with
\[
r_M=
\begin{cases}
r_{M,1} & \text{w.p. }\rho_1,\\
r_{M,2} & \text{w.p. }\rho_2,
\end{cases}
\qquad
r_{M,i}\sim \mathcal N(\mu^i,\Sigma^i),
\]
remain mixtures of Gaussians but are regime mixtures rather than continuous mean-variance mixtures driven by a single positive scalar [1706.03287].

## 6. Portfolio theory, tail risk, and capital allocation

In portfolio applications, MVMN is usually written as
\[
X=\mu+\gamma Z+\sqrt Z\,A N_d,
\]
with \(Z>0\) independent of \(N_d\). For any portfolio \(\omega\),
\[
\omega^T X\overset{d}{=}\omega^T\mu+(\omega^T\gamma)Z+\sqrt Z\,\|\omega\|_\Sigma\,N(0,1),
\qquad
\|\omega\|_\Sigma=\sqrt{\omega^T\Sigma\omega}.
\]
A central stochastic-dominance result states that for
\[
Q=a_1+b_1Z+c_1\sqrt Z\,N,\qquad R=a_2+b_2Z+c_2\sqrt Z\,N,
\]
the conditions
\[
a_1+b_1E Z\ge a_2+b_2E Z,\qquad c_1\le c_2
\]
are sufficient for second-order stochastic dominance under the stated integrability assumptions. This leads to a closed-form frontier theorem: if returns follow the NMVM class and \(\rho\) is any finite-valued, law-invariant, convex risk measure on \(L^k\), then the frontier portfolio for a target return \(r\) is obtained by solving the classical Markowitz mean-variance problem with adjusted mean
\[
\mu_\theta=\mu+\gamma E Z
\]
and covariance
\[
\Sigma_\theta=A^TA=\Sigma.
\]
The optimizer therefore depends on the shifted mean and on the Gaussian covariance component, not on the full covariance
\[
\operatorname{Cov}(X)=\gamma\gamma^T\operatorname{Var}(Z)+\Sigma E Z.
\]
This applies in particular to \(CVaR_\alpha\) and to law-invariant coherent or spectral risk measures of Kusuoka type [2202.02488].

A related portfolio-optimization result concerns mean-risk-skewness criteria. After transforming portfolio weights by \(x^T=\omega^T A\), the latent representation yields
\[
\omega^T X\overset{d}{=}x^T Y,
\]
and, under the condition
\[
m_3(Z)E[Z]\ge 2\operatorname{Var}(Z)^2,
\]
the mean-risk-skewness problem for any law-invariant coherent risk measure reduces to the quadratic program
\[
\min_{x\in\mathbb R^n}x^T x
\quad\text{s.t.}\quad
x^T m=r,\;\;x^T e_A=1.
\]
The same work derives exact integral formulas and approximate closed-form formulas for portfolio VaR and CVaR under general NMVM returns. The approximations use only the one-dimensional risks of
\[
Y_b=bZ+\sqrt Z\,N,
\qquad
Y_{-b}=-bZ+\sqrt Z\,N,
\]
and are reported to be accurate and computationally efficient in the numerical study [2111.04311].

Recent work on capital allocation extends MVMN methodology beyond CTE. For aggregate loss \(S=X_1+\cdots+X_n\), the \(k\)-th order tail central moment is defined by
\[
TCM_k(S)=E[(S-CTE(S))^k\mid S>s_\alpha],
\]
and the corresponding capital allocation rule is
\[
K_i=\operatorname{Cov}\!\left(X_i,(S-CTE(S))^{k-1}\right),\qquad i=1,\dots,n.
\]
Under the NMVM model, the paper derives explicit formulas for tail moments, tail central moments, and the allocation itself:
\[
K_i
=
a_{1,i}\,TCM_k(S)
+
a_{2,i}\big( TCM_{k-1}(S^*)-TCM_{k-1}(S)\big),
\]
with coefficients determined by \(\boldsymbol\mu\), \(\boldsymbol\gamma\), and \(\Sigma\). In the four-asset GH illustration, BA and CVX have fairly stable allocation proportions across CTE, TV, and \(TCM_3\), whereas AXP receives larger TV and \(TCM_3\) allocations than its CTE allocation, and XOM can have a negative \(TCM_3\) contribution [2601.00568].

Across these applications, MVMN is used not merely as a heavier-tailed replacement for the multivariate normal, but as a latent-variable class in which skewness, tail thickness, and dependence are jointly parameterized while much of the Gaussian conditional structure remains available for optimization, conditioning, and likelihood-based inference [2005.06883].

Source: https://www.emergentmind.com/topics/mean-variance-mixture-of-normals-mvmn