diff --git a/src/updateHap.cpp b/src/updateHap.cpp index 4da55d5..38acd5f 100644 --- a/src/updateHap.cpp +++ b/src/updateHap.cpp @@ -206,14 +206,11 @@ void UpdateSingleHap::calcExpectedWsaf( vector & expectedWsaf, vector < void UpdateSingleHap::buildEmission( double missCopyProb ) { - vector noMissProb (this->nLoci_, log(1.0 - missCopyProb)); - vector t1omu = vecSum(llk0_, noMissProb); // t1 one minus u - vector t2omu = vecSum(llk1_, noMissProb); // t2 one minus u - + vector t1omu = vecSum(llk0_, log(1.0 - missCopyProb)); // t1 one minus u + vector t2omu = vecSum(llk1_, log(1.0 - missCopyProb)); // t2 one minus u - vector missProb (this->nLoci_, log(missCopyProb)); - vector t1u = vecSum(llk0_, missProb); - vector t2u = vecSum(llk1_, missProb); + vector t1u = vecSum(llk0_, log(missCopyProb)); + vector t2u = vecSum(llk1_, log(missCopyProb)); assert(emission_.size() == 0 ); for ( size_t i = 0; i < this->nLoci_; i++) {