Mercurial > repos > eric-rasche > circos
comparison unified-histogram.py @ 5:b25b3518c2ce draft
planemo upload for repository https://github.com/TAMU-CPT/galaxy-circos-tool commit 7561690774be81155d79a3c38a77c098c59aa9a2
| author | eric-rasche |
|---|---|
| date | Sat, 10 Jun 2017 12:55:28 -0400 |
| parents | 4ff5ff4c84fa |
| children |
comparison
equal
deleted
inserted
replaced
| 4:4ff5ff4c84fa | 5:b25b3518c2ce |
|---|---|
| 114 else: | 114 else: |
| 115 if MODE == 'histogram': | 115 if MODE == 'histogram': |
| 116 # histogram | 116 # histogram |
| 117 # hs4 0 1999999 5.0000,3.0000,1.0000,19.0000 | 117 # hs4 0 1999999 5.0000,3.0000,1.0000,19.0000 |
| 118 sys.stdout.write(' '.join( | 118 sys.stdout.write(' '.join( |
| 119 (genome, region_start, region_end, ','.join(values)) | 119 (genome, str(region_start), str(region_end), ','.join(map(str, values))) |
| 120 ) + '\n') | 120 ) + '\n') |
| 121 elif MODE == 'heatmap': | 121 elif MODE == 'heatmap': |
| 122 # heatmap | 122 # heatmap |
| 123 # hs1 2000000 3999999 0.0000 id=hs4 | 123 # hs1 2000000 3999999 0.0000 id=hs4 |
| 124 # hs1 4000000 5999999 2.0000 id=hs1 | 124 # hs1 4000000 5999999 2.0000 id=hs1 |
| 127 # hs1 4000000 5999999 0.0000 id=hs4 | 127 # hs1 4000000 5999999 0.0000 id=hs4 |
| 128 # hs1 6000000 7999999 4.0000 id=hs2 | 128 # hs1 6000000 7999999 4.0000 id=hs2 |
| 129 for x in max_idx: | 129 for x in max_idx: |
| 130 if x in data[genome][position]: | 130 if x in data[genome][position]: |
| 131 sys.stdout.write(' '.join( | 131 sys.stdout.write(' '.join( |
| 132 (genome, region_start, region_end, data[genome][position][x], 'id=hm%s' % x) | 132 (genome, str(region_start), str(region_end), data[genome][position][x], 'id=hm%s' % x) |
| 133 ) + '\n') | 133 ) + '\n') |
| 134 else: | 134 else: |
| 135 sys.stdout.write(' '.join( | 135 sys.stdout.write(' '.join( |
| 136 (genome, region_start, region_end, 0.0, 'id=hm%s' % x) | 136 (genome, str(region_start), str(region_end), 0.0, 'id=hm%s' % x) |
| 137 ) + '\n') | 137 ) + '\n') |
| 138 elif MODE == 'line': | 138 elif MODE == 'line': |
| 139 # multiple=False | 139 # multiple=False |
| 140 sys.stdout.write(' '.join( | 140 sys.stdout.write(' '.join( |
| 141 (genome, region_start, region_end, data[genome][position][0]) | 141 (genome, str(region_start), str(region_end), data[genome][position][0]) |
| 142 ) + '\n') | 142 ) + '\n') |
| 143 elif MODE == 'scatter': | 143 elif MODE == 'scatter': |
| 144 # multiple=False | 144 # multiple=False |
| 145 sys.stdout.write(' '.join( | 145 sys.stdout.write(' '.join( |
| 146 (genome, region_start, region_end, data[genome][position][0]) | 146 (genome, str(region_start), str(region_end), data[genome][position][0]) |
| 147 ) + '\n') | 147 ) + '\n') |
| 148 | 148 |
| 149 # Update start of next array | 149 # Update start of next array |
| 150 region_start = position | 150 region_start = position |
| 151 region_end = position | 151 region_end = position |