DANPOS3 wiq does not work

78 views
Skip to first unread message

libin...@gmail.com

unread,
Sep 24, 2020, 2:02:58 PM9/24/20
to DANPOS
Hi

I am trying to normalise my replicates with quantile normalisation but the wig.py seems does not work well and I am not expert in bioinformatics so really need your help here

 python /Users/bingnanli/git/DANPOS3/danpos.py wig2wiq  sacCer3.chrom.size.Roman test.wig

command:
python /Users/git/DANPOS3/danpos.py wig2wiq sacCer3.chrom.size.Roman test.wig

converting test.wig ...
parsing from test.wig
I start from position 1
II start from position 1
Traceback (most recent call last):
  File "/Users/bingnanli/git/DANPOS3/danpos.py", line 1374, in <module>
    elif sys.argv[1]=='wig2wiq':wig2wiq()
  File "/Users/bingnanli/git/DANPOS3/wiq.py", line 185, in wig2wiq
    format_only=True)
  File "/Users/bingnanli/git/DANPOS3/wiq.py", line 94, in qnorwig
    rawsort(ifile=ifile,sort_ofile=opath[:-3]+'sort.wiq',gfile=gfile,format=iformat,step=step,suppress=suppress,buffer=buffer)
  File "/Users/git/DANPOS3/wiq.py", line 13, in rawsort
    wg=Wig(file=ifile,gfile=gfile,step=step,suppress=suppress)
  File "/Users/git/DANPOS3/wig.py", line 26, in __init__
    if file !="": self.load(file,gfile=gfile,suppress=suppress)
  File "/Users/git/DANPOS3/wig.py", line 995, in load
    self.loadFixed(file,gfile=gfile,suppress=suppress)
  File "/Users/git/DANPOS3/wig.py", line 1056, in loadFixed
    if instep<self.step:nplst[functions.div(pos*instep,self.step)]+=lst[pos] ########## add by kaifu on Sep 6, 2012 ##########
IndexError: index 23021 is out of bounds for axis 0 with size 23021

The first part of my test.wig file 

fixedStep chrom=I start=1 step=1 span=1
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0


Chromatin size

I       230218
II      813184
III     316620
IV      1531933
V       576874
VI      270161
VII     1090940
VIII    562643
IX      439888
X       745751
XI      666816
XII     1078177
XIII    924431
XIV     784333
XV      1091291
XVI     948066
M       85779

thanks a lot for you help

Bingnan

Reply all
Reply to author
Forward
0 new messages