I am interested in using your package- 'rseg' for differential peak calling between treated and control read files that I have.
I seem to be having some trouble with it:
Although I'm told that everything installed fine it doesn't seem to work:
I have two sam files (which I converted to bed files)
and I also have a file that specifies the size of the chromosomes as follows:
chrom size
chr1 197195432
chr2 181748087
chr3 159599783
chr4 155630120
chr5 152537259
chr6 149517037
chr7 152524553
chr8 131738871
chr9 124076172
chr10 129993255
chr11 121843856
chr12 121257530
chr13 120284312
chr14 125194864
chr15 103494974
chr16 98319150
chr17 95272651
chr18 90772031
chr19 61342430
chrX 166650296
chrY 15902555
chrM 16299
chr13_random 400311
chr16_random 3994
chr17_random 628739
chr1_random 1231697
chr3_random 41899
chr4_random 160594
chr5_random 357350
chr7_random 362490
chr8_random 849593
NOTE: I should mention that in the documentation, it is not at all clear how the actual input files are specified i.e. the treated and control BED files. Should they be precluded by '-input' and -control'?
However while everything seems to be in order I get an error:
rseg-diff -c ~/mm9.genome -o out.bed -i 20 -v -mode 2 test.bed control.bed #-c the last two arguments are the input treated and control files
Error:
[LOADING_DATA] chromosomes
[LOADING_DATA] reads
Segmentation fault (core dumped)