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
|
#include "gammaDistributionLaguerre.h"
#include "gammaUtilities.h"
#include "logFile.h"
#include <cmath>
gammaDistributionLaguerre::gammaDistributionLaguerre(MDOUBLE alpha,int in_number_of_categories)
: generalGammaDistributionLaguerre(alpha,alpha,in_number_of_categories)
{
}
gammaDistributionLaguerre::gammaDistributionLaguerre(const gammaDistributionLaguerre& other)
: generalGammaDistributionLaguerre(other)
{
}
void gammaDistributionLaguerre::setAlpha(MDOUBLE in_alpha)
{
if (in_alpha == _alpha)
return;
setGammaParameters(categories(), in_alpha);
}
//this function builds the gamma distribution
void gammaDistributionLaguerre::setGammaParameters(int in_number_of_categories, MDOUBLE in_alpha)
{
generalGammaDistributionLaguerre::setGammaParameters(in_number_of_categories, in_alpha, in_alpha);
}
void gammaDistributionLaguerre::change_number_of_categories(int in_number_of_categories)
{
if (in_number_of_categories == categories())
return;
setGammaParameters(in_number_of_categories, _alpha, _alpha);
}
void gammaDistributionLaguerre::setGammaParameters(int numOfCategories ,MDOUBLE alpha, MDOUBLE beta)
{
if (alpha != beta)
errorMsg::reportError("gammaDistributionLaguerre::setGammaParameters : can not set beta because alpha must be equal to beta");
generalGammaDistributionLaguerre::setGammaParameters(numOfCategories, alpha, alpha);
}
|
