| 1648 | } |
| 1649 | |
| 1650 | float GenomeCopyNumber::getMedianCopyNumber() { |
| 1651 | vector<float>selectedValues; |
| 1652 | vector<ChrCopyNumber>::iterator it; |
| 1653 | for ( it=chrCopyNumber_.begin() ; it != chrCopyNumber_.end(); it++ ) { |
| 1654 | vector<float> chr_values = it->getValues(); |
| 1655 | for (int i = 0; i<(int)chr_values.size(); i++) |
| 1656 | //if (chr_values[i]!=0) |
| 1657 | selectedValues.push_back(chr_values[i]); |
| 1658 | } |
| 1659 | float median = get_median(selectedValues); |
| 1660 | //float mean = get_mean(selectedValues); |
| 1661 | selectedValues.clear(); |
| 1662 | return median; |
| 1663 | } |
| 1664 | |
| 1665 | ChrCopyNumber GenomeCopyNumber::getChrCopyNumber(std::string const& chr) { |
| 1666 | int index = findIndex(chr); |
no test coverage detected