bimodality_sarle {microbiome}R Documentation

Sarle's Bimodality Coefficient

Description

Sarle's bimodality coefficient.

Usage

bimodality_sarle(x, bs.iter = 1, type = "Sarle.finite.sample")

Arguments

x

Data vector for which bimodality will be quantified

bs.iter

Bootstrap iterations

type

Score type ('Sarle.finite.sample' or 'Sarle.asymptotic')

Details

The coefficient lies in (0, 1).

The 'Sarle.asymptotic' version is defined as

b=(g^2 + 1) / k

. This is coefficient of bimodality from Ellison AM Am. J. Bot. 1987, for microbiome analysis it has been used for instance in Shade et al. 2014.

The formula for 'Sarle.finite.sample' (SAS 2012):

b=\frac{g^2 + 1}{k + (3(n-1)^2)/((n-2)(n-3))}

where n is sample size and

In both formulas, g is sample skewness and k is the kth standardized moment (also called the sample kurtosis, or excess kurtosis).

Value

Bimodality score

Author(s)

Contact: Leo Lahti microbiome-admin@googlegroups.com

References

See Also

Check the dip.test from the DIP package for a classical test of multimodality.

Examples

## Not run: 
    # b <- bimodality_sarle(rnorm(100), type='Sarle.finite.sample')

## End(Not run)

[Package microbiome version 1.6.0 Index]