https://github.com/cran/robCompositions
Raw File
Tip revision: a4d717b5864005b0632730e50c7ff6f6b79fa30b authored by Matthias Templ on 06 March 2009, 17:46:31 UTC
version 1.2
Tip revision: a4d717b
aDist.R
`aDist` <-
function(x,y){
    d <- 0
    p <- length(x)
    for(i in 1:(p-1)){
      for(j in (i+1):p){
        d <- d + (log(x[i]/x[j]) - log(y[i]/y[j]))^2
      }
    }
    d=d/p
    sqrt(d)
  }

back to top