Skip to contents

Centre of a compositional data set

Source: R/aitcen.R
Compcen.Rd

Formula: xi-hat = C(g_1, ..., g_D), g_i the geometric mean of the ith component

Usage

Compcen(X, total = 1)

Arguments

X

See Usage.

total

See Usage.

Value

List with center, geometric_mean, n, D.

References

Aitchison, A Concise Guide to Compositional Data Analysis, Chapter 2. Verified against the text: the centre estimate is xi-hat = C(g_1, ..., g_D) with g_i the geometric mean of the ith component, and totvar(x) = trace(Gamma) = (1/D) sum_{i<j} var{log(x_i/x_j)}.

Examples

V <- c(1, 2, 3, 4, 5, 6, 7, 8)
Compcen(V)
#> $center
#> [1] 1
#> 
#> $geometric_mean
#> [1] 3.764351
#> 
#> $n
#> [1] 8
#> 
#> $D
#> [1] 1
#> 
#> $method
#> [1] "Compositional centre (closed geometric mean)"
#> 
#> attr(,"class")
#> [1] "morie_rich_result" "list"