-
Notifications
You must be signed in to change notification settings - Fork 2
/
Parameter_cphyperp.h
54 lines (48 loc) · 1.72 KB
/
Parameter_cphyperp.h
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
/*
* DPPDiv version 1.0b source code (git: 9c0ac3d2258f89827cfe9ba2b5038f0f656b82c1)
* Copyright 2009-2011
* Tracy Heath(1,2,3) (NSF postdoctoral fellowship in biological informatics DBI-0805631)
* Mark Holder(1)
* John Huelsenbeck(2)
*
* (1) Department of Ecology and Evolutionary Biology, University of Kansas, Lawrence, KS 66045
* (2) Integrative Biology, University of California, Berkeley, CA 94720-3140
* (3) email: tracyh@berkeley.edu
*
* DPPDiv 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.
*
* This program 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 (the file gpl.txt included with this
* distribution or http://www.gnu.org/licenses/gpl.txt for more
* details.
*
* Some of this code is from publicly available source by John Huelsenbeck
*/
#ifndef PARAMETER_CPHYPERP_H
#define PARAMETER_CPHYPERP_H
class MbRandom;
class Model;
class Cphyperp : public Parameter {
public:
Cphyperp(MbRandom *rp, Model *mp, double ga, double gb, int nn,
int pm, bool fixcp);
~Cphyperp(void);
Cphyperp &operator=(const Cphyperp &c);
void clone(const Cphyperp &c);
double update(double &oldLnL);
void print(std::ostream & o) const;
double lnPrior(void);
std::string writeParam(void);
double getCurrentCP() { return currentCP; }
private:
double gammaAlpha;
double gammaBeta;
double currentCP;
int numNodes;
};
#endif