Method calculateMedianReadCountPerWindow

src/GenomeCopyNumber.cpp:1004–1024

int GenomeCopyNumber::calculateMedianReadCountPerWindow() {
    vector <float> myValues;
	vector<ChrCopyNumber>::iterator it;
	for ( it=chrCopyNumber_.begin() ; it != chrCopyNumber_.end(); it++ ) {
		if (! (it->getChromosome().find("X")!=string::npos || it->getChromosome().find("Y")!=string::npos))
			for (int i = 0; i< it->getLength(); i++) {
				if (it->getValueAt(i) >0) //non-zero values
                    myValues.push_back(it->getValueAt(i));
			}
	}

    if (myValues.size()==0) {
            cerr << "Error: no windows with reads\n";
            cerr << "Unable to proceed..\n";
            exit(-1);
    }

	float median = get_median(myValues);
	myValues.clear();
	return floor(median);
}

double GenomeCopyNumber::calculateMedianAround (float interval, float around) {