CCBR / XAVIER

An easy-to-use, flexible variant calling pipeline for use on the Biowulf cluster at NIH
https://ccbr.github.io/XAVIER/
MIT License
1 stars 3 forks source link

Error in rule freec_exome_somatic #120

Closed samarth8392 closed 1 month ago

samarth8392 commented 1 month ago

Issue

The rule freec_exome_somatic_pass1 fails with the following error:

Loading required package: IRanges
Loading required package: GenomeInfoDb
Error in if (class(resultks) == "try-error") { : 
  the condition has length > 1
Calls: KS
In addition: Warning message:
In ks.test.default(values, score(normals)) :
  p-value will be approximate in the presence of ties
Execution halted

Possible explanation:

In the script assess_signficance.R, on line 51, the code is: if (class(resultks) == "try-error")

but class(resultks) gives the output:

> class(resultks)
[1] "ks.test" "htest"  

Proposed solution:

collapse class(resultks) output into a single string, and then compare.

paste(class(resultks), collapse="_")
[1] "ks.test_htest"