PolySIM.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 _POLYSIM_H_
#define _POLYSIM_H_

#include <Sequence/PolySNP.hpp>

namespace Sequence
  {
  class SimData;
  class PolySIM : public PolySNP
    {
    private:
      //functions for Hudson's Partition Test
      int poly (int *subslist,  int ss,
                int subsize, int subss, int *seq);
      int nextsample (int *subslist, int subsize, int nsam, int seq);
    protected:
      void WallStats(void);
    public:
      explicit PolySIM (const Sequence::SimData * data);
      virtual ~ PolySIM(void);
      //estimators of 4Nu
      double ThetaPi (void);
      double ThetaW (void);
      double ThetaH (void);
      double ThetaL (void);

      //calculate various numbers related to polymorphism
      unsigned NumMutations (void);
      unsigned NumSingletons (void);
      unsigned NumExternalMutations (void);
      //summary statistics of the site frequency spectrum
      double TajimasD (void);
      double Hprime (bool likeThorntonAndolfatto = false);
      double Dnominator (void);
      double FuLiD (void);
      double FuLiF (void);
      double FuLiDStar (void);
      double FuLiFStar (void);
      double WallsB(void);
      unsigned WallsBprime(void);
      double WallsQ(void);
      //Hudson's Haplotype Partition Test
      int HudsonsHaplotypeTest (int subsize, int subss);

      //recombination
      unsigned Minrec (void);
    };
}
#endif

---