| 400 | |
| 401 | |
| 402 | void ChrCopyNumber::removeLowReadCountWindows(ChrCopyNumber control,const int RCThresh) { |
| 403 | if (length_!=control.getLength()) { |
| 404 | cerr << "Warning: control length is not equal to the sample length for chromosome " << chromosome_ << "\n"; |
| 405 | cerr << "Sample: " << length_ << " windows; control: "<< control.getLength()<<" windows\n"; |
| 406 | return; |
| 407 | } |
| 408 | for (int i = 0; i<length_; i++) { |
| 409 | if (control.getValueAt(i) < RCThresh){ |
| 410 | readCount_[i]=NA; |
| 411 | control.setValueAt(i,0); |
| 412 | } |
| 413 | } |
| 414 | } |
| 415 | |
| 416 | void ChrCopyNumber::removeLowReadCountWindows(const int RCThresh) { |
| 417 | for (int i = 0; i<length_; i++) { |
no test coverage detected