source: trunk/lib/utility/Merge.h @ 376

// $Id: Merge.h 376 2005-08-07 22:32:48Z jari $

#ifndef _theplu_utility_merge_  
#define _theplu_utility_merge_ 

#include <c++_tools/gslapi/matrix.h>

#include <iostream>
#include <string.h>
#include <utility>
#include <vector>


namespace theplu {
namespace utility {

  ///
  /// A merger of datapoints on annotations.  Calculates weighted
  /// average (m) of all (i) data points with same annotation:
  /// m=sum(w_i*x_i)sum(w_i) and uncertainty U of new data point:
  /// U=1/(sum(1/u_i))+sum(w_i^2*(x_i-m)^2)/(sum(w_i))^2 w_i is
  /// calculated from errors u_i Missing values have error < 0
  ///

  class Merge
  {
  public:

    ///
    /// The default constructor. Not implemented.
    ///
    Merge(void);

    ///
    /// istream Constructor, arguments: data(matrix), weights(matrix), 
    /// annotations(vector)
    ///
    Merge(std::istream& data, std::istream& error, std::istream& annotations);

    ///
    /// matrix Constructor, arguments: data(matrix), weights(matrix), 
    /// annotations(vector)
    ///
    Merge(gslapi::matrix& data, gslapi::matrix& error,
          std::vector<std::string>& annotations);
    
    ///
    /// Calculate new weighted data and error values.
    ///
    void do_something(void);
    
    ///
    /// @return Const reference to the data matrix.
    ///
    inline const gslapi::matrix& data_matrix(void) const { return data_; }

    ///
    /// @return Const reference to the error matrix.
    ///
    inline const gslapi::matrix& error_matrix(void) const { return error_; }

  private:
    gslapi::matrix data_;
    gslapi::matrix error_;
    std::vector<std::string> annots_;
    u_int nof_assays_;

    void calculate_weights(gslapi::matrix&);
    void merging(std::vector<u_int>, gslapi::matrix&);
  };

}} // of namespace utility and namespace theplu

#endif