source: trunk/yat/normalizer/qQuantileNormalizer.cc @ 1736

// $Id: qQuantileNormalizer.cc 1736 2009-01-19 14:26:49Z peter $

/*
  Copyright (C) 2009 Jari Häkkinen, Peter Johansson

  This file is part of the yat library, http://dev.thep.lu.se/yat

  The yat library 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.

  The yat library 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 yat. If not, see <http://www.gnu.org/licenses/>.
*/

#include "qQuantileNormalizer.h"

#include "yat/regression/CSplineInterpolation.h"
#include "yat/statistics/Averager.h"
#include "yat/utility/Vector.h"
#include "yat/utility/VectorBase.h"

#include <algorithm>
#include <cassert>

namespace theplu {
namespace yat {
namespace normalizer {


  void 
  qQuantileNormalizer::Partitioner::init(const utility::VectorBase& sortedvec,
                                         unsigned int N)
  {
    assert(N>1);
    assert(N<=sortedvec.size());
    double range=static_cast<double>(sortedvec.size())/N;
    assert(range);
    unsigned int start=0;
    for (unsigned int i=0; i<N; ++i) {
      unsigned int end = ( i==(N-1) ? sortedvec.size() :
                           static_cast<unsigned int>((i+1)*range) );
      statistics::Averager av;
      for (unsigned int r=start; r<end; ++r)
        av.add(sortedvec(r));
      average_(i) = av.mean();
      index_(i)   = 0.5*(end+start-1);
      start=end;
    }
  }


  const utility::Vector& qQuantileNormalizer::Partitioner::averages(void) const
  {
    return average_;
  }


  const utility::Vector& qQuantileNormalizer::Partitioner::index(void) const
  {
    return index_;
  }


  size_t qQuantileNormalizer::Partitioner::size(void) const
  {
    return average_.size();
  }

}}} // end of namespace normalizer, yat and thep