Method calculateSDAndMed

src/GenomeCopyNumber.cpp:3266–3303

}

void GenomeCopyNumber::calculateSDAndMed(int ploidy, map <float,float> &sds,map <float,float> &meds) {

	float value,med;
	vector<ChrCopyNumber>::iterator it;

	//fill the map
	map <float,vector <float> > mymap;
	for ( it=chrCopyNumber_.begin() ; it != chrCopyNumber_.end(); it++ ) {
		for (int i = 0; i<it->getLength(); i++) {
			med = it->getMedianProfileAtI(i);
			float level = round_by_ploidy(med, ploidy);
			value = it->getRatioAtBin(i);
			if (isRatioLogged_&& value!=NA)
                value=pow(2,value);
			if (value != NA) {
				if (mymap.count(level) == 0) {
					vector <float> a;
					mymap.insert ( pair<float,vector <float> >(level,a) );
				}
				mymap.find(level)->second.push_back(value);
			}
		}
	}

	//get median and variance for each level
	map<float,vector <float> >::iterator it2;
	for ( it2=mymap.begin() ; it2 != mymap.end(); it2++ ) {
		float level = (*it2).first;
		float median = get_median((*it2).second);
	//	float mean = get_mean((*it2).second);
		float sdev = sd((*it2).second,median);
		meds.insert(pair<float,float>(level,median));
		sds.insert(pair<float,float>(level,sdev));
		(*it2).second.clear();
	}
	mymap.clear();

}

void GenomeCopyNumber::calculateSDs(int ploidy, map <float,float> &sds) {