g.rowMeans {philr} | R Documentation |
Calculates weighted geometric mean (see references). Note
if p=rep(1, nrow(y))
(default) then this is just the
geometric mean of rows.
g.rowMeans(y, p = rep(1, nrow(y)))
y |
shifted data matrix (e.g., output of shiftp) |
p |
weights (should not be closed) |
vector (weighted geometric mean of rows)
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) philr:::g.rowMeans(y, p)