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
|
// $Id: bblEMSeperate.cpp 962 2006-11-07 15:13:34Z privmane $
#include "bblEM.h"
#include "bblEMSeperate.h"
#include "logFile.h"
//#define VERBOS
bblEMSeperate::bblEMSeperate(vector<tree>& et,
const vector<sequenceContainer>& sc,
const vector<stochasticProcess> &sp,
const vector<Vdouble *> * weights,
const int maxIterations,
const MDOUBLE epsilon,
const MDOUBLE tollForPairwiseDist) {
MDOUBLE newL =0;
for (int i=0; i < et.size(); ++i) {
#ifdef VERBOS
LOG(5,<<" OPTIMIZING GENE "<<i<<" ... "<<endl);
#endif
bblEM bblEM1(et[i],sc[i],sp[i],(weights?(*weights)[i]:NULL),maxIterations,epsilon);
MDOUBLE resTmp = bblEM1.getTreeLikelihood();
#ifdef VERBOS
LOG(5,<<" GENE "<<i<<" LOG-L = "<< resTmp<<endl);
#endif
newL += resTmp;
}
_treeLikelihood = newL;
}
|
