Method EstimateSubclonalPopulation

src/SeekSubclones.cpp:197–245

void SeekSubclones::EstimateSubclonalPopulation(vector <float> data, float threshold, int ploidy_)
{
    float sumtmp = 0;
    threshold = threshold*ploidy_;
    int countData = 0;
    for (unsigned int i = 0; i < data.size(); i++) {
        if (data[i]!=NA) {
            sumtmp += data[i];
            countData++;
        }

    }
    if (countData==0) {return;}
    float mean = ploidy_*sumtmp/countData;
    if (mean > threshold)
        {
        int i =1;
        float pop = 1;
        int iter = 0;
        while (pop > minimal_pop_ && iter < 100)
            {
            if ((((mean - threshold)/(i)) > minimal_pop_) && (((mean - threshold)/(i)) < 1) ) // && (threshold + i != ploidy_)
                {
                copynumber_.push_back(threshold + i);
                population_.push_back((mean - threshold)/(i));
                }
            pop = (mean - threshold)/(i);
            i++;
            iter++;
            }
        }
    else if (mean < threshold)
        {
        int i =1;
        float pop = 1;
        int iter = 0;
        while (pop > minimal_pop_ && iter < 100 && (threshold - i >= 0))
            {
            if ((((-mean + threshold)/(i)) > minimal_pop_) && (((-mean + threshold)/(i)) < 1) ) //GOT IT EXCEPT for "&& (threshold - i != ploidy_)" SO I REMOVED IT
                {
                copynumber_.push_back(threshold - i);
                population_.push_back((-mean + threshold)/(i));
                }
            pop = (-mean + threshold)/(i);
            i++;
            iter++;
            }
        }
}