Skip to content

Commit 9b45c5c

Browse files
src/analysis/pmd.cpp: fixing a bug initializing transition matrix
1 parent f2ce2d7 commit 9b45c5c

File tree

1 file changed

+4
-10
lines changed

1 file changed

+4
-10
lines changed

src/analysis/pmd.cpp

Lines changed: 4 additions & 10 deletions
Original file line numberDiff line numberDiff line change
@@ -1083,9 +1083,9 @@ main_pmd(int argc, char *argv[]) { // NOLINT(*-avoid-c-arrays)
10831083
try {
10841084
static constexpr auto min_observations_for_inference{100ul};
10851085
static constexpr auto default_self_transition{0.99};
1086+
static constexpr auto default_transition{0.01};
10861087
static constexpr auto default_start_transition{0.5};
10871088
static constexpr auto default_end_transition{1e-10};
1088-
static constexpr auto default_transition{0.01};
10891089

10901090
// magic numbers from paper: highest jaccard index to wgbs
10911091
static constexpr auto bin_size_for_array{1000ul};
@@ -1100,17 +1100,11 @@ main_pmd(int argc, char *argv[]) { // NOLINT(*-avoid-c-arrays)
11001100
static constexpr auto default_bg_alpha{0.95};
11011101
static constexpr auto default_bg_beta{0.05};
11021102
const auto init_trans = [&] {
1103-
return std::vector{
1104-
std::vector<double>(default_self_transition, default_transition),
1105-
std::vector<double>(default_transition, default_self_transition),
1103+
return std::vector<std::vector<double>>{
1104+
{default_self_transition, default_transition},
1105+
{default_transition, default_self_transition},
11061106
};
11071107
};
1108-
// <std::vector<double>>
1109-
// std::vector<std::vector<double>> t(
1110-
// 2, std::vector<double>(2, default_transition));
1111-
// t[0][0] = t[1][1] = default_self_transition;
1112-
// return t;
1113-
// };
11141108

11151109
// NOLINTBEGIN(*-avoid-magic-numbers)
11161110
std::size_t resolution = 500;

0 commit comments

Comments
 (0)