https://github.com/cran/robCompositions
Raw File
Tip revision: 165037e6629b3969bb2024b9b74bc0bc78b6e1c4 authored by Matthias Templ on 27 February 2013, 00:00:00 UTC
version 1.6.3
Tip revision: 165037e
robVariation.R
`robVariation` <-
function(x, robust=TRUE){
    rvars <- matrix(0, ncol=ncol(x), nrow=ncol(x))
	if(robust){
    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
      }
    }
	} else{
		for( i in 1:ncol(x)){
			for( j in 1:ncol(x)){
				if( i < j ) rvars[i,j] <- (var(log(x[,i]/x[,j])))
			}
		}		
	}
	rvars[lower.tri(rvars)] <- rvars[upper.tri(rvars)]
    return(rvars) 
}

back to top