| 1625 | } |
| 1626 | |
| 1627 | float GenomeCopyNumber::getMedianRatio() { |
| 1628 | vector<float>selectedValues; |
| 1629 | vector<ChrCopyNumber>::iterator it; |
| 1630 | for ( it=chrCopyNumber_.begin() ; it != chrCopyNumber_.end(); it++ ) { |
| 1631 | vector<float> chr_values = it->getRatio(); |
| 1632 | for (int i = 0; i<(int)chr_values.size(); i++) |
| 1633 | if (chr_values[i]!=NA) |
| 1634 | selectedValues.push_back(chr_values[i]); |
| 1635 | else |
| 1636 | cout << chr_values[i] <<"\n"; |
| 1637 | } |
| 1638 | float median = get_median(selectedValues); |
| 1639 | selectedValues.clear(); |
| 1640 | return median; |
| 1641 | } |
| 1642 | |
| 1643 | float GenomeCopyNumber::calculateNormalizationConstant(GenomeCopyNumber & controlCopyNumber) { |
| 1644 | float median,medianControl; |
nothing calls this directly
no test coverage detected