https://github.com/cran/robCompositions
Tip revision: e11c86d2552011aab02cc433769f449e3fd7ce93 authored by Matthias Templ on 06 April 2010, 14:00:48 UTC
version 1.4.0
version 1.4.0
Tip revision: e11c86d
ilr.R
`ilr` <-
function(x){
x.ilr=matrix(NA,nrow=nrow(x),ncol=ncol(x)-1)
D=ncol(x)
for (i in 1:ncol(x.ilr)){
x.ilr[,i]=sqrt((D-i)/(D-i+1))*log(((apply(as.matrix(x[,(i+1):D,drop=FALSE]),1,prod))^(1/(D-i)))/(x[,i]))
}
invisible(x.ilr)
}