Varscan Copy Number Calling Workflow: Using Mergesegments.Pl After Cbc
1
2
Entering edit mode
8.7 years ago
enabieva ▴ 50

I am trying to identify copy number variation in tumor/normal data using the workflow on the VarScan page. I got as far as the last step, where I got stuck on input to mergeSegments.pl mergeSegments.pl says that it requires the following fields in the input:

USAGE: mergeSegments.pl [segments] OPTIONS
segments - A segments file with p-values from the DNAcopy library. This should be tab- or space-delimited
with a header and the following columns: ID, chrom, loc.start, loc.end, num.mark, seg.mean, bstat, pval, lcl, ucl.

Meanwhile, the R script given on the workflow page produces a file with only the following fields:

ID      chrom   loc.start       loc.end num.mark        seg.mean

(even if I write out all columns of segs$output, where segs is the result of calling the segment function in DNAcopy) My question is: where and how do I get the remaining coumns (bstat, pval, lcl, ucl) for input to mergeSegments.pl? varscan copynumber • 4.9k views ADD COMMENT 4 Entering edit mode 8.6 years ago Irsan ★ 7.4k When you run the segments.p method from DNAcopy on a DNACopy-class object (that is returned by segments(...)) a data frame containing p-values is returned segments.p(segs, ngrid=100, tol=1e-6, alpha=0.05, search.range=100, nperm=1000) also see reference manual for DNACopy page 13-14 ADD COMMENT 0 Entering edit mode After I ran the segments.p. I tried to run the mergeSegments.pl. Use of uninitialized value$input in <handle> at ../../mergeSegments.pl line 446. readline() on unopened filehandle at ../../mergeSegments.pl line 446. Can't use an undefined value as a symbol reference at ../../mergeSegments.pl line 456.

Anybody got a similar error?

0
Entering edit mode

https://sourceforge.net/p/varscan/discussion/1073559/thread/663db566/ this helped me, just make a ref-arm-sizes file from the data in the last comment (note that the file should be tab-delimited.) which you can use in the --ref-arm-sizes option in mergeSegments.pl I personally had to add chr before the chrom. number