cpt: Centered planar transform

Description Usage Arguments Details Value Author(s) References See Also Examples

Description

Compute the centered planar transform of a (dataset of) compositions and its inverse.

Usage

1
2
3
          cpt( x,... )
          cptInv( z,... )
          

Arguments

x

a composition or a data.matrix of compositions, not necessarily closed

z

the cpt-transform of a composition or a data matrix of cpt-transforms of compositions. It is checked that the z sum up to 0.

...

generic arguments. not used.

Details

The cpt-transform maps a composition in the D-part real-simplex isometrically to a D-1 dimensional euclidian vector space, identified with a plane parallel to the simplex but passing through the origin. However the transformation is not injective and does not even reach the whole plane. Thus resulting covariance matrices are always singular.

The data can then be analysed in this transformed space by all classical multivariate analysis tools not relying on a full rank of the covariance matrix. See ipt and apt for alternatives. The interpretation of the results is relatively easy since the relation of each transformed component to the original parts is preserved.

The centered planar transform is given by

cpt(x)_i := clo(x)_i - 1/D

Value

cpt gives the centered planar transform, cptInv gives closed compositions with the given cpt-transforms.

Author(s)

K.Gerald v.d. Boogaart http://www.stat.boogaart.de

References

van den Boogaart, K.G. and R. Tolosana-Delgado (2008) "compositions": a unified R package to analyze Compositional Data, Computers & Geosciences, 34 (4), pages 320-338, doi:10.1016/j.cageo.2006.11.017.

See Also

clr,apt,ipt

Examples

1
2
3
4
5
6
(tmp <- cpt(c(1,2,3)))
cptInv(tmp)
cptInv(tmp) - clo(c(1,2,3)) # 0
data(Hydrochem)
cdata <- Hydrochem[,6:19]
pairs(cpt(cdata),pch=".") 


Search within the compositions package
Search all R packages, documentation and source code

Questions? Problems? Suggestions? or email at ian@mutexlabs.com.

Please suggest features or report bugs with the GitHub issue tracker.

All documentation is copyright its authors; we didn't write any of that.