From e959fde98f95f3595e01490b67892678bbcd1b27 Mon Sep 17 00:00:00 2001
From: Johannes Ranke <jranke@uni-bremen.de>
Date: Wed, 7 May 2014 14:47:28 +0200
Subject: Fork the gmkin GUI from mkin. See ChangeLog for details

---
 R/ilr.R | 51 ---------------------------------------------------
 1 file changed, 51 deletions(-)
 delete mode 100644 R/ilr.R

(limited to 'R/ilr.R')

diff --git a/R/ilr.R b/R/ilr.R
deleted file mode 100644
index 389653e..0000000
--- a/R/ilr.R
+++ /dev/null
@@ -1,51 +0,0 @@
-# $Id$
-
-# Copyright (C) 2012 René Lehmann, Johannes Ranke
-# Contact: jranke@uni-bremen.de
-
-# This file is part of the R package mkin
-
-# mkin is free software: you can redistribute it and/or modify it under the
-# terms of the GNU General Public License as published by the Free Software
-# Foundation, either version 3 of the License, or (at your option) any later
-# version.
-
-# This program is distributed in the hope that it will be useful, but WITHOUT
-# ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or FITNESS
-# FOR A PARTICULAR PURPOSE.  See the GNU General Public License for more
-# details.
-
-# You should have received a copy of the GNU General Public License along with
-# this program. If not, see <http://www.gnu.org/licenses/>
-
-ilr <- function(x) {
-  z <- vector()
-  for (i in 1:(length(x) - 1)) {
-    z[i] <- sqrt(i/(i+1)) * log((prod(x[1:i]))^(1/i) / x[i+1])
-  }
-  return(z)
-}
-
-invilr<-function(x) {
-  D <- length(x) + 1
-  z <- c(x, 0)
-  y <- rep(0, D)
-  s <- sqrt(1:D*2:(D+1))
-  q <- z/s
-  y[1] <- sum(q[1:D])
-  for (i in 2:D) {
-    y[i] <- sum(q[i:D]) - sqrt((i-1)/i) * z[i-1]
-  }
-  z <- vector()
-  for (i in 1:D) {
-    z[i] <- exp(y[i])/sum(exp(y))
-  }
-
-  # Work around a numerical problem with NaN values returned by the above
-  # Only works if there is only one NaN value: replace it with 1
-  # if the sum of the other components is < 1e-10
-  if (sum(is.na(z)) == 1 && sum(z[!is.na(z)]) < 1e-10) 
-    z = ifelse(is.na(z), 1, z)
-
-  return(z)
-}
-- 
cgit v1.2.1