search

gavinha / TitanCNA

Analysis of subclonal copy number alterations (CNA) and loss of heterozygosity (LOH) in cancer
GNU General Public License v3.0
93 stars 36 forks source link

Check for NAs during parameter estimation #51

Closed chapmanb closed 5 years ago

chapmanb commented 5 years ago

Avoids errors when running copy number estimation:

titan: Parameter estimation
Optimal state path computation: Using 8 cores.
Error in if (objfun[maxFind] >= objfun[i]) { :
  missing value where TRUE/FALSE needed
Calls: outputTitanResults ... runEMclonalCN -> estimateClonalCNParamsMap -> updateParameters

I couldn't determine the underlying source of NAs when reading the code but put in defensive checks which avoids the issue.