PolySNP.hpp

Go to the documentation of this file.

/*

Copyright (C) 2003-2009 Kevin Thornton, krthornt[]@[]uci.edu

Remove the brackets to email me.

This file is part of libsequence.

libsequence 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.

libsequence 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
long with libsequence.  If not, see <http://www.gnu.org/licenses/>.

*/

#ifndef _POLYSNP_H_
#define _POLYSNP_H_

#include <vector>
#include <memory>
#include <limits>
#include <boost/utility.hpp>
namespace Sequence 
  {
  class PolyTable;
  class _PolySNPImpl;
  class PolySNP : boost::noncopyable
    {
    private:
    protected:
      std::auto_ptr<_PolySNPImpl> rep;
      void DepaulisVeuilleStatistics (void);
      virtual void WallStats(void);
      //various things one needs to know to calculate the summary statistics
      double a_sub_n (void) ;
      double a_sub_n_plus1 (void) ;
      double b_sub_n (void) ;
      double b_sub_n_plus1 (void) ;
      double c_sub_n (void) ;
      double d_sub_n (void) ;
    public:
      explicit PolySNP (const Sequence::PolyTable * data, bool haveOutgroup = false,
                        unsigned outgroup = 0, bool totMuts = true);
      virtual ~ PolySNP (void);
      //estimators of 4Nu
      virtual double ThetaPi (void);                           //Nucleotide diversity (Tajima 1983)
      virtual double ThetaW (void);                            //Watterson's (1975) Theta
      virtual double ThetaH (void);                            //Theta from homozygosity (Fay and Wu 2001)
      virtual double ThetaL (void);                            //A variant on Fay and Wu's H
      //variances of estimators of 4Nu
      double VarPi (void);
      double StochasticVarPi(void);
      double SamplingVarPi (void);
      double VarThetaW (void);
      //calculate various numbers related to polymorphism
      unsigned NumPoly (void);                                 //number of polymorphic sites in data
      virtual unsigned NumMutations (void);                    //number of inferred mutations 
      virtual unsigned NumSingletons (void);                   //number of mutants at frequency 1
      virtual unsigned NumExternalMutations (void);            //number of derived mutations 
      //summary statistics of the site frequency spectrum
      virtual double TajimasD (void);                          //Tajima's (1989) D
      virtual double Hprime (bool likeThorntonAndolfatto = false);  //A normalized statistic 
                                                                    //related to Fay and Wu's H
      virtual double Dnominator (void);                        //Denominator of Tajima's D
      virtual double FuLiD (void);                             //Fu & Li's (1996) D
      virtual double FuLiF (void);                             //Fu & Li's (1996) F
      virtual double FuLiDStar (void);                         //Fu & Li's (1996) D*
      virtual double FuLiFStar (void);                         //Fu & Li's (1996) F*
      //summary statistics of haplotypes
      double DandVH (void);                                    //Depaulis & Veuille (1998) Haplotype diversity
      unsigned DandVK (void);                                  //Depaulis & Veuille (1998) number of haplotypes
      virtual double WallsB(void);                             //Jeff Wall's (2000) B  statistic
      virtual unsigned WallsBprime(void);                      //Jeff Wall's (2000) B' statistic
      virtual double WallsQ(void);                             //Jeff Wall's (2000) Q  statistic
      //recombination
      double HudsonsC (void);                                  //Dick Hudson's (1987) Chat = 4Nr
      virtual unsigned Minrec (void);                          //Hudson & Kaplan's (1985) min # of recombination events
      std::vector < std::vector < double > >
      Disequilibrium( const unsigned & mincount = 1,
                      const double & max_marker_distance = std::numeric_limits<double>::max()) ;                 //summary stats of pairwise LD
    };
}
#endif

---