Commit 6e5435aa authored by Corson N. Areshenkoff's avatar Corson N. Areshenkoff

Logm2 for riemann mean

parent e24620ba
......@@ -23,7 +23,7 @@ spd.mean.p.est <- function(x, p = .5, tol = .01, max.iter = 20){
# Update
H <- Reduce(`+`,
lapply(C, function(i) {
expm2(abs(p) * logm(X %*% i %*% t(X)))
expm2(abs(p) * logm2(X %*% i %*% t(X)))
})) / n
X <- expm2(-phi * logm2(H)) %*% X
......
Markdown is supported
0% or
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment