Question: Question: Question About Mergesegment.Pl In Varscan 2
gravatar for luxin7091012
6.4 years ago by
luxin70910120 wrote:

Hi there, I'm using Varscan2 CNV pipeline to deal with my exome data, it works fine until the last step which needed a perl script "". my command is like:

perl --ref-arm-sizes ./chr_arm_coordinate.txt ./blahblahblah.copynumber The copynumber file was generated by R package DNAcopy, with 10 columns as requires, here's the format:

ID chrom loc.start loc.end num.mark seg.mean bstat pval lcl ucl MetS6 1 762098 1635855 1040 -0.8484 7.868273154 6.30E-13 1635609 1646602 MetS6 1 1645759 1654150 11 -1.402 6.795998163 3.43E-10 1647851 1663751 MetS6 1 1663751 1690511 33 -0.8588 4.428916611 0.000268018 1688073 1721814 Then I got error below:

Use of uninitialized value $stats{"num_merged_events"} in array element at line 182, <gen1> line 3. Use of uninitialized value in split at line 255. Use of uninitialized value $chrom_size in division (/) at line 271. Illegal division by zero at line 271.

two questions: 1: Can somebody give me a correct arm size format.(I use a format gived by "Questions about"like 1 0 125000000 p 1 125000000 249250621 q I got error too ) 2: Can somebody who succeeded run this perl script share all your input and output. If the author can give some detail format, I will be greatly appreciate. Thanks very much for any help.

ADD COMMENTlink written 6.4 years ago by luxin70910120

Do a search on Biostar for "mergeSegment" I'm almost positive that this has been asked and answered before.

ADD REPLYlink written 6.4 years ago by Chris Miller21k
Please log in to add an answer.


Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 2253 users visited in the last hour