Method calculateXiSum

src/ChrCopyNumber.cpp:1186–1213

}

double ChrCopyNumber::calculateXiSum(int ploidy, map <float,float> &sds, map <float,float> &meds) {
	double sum = 0;
	const double PI = 3.141592;

	for (int i  = 0; i <= (int)bpfinal_.size();i++) {

		float level = getLevelAt(i,ploidy);
		if (level != NA) {
			double theta_i = meds.find(level)->second;
			//bpLengths_[i];
			//sd_[i];
			double d = medianValues_[i]-theta_i;
			if (fragmentNotNA_lengths_[i] != 0) {
				double sqrtVar = sqrt(float(PI/2./fragmentNotNA_lengths_[i]))*max(sds.find(level)->second,sd_[i]);
				double term = pow(d/sqrtVar,2);
				sum += term;
			}
		}
	}
		/*ofstream myfile;
		 myfile.open ("example.txt");
		 myfile <<"ratio\tmedianValues\n";
		for (int i  = 0; i < (int)ratio_.size();i++) {
			myfile << ratio_[i] << "\t" <<medianProfile_[i] <<"\n";
		}
		myfile.close();*/
	return sum;
}

double ChrCopyNumber::calculateXiSum(int ploidy, map <float,float> &sds) {
	double sum = 0;