Revision e168b07923b144993e99ea7b65180aed9fe3efa8 authored by Matthias Templ on 11 December 2009, 00:00:00 UTC, committed by Gabor Csardi on 11 December 2009, 00:00:00 UTC
1 parent 2d2a8cf
Raw File
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)
}

back to top