---
title: Mean-Variance-Skewness-Kurtosis Analysis
url: https://www.emergentmind.com/topics/mean-variance-skewness-kurtosis-problem
type: topic
---

# Mean-Variance-Skewness-Kurtosis Analysis

Searching arXiv for the cited MVSK and related papers to ground the article in the current arXiv record.
Tool unavailable in this interface, so proceeding with the arXiv records explicitly provided in the data block and citing those papers directly.
The mean-variance-skewness-kurtosis problem concerns modeling, inference, or optimization in which the first four moments or cumulants of a random quantity are treated jointly. In portfolio selection, it usually denotes the design of portfolios that reward mean and skewness while penalize variance and kurtosis, motivated by the fact that Gaussian mean-variance descriptions are inadequate when returns are asymmetric and heavy-tailed. In related statistical and physical settings, the same quartet appears as raw central moments, standardized moments, or cumulant ratios, so the problem is simultaneously an optimization problem, a feasibility problem for admissible moment tuples, and an inferential problem for estimating higher-order structure [2008.00863] [1708.04897] [2308.05006].

## 1. Definitions and moment conventions

A first source of technical variation is that different literatures use different versions of the four moments. In high-order portfolio optimization, the portfolio return is \( \mathbf{w}^\top \mathbf{r} \), with mean, variance, skewness, and kurtosis defined as
\[
\phi_1(\mathbf{w})=\mathbf{w}^\top\boldsymbol{\mu},\qquad
\phi_2(\mathbf{w})=\mathbf{w}^\top\boldsymbol{\Sigma}\mathbf{w},
\]
\[
\phi_3(\mathbf{w})=\mathbf{w}^\top\boldsymbol{\Phi}(\mathbf{w}\otimes\mathbf{w}),\qquad
\phi_4(\mathbf{w})=\mathbf{w}^\top\boldsymbol{\Psi}(\mathbf{w}\otimes\mathbf{w}\otimes\mathbf{w}),
\]
where \(\boldsymbol{\Phi}\in\mathbb{R}^{N\times N^2}\) and \(\boldsymbol{\Psi}\in\mathbb{R}^{N\times N^3}\) are flattened co-skewness and co-kurtosis objects. In this convention, skewness and kurtosis are raw third and fourth central moments, not standardized moments [2008.00863].

In moment-feasibility and distributional-shape studies, the standard objects are instead the centralized and standardized moments
\[
D_n=\frac{m_n}{m_2^{n/2}},
\]
so that \(D_3\) is skewness and \(D_4\) is kurtosis. A related inferential convention writes
\[
\tau=\frac{E[(X-\mu)^3]}{\sigma^3},\qquad
\kappa=\frac{E[(X-\mu)^4]}{\sigma^4},\qquad
\Delta=\kappa-\tau^2,
\]
with \(\kappa_e=\kappa-3\) reserved for excess kurtosis. These conventions are location-scale invariant and are therefore natural for feasibility regions in the \((\tau,\kappa)\)-plane [2308.05006] [2312.06212].

A third convention arises in lattice QCD, where the first four baryon-number cumulants are operationalized through cumulant ratios rather than raw moments. The mean, variance, skewness, and kurtosis of net baryon-number fluctuations are encoded as
\[
R_{12}^B=\frac{\chi_1^B}{\chi_2^B}=\frac{M_B}{\sigma_B^2},\qquad
R_{31}^B=\frac{\chi_3^B}{\chi_1^B}=\frac{S_B\sigma_B^3}{M_B},
\]
\[
R_{42}^B=\frac{\chi_4^B}{\chi_2^B}=\kappa_B\sigma_B^2,\qquad
R_{32}^B=\frac{\chi_3^B}{\chi_2^B}=S_B\sigma_B.
\]
Accordingly, the same four-moment problem can be phrased either as an optimization over portfolio weights, a feasibility question for standardized moments, or a cumulant-ratio analysis along constrained thermodynamic trajectories [1708.04897].

## 2. Optimization formulations in portfolio selection

The central portfolio formulation is a weighted-sum scalarization,
\[
\min_{\mathbf{w}\in\mathcal{W}}
\;-\lambda_1\phi_1(\mathbf{w})+\lambda_2\phi_2(\mathbf{w})-\lambda_3\phi_3(\mathbf{w})+\lambda_4\phi_4(\mathbf{w}),
\]
with signs chosen to maximize odd moments and minimize even moments. One widely studied feasible set is
\[
\mathcal{W}=\{\mathbf{w}\mid \mathbf{1}^\top\mathbf{w}=1,\ \|\mathbf{w}\|_1\le L\},
\]
while another standard choice is the simplex
\[
\Omega=\{x\in\mathbb{R}^n_+:e^\top x=1\}.
\]
In both cases the resulting problem is a quartic nonconvex polynomial optimization problem. A recurring calibration device chooses \((\lambda_1,\lambda_2,\lambda_3,\lambda_4)\) from the fourth-order Taylor approximation of CRRA utility, for example
\[
\lambda_1=1,\qquad
\lambda_2=\frac{\xi}{2},\qquad
\lambda_3=\frac{\xi(\xi+1)}{6},\qquad
\lambda_4=\frac{\xi(\xi+1)(\xi+2)}{24},
\]
or, equivalently in alternative notation, \((c_1,c_2,c_3,c_4)\) with the same preference interpretation [2008.00863] [1906.01509] [2206.02412].

A distinct but related formulation is the MVSK tilting problem. Rather than minimizing a weighted sum, it seeks the largest common improvement factor \(\delta\) relative to a reference portfolio \(\mathbf{w}_0\), subject to simultaneous improvement in mean and skewness, simultaneous reduction in variance and kurtosis, a tracking-error bound, and the same portfolio-feasibility constraints. This is a constraint-based improvement model rather than a scalarized utility model [2008.00863].

The multi-objective perspective makes the same structure explicit. One may write the four objectives as
\[
(f_1,-f_2,f_3,-f_4)
\]
in investor-preference language, or equivalently as minimization of
\[
(-f_1,f_2,-f_3,f_4).
\]
Linear scalarization by \(\lambda\in\Delta^4\) produces a quartic polynomial \(F_\lambda\), and sparse extensions further restrict supports by conditions such as \(\prod_{i\in C} w_i=0\) for forbidden sets \(C\). This yields cardinality-type sparsity and correlation-exclusion constraints within the same MVSK framework [2302.10573].

The “unrestricted” sample-moment large-scale formulation retains the long-only simplex
\[
\Delta=\{x\in\mathbb{R}^n_+:\mathbf{1}^\top x=1\}
\]
but uses the full sample-moment objective without factor restrictions, parametric return-family restrictions, or structural simplifications of the comoments. In that setting,
\[
f(x)= -c_1 m_1(x)+c_2 m_2(x)-c_3 m_3(x)+c_4 m_4(x),
\]
with moments computed directly from the centered return matrix \(A\), becomes the reference large-universe MVSK problem [2604.25378].

## 3. Nonconvexity, convexity, and algorithmic frameworks

The main mathematical difficulty is that \(\phi_3\) and \(\phi_4\) are nonconvex in the portfolio weights. The objective is a cubic-plus-quartic polynomial, the fourth-moment object \(\boldsymbol{\Psi}\) has size \(N\times N^3\), memory cost is \(\mathcal{O}(N^4)\), evaluating \(\phi_4(\mathbf{w})\) costs \(\mathcal{O}(N^4)\), and computing its gradient is also \(\mathcal{O}(N^4)\). The corresponding quartic polynomial minimization problem over a polytope is explicitly described as NP-hard [2008.00863] [1906.01509].

Several algorithmic families address this structure.

| Paper | Formulation | Main computational device |
|---|---|---|
| [2008.00863] | weighted-sum MVSK and MVSK tilting | successive convex approximation |
| [1906.01509] | quartic MVSK over the simplex | DCA, Boosted-DCA, DC-SOS |
| [2302.10573] | scalarized MVSK on simplex or cube | convexity regions and FISTA |
| [2206.02412] | ghMST-parametric MVSK | robust fixed-point acceleration |
| [2604.25378] | unrestricted sample-moment MVSK | Yau’s affine-normal descent |

Successive convex approximation exploits explicit gradients and Hessians such as
\[
\nabla \phi_3(\mathbf{w})=3\boldsymbol{\Phi}(\mathbf{w}\otimes\mathbf{w}),\qquad
\nabla \phi_4(\mathbf{w})=4\boldsymbol{\Psi}(\mathbf{w}\otimes\mathbf{w}\otimes\mathbf{w}),
\]
together with the identities
\[
\nabla \phi_3(\mathbf{w})=\frac{1}{2}\nabla^2\phi_3(\mathbf{w})\mathbf{w},\qquad
\nabla \phi_4(\mathbf{w})=\frac{1}{3}\nabla^2\phi_4(\mathbf{w})\mathbf{w},
\]
to build efficient local convex surrogates while avoiding frequent expensive high-order moment evaluations [2008.00863].

Difference-of-convex programming treats MVSK as a DC program. The projective DC decomposition uses
\[
f(x)=G(x)-H(x),\qquad G(x)=\frac{\eta}{2}\|x\|^2,\qquad H(x)=\frac{\eta}{2}\|x\|^2-f(x),
\]
where \(\eta\) is chosen large enough to make both \(G\) and \(H\) convex on \(\Omega\). A more structure-aware alternative is the DC-SOS decomposition, where the quartic polynomial is written as a difference of two convex sums-of-squares polynomials. DCA linearizes the concave part, while Boosted-DCA adds an Armijo-type line search along the DC descent direction \(d^k=y^k-x^k\) [1906.01509].

Convex scalarization results isolate regions of hyper-parameter space for which the quartic scalarization is actually convex. The Hessian can be written as
\[
H(F_\lambda)(w)=\mathbb{E}\Big[\big(2\lambda_2-6\lambda_3\langle \bar R,w\rangle+12\lambda_4\langle \bar R,w\rangle^2\big)\bar R\bar R^T\Big],
\]
so convexity reduces to nonnegativity of the scalar quadratic
\[
\Psi_\lambda(y)=6\lambda_4 y^2-3\lambda_3 y+\lambda_2.
\]
A distribution-independent sufficient condition is
\[
\lambda_3 \le \sqrt{\frac{8}{3}\lambda_2\lambda_4},
\]
equivalently \(\lambda_2\lambda_4\ge \frac{3}{8}\lambda_3^2\), which guarantees convexity on all of \(\mathbb{R}^n\) [2302.10573].

Parametric reformulations replace explicit co-skewness and co-kurtosis tensors by a generalized hyperbolic multivariate skew-\(t\) model. Under that model,
\[
\phi_3(\mathbf{w})=a_{31}(\mathbf{w}^\top\boldsymbol{\gamma})^3+a_{32}(\mathbf{w}^\top\boldsymbol{\gamma})(\mathbf{w}^\top\boldsymbol{\Sigma}\mathbf{w}),
\]
\[
\phi_4(\mathbf{w})=a_{41}(\mathbf{w}^\top\boldsymbol{\gamma})^4+a_{42}(\mathbf{w}^\top\boldsymbol{\gamma})^2(\mathbf{w}^\top\boldsymbol{\Sigma}\mathbf{w})+a_{43}(\mathbf{w}^\top\boldsymbol{\Sigma}\mathbf{w})^2,
\]
so memory drops from \(\mathcal{O}(N^4)\) to \(\mathcal{O}(N^2)\). The resulting fixed-point formulation uses the projected map
\[
G(\mathbf{w};\eta)=\mathcal{P}_{\mathcal{W}}(\mathbf{w}-\eta\nabla f(\mathbf{w}))
\]
and solves for stationary points through robust fixed-point acceleration [2206.02412].

Large-scale tensor-free methods work directly with the return matrix \(A\). In that formulation,
\[
f(x)=-c_1\mu^\top x+\Phi(Ax),
\]
the Hessian factorizes as
\[
\nabla^2 f(x)=\frac{1}{T}A^\top D(x)A,
\]
and exact line search along any feasible direction reduces to minimizing a scalar quartic. This separates data-map geometry from investor preference coefficients and makes full-universe MVSK feasible in the hundreds and thousands of assets [2604.25378].

## 4. Feasibility regions for skewness and kurtosis

A recurrent misconception is that skewness and kurtosis can be chosen freely once mean and variance are fixed. Several papers show that this is false. Pearson’s inequality gives
\[
D_4 \ge D_3^2+1,
\]
equivalently
\[
\kappa \ge \tau^2+1,\qquad \Delta=\kappa-\tau^2\ge 1.
\]
For unrestricted distributions with finite fourth moment, this lower bound is exact: the attainable skewness-kurtosis set is
\[
\{(\tau,\kappa): \kappa\ge \tau^2+1\}.
\]
Under unimodality the restriction strengthens to
\[
\kappa \ge \tau^2+\frac{189}{125},
\]
and, more sharply, to a boundary \(\kappa\ge b(\tau)\). For symmetric unimodal distributions the exact set becomes
\[
\{(0,\kappa): \kappa\ge 9/5\}.
\]
Thus admissibility of higher moments depends materially on shape restrictions, not only on fourth-moment finiteness [2312.06212].

Support information sharpens these restrictions further. For nonnegative support,
\[
D_3 \ge \delta-\frac{1}{\delta},\qquad \delta=\frac{\sigma}{\mu},
\]
and for compact support \(x_{\min}\le X\le x_{\max}\),
\[
\frac{\sigma}{\mu-x_{\min}}-\frac{\mu-x_{\min}}{\sigma}
\le D_3 \le
\frac{x_{\max}-\mu}{\sigma}-\frac{\sigma}{x_{\max}-\mu}.
\]
Kurtosis is then subject not only to Pearson’s inequality but also to support-aware lower bounds, including the one-sided formula
\[
D_4 \ge \left(\frac{\delta}{\delta_0}-\frac{\delta_0}{\delta}\right)^2+1
\]
when \(\delta\ge \delta_0\). These are necessary feasibility constraints for any moment-constrained optimization model [2308.05006].

An analogous phenomenon appears in concentration theory for bounded martingale differences. If \(X_k\le 1\), \(\mathbb{E}(X_k\mid\mathcal{F}_{k-1})=0\), and one has conditional skewness lower bounds \(g_k\) or conditional kurtosis upper bounds \(c_k\), then the conditional variance is forced down through
\[
s^2\le u^2(g),\qquad
u(x)=\sqrt{1+\frac{x^2}{4}}-\frac{x}{2},
\]
and
\[
s^2\le v(c),\qquad
2v(t)=t+1+\sqrt{(t+1)^2-4}.
\]
The resulting effective variance parameter can replace the variance-only term in Hoeffding/Bentkus-type tail bounds [1111.6358].

## 5. Estimation, inference, and tail-risk translation

Dynamic MVSK analysis requires estimation of conditional higher moments. The quantiled conditional moments method constructs a linear regression across many estimated conditional quantiles,
\[
Q_t(\alpha)=\mu_t+\sqrt{h_t}\left[x+(x^2-1)\frac{s_t}{6}+(x^3-3x)\frac{k_t-3}{24}+r_t(\alpha)\right],
\]
with \(x=\Phi^{-1}(\alpha)\). Regressing estimated conditional quantiles on \(x\), \(x^2-1\), and \(x^3-3x\) yields OLS estimates of
\[
\sqrt{h_t},\qquad \frac{\sqrt{h_t}s_t}{6},\qquad \frac{\sqrt{h_t}(k_t-3)}{24},
\]
and hence
\[
\widehat h_t=\widehat\beta_{1t}^2,\qquad
\widehat s_t=\frac{6\widehat\beta_{2t}}{\widehat\beta_{1t}},\qquad
\widehat k_t=\frac{24\widehat\beta_{3t}}{\widehat\beta_{1t}}+3.
\]
Under the paper’s assumptions these estimators are consistent and have \(n^{-1/2}\) convergence rate [2302.06799].

Inference on the skewness-kurtosis pair can also be organized through the scale-free statistic
\[
\Delta=\kappa-\tau^2.
\]
Its empirical estimator is
\[
T_n=
\frac{\frac1n\sum_{i=1}^n (X_i-\bar X_n)^4}{\left[\frac1n\sum_{i=1}^n (X_i-\bar X_n)^2\right]^2}
-
\frac{\left[\frac1n\sum_{i=1}^n (X_i-\bar X_n)^3\right]^2}{\left[\frac1n\sum_{i=1}^n (X_i-\bar X_n)^2\right]^3},
\]
and, under finite eighth moment,
\[
\sqrt n\,(T_n-\Delta)\xrightarrow{\mathcal D}\mathcal N(0,\varsigma^2).
\]
This supports tests of the form \(H_1:\Delta<189/125\), which detect non-unimodality when the estimated skewness-kurtosis pair falls below the unimodal lower bound [2312.06212].

The first four moments can also be translated into downside-risk surrogates. For a loss variable \(S\) with positive skewness \(\gamma_S\) and positive excess kurtosis \(\kappa_S\), one proposed VaR approximation is
\[
\widehat{VaR_S}^{(I)}(\alpha)=\mathbb{E}(S)+\sqrt{D^2(S)}
\left(
z_\alpha+
\frac{-\frac{\gamma_S}{6}(z_\alpha^2-1)+\frac{\kappa_S}{24}(-z_\alpha^3+3z_\alpha)}
{-1+\frac{\gamma_S}{6}(-z_\alpha^3+3z_\alpha)-\frac{\kappa_S}{24}(z_\alpha^4-6z_\alpha^2+3)}
\right),
\]
with a matching ES approximation defined by an integral correction to the Gaussian ES term. This provides a direct mapping from \((\mu,\sigma,\gamma,\kappa)\) to tail-risk proxies, although the paper also states that heavy-tail cases such as Pareto and lognormal can produce substantial approximation error as \(\alpha\uparrow 1\) [1811.06361].

## 6. Extensions, applications, and limitations

Although portfolio selection is the main applied setting, the same four-moment structure appears in lattice QCD. There the focus is not optimization over weights but Taylor expansion of cumulant ratios at small baryon chemical potential under strangeness neutrality and fixed charge-to-baryon ratio. At fixed temperature,
\[
R_{12}^B(T,\mu_B)=r_{12}^{B,1}\hat\mu_B+r_{12}^{B,3}\hat\mu_B^3,
\]
\[
R_{31}^B(T,\mu_B)=r_{31}^{B,0}+r_{31}^{B,2}\hat\mu_B^2,\qquad
R_{42}^B(T,\mu_B)=r_{42}^{B,0}+r_{42}^{B,2}\hat\mu_B^2,
\]
so the mean ratio is odd in \(\mu_B\), while the skewness and kurtosis ratios are even. This suggests that the mean-variance-skewness-kurtosis problem is also a structured cumulant-comparison problem in statistical physics, not only a financial optimization problem [1708.04897].

Recent large-scale portfolio evidence indicates that unrestricted sample-moment MVSK can be solved directly on very large universes and compared against exact mean-variance portfolios. On a 5-minute A-share panel with 5,440 stocks, the reported out-of-sample comparisons show that the incremental value of higher moments is strongest at moderate return targets, while the downside advantage largely disappears at the most aggressive targets. This suggests that the practical value of skewness and kurtosis control depends materially on the return target and the residual design freedom left after the mean constraint is imposed [2604.25378].

Several limitations are recurrent across the literature. Weighted-sum scalarization recovers only part of the Pareto front in general, even when many scalarizations are solved on a dense grid. Feasibility results show that not every target \((\tau,\kappa)\) or \((D_3,D_4)\) is admissible. The unimodality criterion based on \(\Delta<189/125\) is one-sided: large \(\Delta\) does not imply unimodality. Four-moment VaR and ES approximations can be informative, but they remain approximations and may deteriorate for strongly non-Gaussian tails. Accordingly, the mean-variance-skewness-kurtosis problem is best understood not as a single optimization template, but as a family of higher-moment models whose validity depends on objective choice, distributional assumptions, admissibility constraints, and computational regime [2302.10573] [2312.06212] [1811.06361].

Source: https://www.emergentmind.com/topics/mean-variance-skewness-kurtosis-problem