Additive logistic transform (unconstrained coordinates to the simplex)

Implements the reparametrisation $\boldsymbol{\psi}:\boldsymbol{\rho}\mapsto\boldsymbol{p}$ used in the article, sending unconstrained coordinates $\boldsymbol{\rho}\in\mathbb{R}^{J-1}$ to cellular proportions $\boldsymbol{p}\in\Delta^{J-1}$. This is the additive logistic transform of Aitchison, i.e. the inverse additive log-ratio map ($\mathrm{alr}^{-1}$), equivalently a softmax with the last category $J$ pinned as reference ($\rho_J\equiv 0$).

Recovers the unconstrained additive log-ratio coordinates $\rho_j=\ln(p_j/p_J)$ for $j=1,\ldots,J-1$, with the last part $p_J$ as reference. This is Aitchison's additive log-ratio ($\mathrm{alr}$) transform, equivalently the multinomial-logit link with reference category $J$ (see compositions::alr() and vignette("softmax-alr-derivatives", package = "DeCovarT")).

Usage

additive_logistic(rho)

additive_log_ratio(p)

Arguments

rho: Numeric vector $\boldsymbol{\rho}\in\mathbb{R}^{J-1}$ of unconstrained additive log-ratio coordinates (reference cell type $J$).
p: Numeric vector $\boldsymbol{p}$ on the open simplex.

Value

Numeric vector $\boldsymbol{p}\in\mathbb{R}^{J}$ on the unit simplex ($\mathbf{1}^{\mathsf{T}}\boldsymbol{p}=1$, $\boldsymbol{p}\ge\mathbf{0}$).

Numeric vector $\boldsymbol{\rho}\in\mathbb{R}^{J-1}$.

Details

With $A=\sum_{k=1}^{J-1}\mathrm{e}^{\rho_k}+1$, $$ p_j=\frac{\mathrm{e}^{\rho_j}}{A}\quad(j<J),\qquad p_J=\frac{1}{A}. $$ Equivalently, $\boldsymbol{p}=\boldsymbol{\psi}(\boldsymbol{\rho})$ with $\boldsymbol{\psi}(\boldsymbol{\rho})\propto (\mathrm{e}^{\rho_1},\ldots,\mathrm{e}^{\rho_{J-1}},1)^{\mathsf{T}}$. Jacobians and Hessians of both maps are derived in the package vignette vignette("softmax-alr-derivatives", package = "DeCovarT"). See also compositions::alrInv().

Examples

rho <- c(0.2, -0.5)
p <- additive_logistic(rho)
sum(p)
#> [1] 1
additive_log_ratio(p)
#> [1]  0.2 -0.5