source: trunk/src/Merge.h @ 81

// $Id: Merge.h 81 2004-05-26 13:30:01Z cecilia $

#ifndef _theplu_cpptools_merge_  //vad är detta?
#define _theplu_cpptools_merge_ 

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

#include "vector.h"
#include "matrix.h"

namespace theplu {
namespace cpptools {

  using namespace std;

/**
   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
*/
  class Merge
  {
  public:

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

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

    ///
    /// Calculate new weighted data values.
    ///
    // modify data_ and error_
    void do_something(void);
    
    ///
    /// Return matrix
    ///
    /// @return Const reference to the data matrix.
    ///
    inline const gslapi::matrix& data_matrix(void) const { return data_; }
    inline const gslapi::matrix& error_matrix(void) const { return error_; }

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

    // or gslapi::matrix& ?
    gslapi::matrix calculate_weights(void) const;
    void Merge::merging(vector<u_int>, gslapi::matrix&);
  };

}} // of namespace cpptools and namespace theplu

#endif