BIE Overview Quick Start Theoretical overview Available methods Parallel chains Computation engine Command Line Interface Graphical User Interface Assigning Output Visualization tool Data handling Software technology Parallel debugging Results to date Recent developments Links User guide Future features Project goals Project Team Copyright and license

/home/weinberg/src/BIE/include/MixturePrior.h // This is really -*- C++ -*- #ifndef MixturePrior_h #define MixturePrior_h #include <vector> #include <string> #include "Prior.h" #include "MHWidth.h" #include "Dirichlet.h" #include "Chain.h" #include "Poisson.h" #include "Serializable.h" namespace BIE { //+ CLICLASS MixturePrior class MixturePrior : public Prior { protected: bool Ordered; int Iorder; int Nmix; vector<double> ret; double Alpha; vector<Dirichlet*> mix; Uniform * unit; Uniform * useg; Normal * normal; double pmean; Poisson * pois; double GenGamma(double a); double GenBeta(double a, double b); double minW, maxW; public: static unsigned ITMAX; static double min_prob; MixturePrior(); MixturePrior(MixturePrior*); MixturePrior(int nmix, string pfile); MixturePrior(int nmix, double alpha, string pfile); MixturePrior(int nmix, double alpha, double pmean, string pfile); ~MixturePrior(); //+ CLIMETHOD void setMaxIter int void setMaxIter(int i) { ITMAX = i; } //+ CLIMETHOD void setMinMaxWeight double double void setMinMaxWeight(double minw, double maxw) { minW = minw; maxW = maxw; } void getMinMaxWeight(double& minw, double& maxw) { minw = minW; maxw = maxW; } virtual MixturePrior* New() = 0; virtual double PDF(State&) = 0; virtual double logPDF(State&) = 0; virtual double logPDFMarginal(int M, int n, vector<double>& V) = 0; virtual vector<double> lower(void) = 0; virtual vector<double> upper(void) = 0; /* virtual vector<double> Mean(int m) = 0; virtual vector<double> StdDev(int m) = 0; virtual vector<double> Moments(int, int) = 0; */ virtual int SampleM(void) { int m = Nmix; if (pois) { do { m = static_cast<int>(floor((*pois)()+1.0e-6)); } while (m>Nmix || m<1); } return m; } virtual vector<double> Sample(void) { return Sample(SampleM()); } virtual vector<double> Sample(int m) = 0; virtual void SamplePrior(int&, vector<double>&, vector< vector<double> >&) = 0; virtual void SamplePrior(int, int, vector<double>&) = 0; virtual void SampleProposal(Chain &ch, MHWidth* width); virtual double BirthWeight(int M); virtual double BirthWeightPDF(int M, double p); virtual int MinC() { return 1; } virtual int MaxC() { return Nmix; } //+ CLIMETHOD void setOrdering int virtual void setOrdering(int icomp) { Ordered = true; Iorder = icomp; } //+ CLIMETHOD void unsetOrdering virtual void unsetOrdering() { Ordered = false; } virtual bool Order() { return Ordered; } virtual int Oindx() { return Iorder; } // AUTO GENERATED BY ../persistence/autopersist.py private: friend class boost::serialization::access; template<class Archive> void serialize(Archive & ar, const unsigned int version) { this->pre_serialize(ar, version); try { ar & BOOST_SERIALIZATION_BASE_OBJECT_NVP(Prior); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar & BOOST_SERIALIZATION_NVP(Ordered); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar & BOOST_SERIALIZATION_NVP(Iorder); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar & BOOST_SERIALIZATION_NVP(Nmix); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar & BOOST_SERIALIZATION_NVP(ret); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar & BOOST_SERIALIZATION_NVP(Alpha); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar & BOOST_SERIALIZATION_NVP(mix); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar & BOOST_SERIALIZATION_NVP(unit); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar & BOOST_SERIALIZATION_NVP(useg); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar & BOOST_SERIALIZATION_NVP(normal); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar & BOOST_SERIALIZATION_NVP(pmean); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar & BOOST_SERIALIZATION_NVP(pois); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar & BOOST_SERIALIZATION_NVP(ITMAX); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar & BOOST_SERIALIZATION_NVP(min_prob); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } this->post_serialize(ar, version); } }; } // namespace BIE BIE_CLASS_EXPORT_KEY(BIE::MixturePrior) #endif