relatedness_data.h

#ifndef _RELATEDNESS_DATA_H_
#define _RELATEDNESS_DATA_H_

#include <cstring>

#include <iostream>
#include <cfloat>
#include <iomanip>
#include <vector>
#include <sstream>

#include "typedef.h"
#include "data.h"

#define E_LIM 25

class Relatedness : public Data{ 
private:
    void write (std::ostream&) const;   
    void read (std::istream&);      
    static const Registration registered;
    static Data * create(const std::vector <std::string> & Columns){
        return new Relatedness(Columns);
    }
    bool likelihoods_;          
public:
    char delim; 

    std::string X_str_; 
    std::string Y_str_; 

    id0_t X_, Y_;   

    id1_t sites;    

    float_t e_X_[E_LIM], e_X_ll;
    float_t e_Y_[E_LIM], e_Y_ll;
    /*
     * See ?
     */
    float_t f_X_, f_X_ll;
    float_t f_Y_, f_Y_ll;
    float_t theta_XY_, theta_XY_ll;
    float_t gamma_XY_, gamma_XY_ll;
    float_t gamma_YX_, gamma_YX_ll;
    float_t delta_XY_, delta_XY_ll;
    float_t Delta_XY_, Delta_XY_ll;
    float_t null_ll_, max_ll_;

    Relatedness();  
    Relatedness(const std::vector <std::string> &);
    Relatedness(const std::string &, const std::string &);        

    std::string header(void) const;
    size_t size(void) const;

    void set_X_name(const std::string &);
    void set_X_name(const id0_t &);
    void set_Y_name(const std::string &);
    void set_Y_name(const id0_t &);

    void clear(void); 
    void zero(void);  

    static const std::string file_name; 
    static const std::string table_name;    

    const std::string get_file_name(void) const;
    const std::string get_table_name(void) const;

    static const bool binary;

    const bool get_binary() const;

    Relatedness& operator=(const Relatedness &rhs);
};

#endif