lnmultinomial.h

//General likelihood functions

#ifndef _LNMULTINOMIAL_H_
#define _LNMULTINOMIAL_H_

#include <map>
#include <string>
#include <vector>
#include <cfloat>
#include <algorithm>
#include <iostream>
#include <string.h>
#include <limits>   // std::numeric_limits

#include "typedef.h"
#include "data_types/allele.h"

/* \brief information that should be passed to vcf files.
 */

//* This is pathological. Working to kill it.*/
/* \brief returns the log... 
 */
class lnmultinomial {
private:
    std::vector <float_t> lnfact_vector;        
    float_t *lnp_;                  
    count_t size_;                  
public: 
    lnmultinomial (float_t*, const count_t&);   
    lnmultinomial (const count_t&);         
    lnmultinomial (void);               
    ~lnmultinomial (void);              
    void set (float_t*);                
    void set (float_t, float_t, float_t, float_t);  
    void set(void (*)(const Allele&, float_t *), const Allele&);

    float_t lnprob(const count_t*) ;        
    float_t lnprob_approx(const count_t*);      
    float_t lnfact(const count_t&);         
    float_t lnmultinomcoef(const count_t*);     

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

std::vector <std::pair <count_t, float_t> > sort (const float_t *, const count_t &);
std::vector <std::pair <count_t, count_t> > sort (const count_t *, const count_t &);
#endif