(7) wiq

WIQ is a function in DANPOS to do genome wide quantile normalization for wiggle format protein occupancy data.

Example 1. normalize two wiggle format data sets to have the same quantile:

Test data (right click to download and put files in directories as indicated below):

Command:

python danpos.py wiq --buffer_size 50 mm9.chr.sizes.xls rawA.wig:rawB.wig

Example 2. normalize wiggle format data rawA.wig to have the same quantile as in rawB.wig:

Test data (right click to download and put files in directories as indicated below):

Command:

python danpos.py wiq --buffer_size 50 mm9.chr.sizes.xls rawA.wig --reference rawB.wig

Example 3. normalize 10 wiggle format data sets raw1.wig, raw2.wig, … , raw10.wig to have the same quantile as in rawA.wig:

Test data (right click to download and put files in directories as indicated below):

(…) ...

Command:

step1: python danpos.py wig2wiq --buffer_size 50 mm9.chr.sizes.xls rawA.wig

(now we will have a result file rawA.sorted.wiq)

step2: python danpos.py wiq --buffer_size 50 mm9.chr.sizes.xls raw1.wig --reference rawA.sorted.wiq --rformat wiq --rsorted 1

(using a sorted .wiq format data as reference will be 2 fold faster than using a .wig format reference data)

Page updated

Google Sites

Report abuse