analysis.add_a_seq(s0);
analysis.add_a_seq(s1);
analysis.add_a_seq(s2);
- analysis.analyze(4,3);
+ analysis.set_window(4);
+ analysis.set_threshold(3);
+ analysis.analyze();
NwayPaths npath = analysis.paths();
// there should be no paths for these sequences
for (std::list<ConservedPath >::iterator pathz_i = npath.pathz.begin();
Mussa analysis;
analysis.add_a_seq(s0);
analysis.add_a_seq(s1);
- analysis.analyze(4,3);
+ analysis.set_window(4);
+ analysis.set_threshold(3);
+ analysis.analyze();
NwayPaths npath = analysis.paths();
for (std::list<ConservedPath >::iterator pathz_i = npath.pathz.begin();
pathz_i != npath.pathz.end();
mupa_path /= "mck3test.mupa";
Mussa m1;
m1.load_mupa_file( mupa_path );
- m1.analyze(0, 0);
+ m1.analyze();
const NwayPaths& npath = m1.paths();
BOOST_CHECK_EQUAL (npath.path_size(), npath.refined_path_size());
size_t first_refined_size = npath.refined_path_size();
m1.add_a_seq(seq2);
m1.add_a_seq(seq3);
- m1.analyze(10, 8);
+ m1.set_window(10);
+ m1.set_threshold(8);
+ m1.analyze();
m1.set_soft_threshold(10);
m1.nway();
}