| 2068 | |
| 2069 | |
| 2070 | void getCopyNumbers (float copyNumber, vector <int> & copyNumbers) { |
| 2071 | int lowCopy = int(floor(copyNumber)); |
| 2072 | int highCopy = int(ceil(copyNumber)); |
| 2073 | float minDiff = 0.45; minDiff = 0.23; //TODO : external parameter!! |
| 2074 | if (highCopy==1) { |
| 2075 | copyNumbers.push_back(round_f(copyNumber)); |
| 2076 | return; |
| 2077 | } |
| 2078 | if (copyNumber-lowCopy<minDiff) |
| 2079 | copyNumbers.push_back(lowCopy); |
| 2080 | else if (highCopy-copyNumber<minDiff) |
| 2081 | copyNumbers.push_back(highCopy); |
| 2082 | else { |
| 2083 | copyNumbers.push_back(lowCopy); |
| 2084 | copyNumbers.push_back(highCopy); |
| 2085 | } |
| 2086 | } |
| 2087 | |
| 2088 | void getCopyNumbers (float copyNumber, vector <int> & copyNumbers,int ploidy, bool noisyData) { |
| 2089 | int lowCopy = int(floor(copyNumber)); |