https://github.com/cran/robCompositions
Tip revision: 54cbcf1d6d3ecbca1153e8f00a804e837237f09d authored by Matthias Templ on 07 February 2014, 00:00:00 UTC
version 1.7.0
version 1.7.0
Tip revision: 54cbcf1
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)
}