fix polarity of polarization error :)

terhorst · terhorst · commit a22941d0b8b1 · 2017-01-26T11:14:02.000-08:00
diff --git a/src/inference_manager.cpp b/src/inference_manager.cpp
@@ -334,9 +334,9 @@ NPopInferenceManager<P>::construct_bins(const double polarization_error)
                 block_key_prob_map m, new_keys;
                 for (const std::pair<block_key, double> &kp : bin_key<P>::run(bk, na))
                 {
-                    new_keys[kp.first] += polarization_error * kp.second;
+                    new_keys[kp.first] += (1. - polarization_error) * kp.second;
                     if (polarization_error > 0)
-                        new_keys[folded_key(kp.first)] += (1. - polarization_error) * kp.second;
+                        new_keys[folded_key(kp.first)] += polarization_error * kp.second;
                 }
                 new_keys = merge_monomorphic(new_keys);
                 for (const std::pair<block_key, double> &kp : new_keys)

Original file line number	Diff line number	Diff line change
`@@ -334,9 +334,9 @@ NPopInferenceManager<P>::construct_bins(const double polarization_error)`
`334`	`334`	`block_key_prob_map m, new_keys;`
`335`	`335`	`for (const std::pair<block_key, double> &kp : bin_key<P>::run(bk, na))`
`336`	`336`	`{`
`337`		`- new_keys[kp.first] += polarization_error * kp.second;`
	`337`	`+ new_keys[kp.first] += (1. - polarization_error) * kp.second;`
`338`	`338`	`if (polarization_error > 0)`
`339`		`- new_keys[folded_key(kp.first)] += (1. - polarization_error) * kp.second;`
	`339`	`+ new_keys[folded_key(kp.first)] += polarization_error * kp.second;`
`340`	`340`	`}`
`341`	`341`	`new_keys = merge_monomorphic(new_keys);`
`342`	`342`	`for (const std::pair<block_key, double> &kp : new_keys)`