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