Method calculatePloidy

src/GenomeCopyNumber.cpp:3227–3240

}

void GenomeCopyNumber::calculatePloidy(int minCNAlength) {
	//check if median for copy numbers is calculated
	if (chrCopyNumber_[0].getMedianValues().size()==0)
		calculateCopyNumberMedians(minCNAlength,0,0);
	cout << "..Getting ploidy information\n";
	int ploidy[] = {1,2,3,4,5,6,7,8,9,10};
	vector <double> scores;
	for (int i = 0; i < 6 ; i++) {
		scores.push_back(calculateXiSum(ploidy[i]));
	}
	int i = argmin(scores);
	ploidy_ = ploidy[i];
	cout << "..selected ploidy: " <<ploidy_<<"\n";
}

double GenomeCopyNumber::calculateXiSum(int ploidy) {
	int totalNumberOfBP = 0;