properly order the jointseg

mskcc / facets

Algorithm to implement Fraction and Copy number Estimate from Tumor/normal Sequencing.

144 stars 67 forks source link

properly order the jointseg #190

Open ShixiangWang opened 1 year ago

ShixiangWang commented 1 year ago

An error occurred when I run facets with mouse genome data Error in seg.start[seg]:seg.end[seg] : NA/NaN argument.

After I debug, I found that the jointseg in the modified function is not properly ordered. This makes seg.start <- which(diff(c(0, jointseg$seg)) == 1) failed.

Browse[3]> nsegs
[1] 2487
Browse[3]> length(seg.start )
[1] 2484