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
robVariation.R
`robVariation` <-
function(x){
    rvars <- matrix(0, ncol=ncol(x), nrow=ncol(x))
    for( i in 1:ncol(x)){
      for( j in 1:ncol(x)){
        if( i < j ) rvars[i,j] <- (mad(log(x[,i]/x[,j])))^2
      }
    }
    invisible(rvars) 
}

back to top