[Paper Review] The Shape of the Noncentral Chi-square Density

Yaming Yu|arXiv (Cornell University)|Jun 26, 2011

Matrix Theory and Algorithms16 references3 citations

TL;DR

This paper provides a rigorous mathematical characterization of the shape of the noncentral chi-square density, showing that for degrees of freedom $\nu < 2$, the density transitions from decreasing to bi-modal as the noncentrality parameter $\lambda$ increases beyond a critical threshold $\lambda_\nu$. The critical value $\lambda_\nu$ is defined by a transcendental equation involving a ratio of modified Bessel functions, and the paper derives precise bounds on the location of the interior mode when bi-modality occurs.

ABSTRACT

A noncentral chi-square density is log-concave if the degree of freedom is nu&gt;=2. We complement this known result by showing that, for each 00 such that the chi-square with nu degrees of freedom and noncentrality parameter lambda has a decreasing density if lambda &lt;= lambda_nu, and is bi-modal otherwise. The critical lambda_nu is characterized by an equation involving a ratio of modified Bessel functions. When an interior mode exists we derive precise bounds on its location.

Motivation & Objective

• To formally establish the shape behavior of the noncentral chi-square density, particularly the emergence of bi-modality for $\nu < 2$ and large $\lambda$, which had been previously considered folklore or supported only numerically.
• To derive a precise characterization of the critical noncentrality parameter $\lambda_\nu$ that separates decreasing and bi-modal regimes for $0 < \nu < 2$.
• To provide analytical bounds on the location of the interior mode when the density is bi-modal, using properties of modified Bessel functions and their ratios.
• To establish that log-concavity holds for $\nu \geq 2$, extending known results and clarifying the boundary case at $\nu = 2$.

Proposed method

• The analysis relies on the logarithmic derivative of the density, expressed through the ratio $r_\nu(x) = I_\nu(x)/I_{\nu-1}(x)$ of modified Bessel functions of the first kind.
• A key function $g_\nu(\lambda)$ is defined involving $r_{\nu/2}(\sqrt{\lambda(\lambda + \nu - 4)})$ and a rational term, whose sign determines the shape behavior.
• The existence and uniqueness of the critical $\lambda_\nu$ is proven via intermediate value theorem and monotonicity analysis of $g_\nu(\lambda)$, which crosses zero exactly once.
• The location of the interior mode is bounded using the second derivative of the log-density and properties of convexity/concavity, with critical points analyzed via implicit differentiation.
• The method leverages known asymptotic and differential properties of modified Bessel functions, particularly the differential equation satisfied by $r_\nu(x)$.
• A sign analysis technique is applied to functions like $h(\lambda)$, showing that sign changes are impossible under certain conditions, thereby proving inequalities on the log-derivative.

Experimental results

Research questions

• RQ1 For $0 < \nu < 2$, what is the threshold value $\lambda_\nu$ beyond which the noncentral chi-square density becomes bi-modal?
• RQ2 How does the location of the interior mode in the bi-modal case depend on $\nu$ and $\lambda$?
• RQ3 Under what conditions is the noncentral chi-square density log-concave, and how does this relate to the value of $\nu$?
• RQ4 Can the critical $\lambda_\nu$ be characterized analytically, and is it expressible in terms of special functions?
• RQ5 What are the precise bounds on the position of the interior mode when the density is bi-modal?

Key findings

• For $0 < \nu < 2$, there exists a unique critical noncentrality parameter $\lambda_\nu \in (4 - \nu, \infty)$ such that the density is decreasing if $\lambda \leq \lambda_\nu$ and bi-modal if $\lambda > \lambda_\nu$.
• $\lambda_\nu$ is characterized as the unique solution to the equation $r_{\nu/2}(\sqrt{\lambda(\lambda + \nu - 4)}) = \frac{\lambda - 2}{\sqrt{\lambda(\lambda + \nu - 4)}}$.
• As $\nu \downarrow 0$, $\lambda_\nu \to 4$, and as $\nu \uparrow 2$, $\lambda_\nu \to 2$, showing continuity at the boundary.
• When a bi-modal shape occurs, the interior mode lies strictly to the right of $\lambda + \nu - 4$, and the paper provides bounds on its location using the second derivative of the log-density.
• The mode location $M(\nu, \lambda)$ is strictly increasing in $\lambda$ for fixed $\nu \in (0,2)$, and the function $x r_\nu(x)$ is strictly increasing in $x$ for all $\nu > 0$.
• The paper confirms that the noncentral chi-square density is log-concave for $\nu \geq 2$, and this property holds with equality at $\nu = 2$ for $\lambda > 2$.