medcouple: A robust measure of skewness for univariate data

In phylter: Detect and Remove Outliers in Phylogenomics Datasets

A robust measure of skewness for univariate data

Description

Computes the medcouple, a robust measure of skewness for univariate data. For multivariate data the medcouple is calculated on each column of the data matrix.

Usage

medcouple(x, do.reflect = NULL)

Arguments

x	An n by p data matrix.
do.reflect	Logical indicating whether the medcouple should also be computed on the reflected sample -x, with final result (mc(x)-mc(-x))/2. Defaults TRUE when n<=100 and to FALSE otherwise.

Details

The medcouple is a robust measure of skewness yielding values between -1 and 1. For left- and right-skewed data the medcouple is negative and positive respectively.
The medcouple is defined as the median of the kernel function h(x_i,x_j) = \frac{(x_j - med(x)) - (med(x)-x_i)}{x_j-x_i} evaluated over all couples (x_i,x_j) where x_i is smaller than the median of x and x_j larger than the median of x. When there are multiple observations tied to the median, the kernel is defined separately as the denominator is not defined for these observations. Let m_1 < ... < m_k denote the indices of the observations which are tied to the median. Then h(x_{m_i},x_{m_j}) is defined to equal -1 if i + j - 1 < k, 0 when i + j - 1 = k and +1 if i + j - 1 > k. To compute the medcouple an algorithm with time complexity O(n log(n)) is applied. For details, see https://en.wikipedia.org/wiki/Medcouple. For numerical accuracy it is advised, for small data sets, to compute the medcouple on both x and -x. The final value of the medcouple may then be obtained as a linear combination of both calculations. This procedure is warranted by the properties of the medcouple. Indeed the medcouple of the distribution X equals minus the medcouple of the reflected distribution -X. Moreover the medcouple is location and scale invariant. Note that missing values are not allowed.

Value

mc A p-vector containing the medcouple of each column of the data matrix x.

Note

This function is extracted from the package mrfDepth - 07/2023

Author(s)

P. Segaert with original code from M. Maechler and G. Brys.

References

Brys G., Hubert M., Struyf A. (2004). A robust measure of skewness. Journal of Computational and Graphical Statistics, 13, 996–1017.

Examples

# Calculate the medcouple of univariate data sets.
# For 2000 normally distributed values
# the medcouple value is close to 0 because 
# data are not skewed
x<-rnorm(2000)
medcouple(x) 
# For 2000 values following a lognormal
# distribution (mean 0,sd 1), medcouple is close to 1
# because values are right-skewed
y<-rnorm(2000)
medcouple(y) 
# Use the option do.reflect to increase expected accuracy. 
medcouple(y, do.reflect = TRUE)

phylter documentation built on Aug. 24, 2023, 9:10 a.m.