Question: varscan2: merge the CNA segments using mergeSegments.pl
1
gravatar for Chirag Nepal
3.7 years ago by
Chirag Nepal2.2k
Copenhagen
Chirag Nepal2.2k wrote:

Hi There,

I did copynumber call and performed segmentation using DNAcopy.

 

cnP1 = read.table ("copyNumber_varscanIdhP1.copynumber", header=TRUE)
CNA.object <- CNA( genomdat = cnP1[,7], chrom = cnP1[,1], maploc = cnP1 [,2], data.type = 'logratio')

CNA.object.smoothed     <- smooth.CNA(CNA.object)
CNA.object.smoothed.seg <- segment(CNA.object.smoothed, verbose=0, min.width=2)
segmented = CNA.object.smoothed.seg$output
 
write.table (segmented[,2:6], file="cbsP1", row.names=F, col.names=F, quote=F, sep="\t")

 

Now, i want to merge these output segments:

perl mergeSegments.pl
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.
        OPTIONS:

 

WHAT is the input file here ?

Is the input to the script is output from DNAcopy output.

But the columns are different:

head cbs_P1
chr1    567331    569882    6    -2.3721
chr1    761847    905536    162    -0.1232
chr1    905636    3549990    2846    -0.0335

 

Can you please suggest me to make the input file in compatible format.

thanks !

Edit:  This is how you could get input file (output from DNACopy) to mergeSegments.pl

seg.pvalue <- segments.p(CNA.object.smoothed.seg, ngrid=100, tol=1e-6, alpha=0.05, search.range=100, nperm=1000)
write.table (seg.pvalue, file="blah", row.names=F, col.names=F, quote=F, sep="\t")

 

 

 

 

 

 

cnv mergesegments.pl varscan • 2.1k views
ADD COMMENTlink modified 3.6 years ago • written 3.7 years ago by Chirag Nepal2.2k
2
gravatar for poisonAlien
3.7 years ago by
poisonAlien2.7k
Asgard
poisonAlien2.7k wrote:

Its the output from segments.p function from DNACopy.

ADD COMMENTlink written 3.7 years ago by poisonAlien2.7k
Please log in to add an answer.

Help
Access

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