Sensitivity index

From HandWiki
Figure 1: Bayes-optimal classification error probability [math]\displaystyle{ e_b }[/math] and Bayes discriminability index [math]\displaystyle{ d'_b }[/math] between two univariate histograms computed from their overlap area. Figure 2: Same computed from the overlap volume of two bivariate histograms. Figure 3: discriminability indices of two univariate normal distributions with unequal variances. The classification boundary is in black. Figure 4: discriminability indices of two bivariate normal distributions with unequal covariance matrices (ellipses are 1 sd error ellipses). Color-bar shows the relative contribution to the discriminability by each dimension. These are computed by numerical methods[1].

The sensitivity index or discriminability index or detectability index is a dimensionless statistic used in signal detection theory. A higher index indicates that the signal can be more readily detected.

Definition

The discriminability index is the separation between the means of two distributions (typically the signal and the noise distributions), in units of the standard deviation.

Equal variances/covariances

For two univariate distributions [math]\displaystyle{ a }[/math] and [math]\displaystyle{ b }[/math] with the same standard deviation, it is denoted by [math]\displaystyle{ d' }[/math] ('dee-prime'):

[math]\displaystyle{ d' = \frac{\left\vert \mu_a - \mu_b \right\vert}{\sigma} }[/math].

In higher dimensions, i.e. with two multivariate distributions with the same variance-covariance matrix [math]\displaystyle{ \mathbf{\Sigma} }[/math], (whose symmetric square-root, the standard deviation matrix, is [math]\displaystyle{ \mathbf{S} }[/math]), this generalizes to the Mahalanobis distance between the two distributions:

[math]\displaystyle{ d'=\sqrt{(\boldsymbol{\mu}_a-\boldsymbol{\mu}_b)'\mathbf{\Sigma}^{-1}(\boldsymbol{\mu}_a-\boldsymbol{\mu}_b)} = \lVert \mathbf{S}^{-1}(\boldsymbol{\mu}_a-\boldsymbol{\mu}_b) \rVert = \lVert \boldsymbol{\mu}_a-\boldsymbol{\mu}_b \rVert /\sigma_{\boldsymbol{\mu}} }[/math],

where [math]\displaystyle{ \sigma_{\boldsymbol{\mu}} = 1/ \lVert\mathbf{S}^{-1}\boldsymbol{\mu}\rVert }[/math] is the 1d slice of the sd along the unit vector [math]\displaystyle{ \boldsymbol{\mu} }[/math] through the means, i.e. the [math]\displaystyle{ d' }[/math] equals the [math]\displaystyle{ d' }[/math] along the 1d slice through the means.[1]

For two bivariate distributions with equal variance-covariance, this is given by:

[math]\displaystyle{ {d'}^2 =\frac{1}{1-\rho^2} \left({d'}^2_x+{d'}^2_y-2\rho {d'}_x {d'}_y \right) }[/math],

where [math]\displaystyle{ \rho }[/math] is the correlation coefficient, and here [math]\displaystyle{ d'_x=\frac{{\mu_b}_x-{\mu_a}_x}{\sigma_x} }[/math] and [math]\displaystyle{ d'_y=\frac{{\mu_b}_y-{\mu_a}_y}{\sigma_y} }[/math], i.e. including the signs of the mean differences instead of the absolute.[1]

[math]\displaystyle{ d' }[/math] is also estimated as [math]\displaystyle{ Z(\text{hit rate})-Z(\text{false alarm rate}) }[/math].[2]:8

Unequal variances/covariances

When the two distributions have different standard deviations (or in general dimensions, different covariance matrices), there exist several contending indices, all of which reduce to [math]\displaystyle{ d' }[/math] for equal variance/covariance.

Bayes discriminability index

This is the maximum (Bayes-optimal) discriminability index for two distributions, based on the amount of their overlap, i.e. the optimal (Bayes) error of classification [math]\displaystyle{ e_b }[/math] by an ideal observer, or its complement, the optimal accuracy [math]\displaystyle{ a_b }[/math]:

[math]\displaystyle{ d'_b=-2Z\left(\text{Bayes error rate } e_b\right)=2Z\left(\text{best accuracy rate } a_b\right) }[/math],[1]

where [math]\displaystyle{ Z }[/math] is the inverse cumulative distribution function of the standard normal. The Bayes discriminability between univariate or multivariate normal distributions can be numerically computed [1] (Matlab code), and may also be used as an approximation when the distributions are close to normal.

[math]\displaystyle{ d'_b }[/math] is a positive-definite statistical distance measure that is free of assumptions about the distributions, like the Kullback-Leibler divergence [math]\displaystyle{ D_\text{KL} }[/math]. [math]\displaystyle{ D_\text{KL}(a,b) }[/math] is asymmetric, whereas [math]\displaystyle{ d'_b(a,b) }[/math] is symmetric for the two distributions. However, [math]\displaystyle{ d'_b }[/math] does not satisfy the triangle inequality, so it is not a full metric. [1]

In particular, for a yes/no task between two univariate normal distributions with means [math]\displaystyle{ \mu_a,\mu_b }[/math] and variances [math]\displaystyle{ v_a\gt v_b }[/math], the Bayes-optimal classification accuracies are:[1]

[math]\displaystyle{ p(A|a)=p({\chi'}^2_{1,v_a \lambda} \gt v_b c), \; \; p(B|b)=p({\chi'}^2_{1,v_b \lambda} \lt v_a c) }[/math],

where [math]\displaystyle{ \chi'^2 }[/math] denotes the non-central chi-squared distribution, [math]\displaystyle{ \lambda=\left(\frac{\mu_a-\mu_b}{v_a-v_b}\right)^2 }[/math], and [math]\displaystyle{ c=\lambda+\frac{\ln v_a -\ln v_b}{v_a-v_b} }[/math]. The Bayes discriminability [math]\displaystyle{ d'_b=2Z\left(\frac{p\left(A|a\right)+p\left(B|b\right)}{2} \right). }[/math]

[math]\displaystyle{ d'_b }[/math] can also be computed from the ROC curve of a yes/no task between two univariate normal distributions with a single shifting criterion. It can also be computed from the ROC curve of any two distributions (in any number of variables) with a shifting likelihood-ratio, by locating the point on the ROC curve that is farthest from the diagonal. [1]

For a two-interval task between these distributions, the optimal accuracy is [math]\displaystyle{ a_b=p \left( \tilde{\chi}^2_{\boldsymbol{w}, \boldsymbol{k}, \boldsymbol{\lambda},0,0}\gt 0 \right) }[/math] ([math]\displaystyle{ \tilde{\chi}^2 }[/math] denotes the generalized chi-squared distribution), where [math]\displaystyle{ \boldsymbol{w}=\begin{bmatrix} \sigma_s^2 & -\sigma_n^2 \end{bmatrix}, \; \boldsymbol{k}=\begin{bmatrix} 1 & 1 \end{bmatrix}, \; \boldsymbol{\lambda}=\frac{\mu_s-\mu_n}{\sigma_s^2-\sigma_n^2} \begin{bmatrix} \sigma_s^2 & \sigma_n^2 \end{bmatrix} }[/math].[1] The Bayes discriminability [math]\displaystyle{ d'_b=2Z\left(a_b\right) }[/math].

RMS sd discriminability index

A common approximate (i.e. sub-optimal) discriminability index that has a closed-form is to take the average of the variances, i.e. the rms of the two standard deviations: [math]\displaystyle{ d'_a=\left\vert \mu_a -\mu_b \right\vert/\sigma_\text{rms} }[/math] [3] (also denoted by [math]\displaystyle{ d_a }[/math]). It is [math]\displaystyle{ \sqrt{2} }[/math] times the [math]\displaystyle{ z }[/math]-score of the area under the receiver operating characteristic curve (AUC) of a single-criterion observer. This index is extended to general dimensions as the Mahalanobis distance using the pooled covariance, i.e. with [math]\displaystyle{ \mathbf{S}_\text{rms}=\left[\left(\mathbf{\Sigma}_a+\mathbf{\Sigma}_b\right)/2 \right]^\frac{1}{2} }[/math] as the common sd matrix.[1]

Average sd discriminability index

Another index is [math]\displaystyle{ d'_e=\left\vert \mu_a -\mu_b \right\vert/\sigma_\text{avg} }[/math], extended to general dimensions using [math]\displaystyle{ \mathbf{S}_\text{avg}=\left(\mathbf{S}_a+\mathbf{S}_b\right)/2 }[/math] as the common sd matrix.[1]

Comparison of the indices

It has been shown that for two univariate normal distributions, [math]\displaystyle{ d'_a \leq d'_e \leq d'_b }[/math], and for multivariate normal distributions, [math]\displaystyle{ d'_a \leq d'_e }[/math] still.[1]

Thus, [math]\displaystyle{ d'_a }[/math] and [math]\displaystyle{ d'_e }[/math] underestimate the maximum discriminability [math]\displaystyle{ d'_b }[/math] of univariate normal distributions. [math]\displaystyle{ d'_a }[/math] can underestimate [math]\displaystyle{ d'_b }[/math] by a maximum of approximately 30%. At the limit of high discriminability for univariate normal distributions, [math]\displaystyle{ d'_e }[/math] converges to [math]\displaystyle{ d'_b }[/math]. These results often hold true in higher dimensions, but not always.[1] Simpson and Fitter [3] promoted [math]\displaystyle{ d'_a }[/math] as the best index, particularly for two-interval tasks, but Das and Geisler [1] have shown that [math]\displaystyle{ d'_b }[/math] is the optimal discriminability in all cases, and [math]\displaystyle{ d'_e }[/math] is often a better closed-form approximation than [math]\displaystyle{ d'_a }[/math], even for two-interval tasks.

The approximate index [math]\displaystyle{ d'_{gm} }[/math], which uses the geometric mean of the sd's, is less than [math]\displaystyle{ d'_b }[/math] at small discriminability, but greater at large discriminability.[1]

Contribution to discriminability by each dimension

In general, the contribution to the total discriminability by each dimension or feature may be measured using the amount by which the discriminability drops when that dimension is removed. If the total Bayes discriminability is [math]\displaystyle{ d' }[/math] and the Bayes discriminability with dimension [math]\displaystyle{ i }[/math] removed is [math]\displaystyle{ d'_{-i} }[/math], we can define the contribution of dimension [math]\displaystyle{ i }[/math] as [math]\displaystyle{ \sqrt{d'^2-{d'_{-i}}^2} }[/math]. This is the same as the individual discriminability of dimension [math]\displaystyle{ i }[/math] when the covariance matrices are equal and diagonal, but in the other cases, this measure more accurately reflects the contribution of a dimension than its individual discriminability.[1]

See also

References

  1. 1.00 1.01 1.02 1.03 1.04 1.05 1.06 1.07 1.08 1.09 1.10 1.11 1.12 1.13 1.14 1.15 Das, Abhranil; Wilson S Geisler (2020). "Methods to integrate multinormals and compute classification measures". arXiv:2012.14331 [stat.ML].
  2. MacMillan, N.; Creelman, C. (2005). Detection Theory: A User's Guide. Lawrence Erlbaum Associates. ISBN 9781410611147. https://books.google.com/books?id=hDX65v9bReYC. 
  3. 3.0 3.1 Simpson, A. J.; Fitter, M. J. (1973). "What is the best index of detectability?". Psychological Bulletin 80 (6): 481–488. doi:10.1037/h0035203. 

External links