Error in rule freec_exome_somatic

[samarth8392]

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"