cpt | R Documentation |
Compute the centered planar transform of a (dataset of) compositions and its inverse.
cpt( x,... )
cptInv( z,...,orig=gsi.orig(z) )
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. |
orig |
a compositional object which should be mimicked by the inverse transformation. It is especially used to reconstruct the names of the parts. |
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 - \frac1D
cpt
gives the centered planar transform,
cptInv
gives closed compositions with the given cpt-transforms.
K.Gerald v.d. Boogaart http://www.stat.boogaart.de
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, \Sexpr[results=rd]{tools:::Rd_expr_doi("10.1016/j.cageo.2006.11.017")}.
clr
,apt
,ipt
(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=".")
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.