clrp {philr} | R Documentation |
Weighted CLR Transform
clrp(y, p) clrpInv(y.star)
y |
shifted data matrix (e.g., output of shiftp) |
p |
weights (should not be closed) |
y.star |
a data matrix that represents data transformed by |
Note that this function will close the dataset y
to 1.
Inverting clrp
transform should be followed by shiftpInv
to return to unshifted original compositoin (see examples).
matrix
Justin Silverman
J. J. Egozcue, V. Pawlowsky-Glahn (2016) Changing the Reference Measure in the Simplex and its Weighting Effects. Austrian Journal of Statistics 45(4):25-44
p <- seq(.1,1,by=.2) c <- t(rmultinom(10,100,c(.1,.6,.2,.3,.2))) + 0.65 # add a small pseudocount x <- miniclo(c) y <- shiftp(x, p) y.star <- clrp(y, p) y.star # Untransform data (note use of shiftp and miniclo to return to x) y.closed <- clrpInv(y.star) all.equal(miniclo(shiftpInv(y.closed, p)), x)