Mercurial > repos > greg > repmatch_gff3
comparison repmatch_gff3_util.py @ 5:53cbf79396d7 draft
Uploaded
| author | greg |
|---|---|
| date | Sat, 28 Nov 2015 14:54:22 -0500 |
| parents | 94e374d13037 |
| children | d10ae3aeebc8 |
comparison
equal
deleted
inserted
replaced
| 4:1cea96c461a2 | 5:53cbf79396d7 |
|---|---|
| 19 pyplot.rc('ytick.major', size=10.00) | 19 pyplot.rc('ytick.major', size=10.00) |
| 20 pyplot.rc('lines', linewidth=4.00) | 20 pyplot.rc('lines', linewidth=4.00) |
| 21 pyplot.rc('axes', linewidth=3.00) | 21 pyplot.rc('axes', linewidth=3.00) |
| 22 pyplot.rc('font', family='Bitstream Vera Sans', size=32.0) | 22 pyplot.rc('font', family='Bitstream Vera Sans', size=32.0) |
| 23 | 23 |
| 24 PLOT_FORMATS = ['png', 'pdf', 'svg'] | |
| 25 COLORS = 'krb' | 24 COLORS = 'krb' |
| 26 | 25 |
| 27 | 26 |
| 28 class Replicate(object): | 27 class Replicate(object): |
| 29 | 28 |
| 252 | 251 |
| 253 def gff_row(cname, start, end, score, source, type='.', strand='.', phase='.', attrs={}): | 252 def gff_row(cname, start, end, score, source, type='.', strand='.', phase='.', attrs={}): |
| 254 return (cname, source, type, start, end, score, strand, phase, gff_attrs(attrs)) | 253 return (cname, source, type, start, end, score, strand, phase, gff_attrs(attrs)) |
| 255 | 254 |
| 256 | 255 |
| 257 def get_temporary_plot_path(plot_format): | 256 def get_temporary_plot_path(): |
| 258 """ | 257 """ |
| 259 Return the path to a temporary file with a valid image format | 258 Return the path to a temporary file with a valid image format |
| 260 file extension that can be used with bioformats. | 259 file extension that can be used with bioformats. |
| 261 """ | 260 """ |
| 262 tmp_dir = tempfile.mkdtemp(prefix='tmp-repmatch-') | 261 tmp_dir = tempfile.mkdtemp(prefix='tmp-repmatch-') |
| 263 fd, name = tempfile.mkstemp(suffix='.%s' % plot_format, dir=tmp_dir) | 262 fd, name = tempfile.mkstemp(suffix='.pdf', dir=tmp_dir) |
| 264 os.close(fd) | 263 os.close(fd) |
| 265 return name | 264 return name |
| 266 | 265 |
| 267 | 266 |
| 268 def process_files(dataset_paths, galaxy_hids, method, distance, step, replicates, up_limit, low_limit, output_files, | 267 def process_files(dataset_paths, galaxy_hids, method, distance, step, replicates, up_limit, low_limit, output_files, |
| 269 plot_format, output_summary, output_orphan, output_detail, output_key, output_histogram): | 268 output_summary, output_orphan, output_detail, output_key, output_histogram): |
| 270 output_histogram_file = output_files in ["all"] and method in ["all"] | 269 output_histogram_file = output_files in ["all"] and method in ["all"] |
| 271 if len(dataset_paths) < 2: | 270 if len(dataset_paths) < 2: |
| 272 return | 271 return |
| 273 if method == 'all': | 272 if method == 'all': |
| 274 match_methods = METHODS.keys() | 273 match_methods = METHODS.keys() |
| 282 step, | 281 step, |
| 283 replicates, | 282 replicates, |
| 284 up_limit, | 283 up_limit, |
| 285 low_limit, | 284 low_limit, |
| 286 output_files, | 285 output_files, |
| 287 plot_format, | |
| 288 output_summary, | 286 output_summary, |
| 289 output_orphan, | 287 output_orphan, |
| 290 output_detail, | 288 output_detail, |
| 291 output_key, | 289 output_key, |
| 292 output_histogram) | 290 output_histogram) |
| 293 if output_histogram_file: | 291 if output_histogram_file: |
| 294 tmp_histogram_path = get_temporary_plot_path(plot_format) | 292 tmp_histogram_path = get_temporary_plot_path() |
| 295 frequency_histogram([stat['distribution'] for stat in [statistics]], | 293 frequency_histogram([stat['distribution'] for stat in [statistics]], |
| 296 tmp_histogram_path, | 294 tmp_histogram_path, |
| 297 METHODS.keys()) | 295 METHODS.keys()) |
| 298 shutil.move(tmp_histogram_path, output_histogram) | 296 shutil.move(tmp_histogram_path, output_histogram) |
| 299 | 297 |
| 300 | 298 |
| 301 def perform_process(dataset_paths, galaxy_hids, method, distance, step, num_required, up_limit, low_limit, output_files, | 299 def perform_process(dataset_paths, galaxy_hids, method, distance, step, num_required, up_limit, low_limit, output_files, |
| 302 plot_format, output_summary, output_orphan, output_detail, output_key, output_histogram): | 300 output_summary, output_orphan, output_detail, output_key, output_histogram): |
| 303 output_detail_file = output_files in ["all"] and output_detail is not None | 301 output_detail_file = output_files in ["all"] and output_detail is not None |
| 304 output_key_file = output_files in ["all"] and output_key is not None | 302 output_key_file = output_files in ["all"] and output_key is not None |
| 305 output_orphan_file = output_files in ["all", "simple_orphan"] and output_orphan is not None | 303 output_orphan_file = output_files in ["all", "simple_orphan"] and output_orphan is not None |
| 306 output_histogram_file = output_files in ["all"] and output_histogram is not None | 304 output_histogram_file = output_files in ["all"] and output_histogram is not None |
| 307 replicates = [] | 305 replicates = [] |
| 456 orphan.midpoint+1, | 454 orphan.midpoint+1, |
| 457 orphan.value, | 455 orphan.value, |
| 458 orphan.distance, | 456 orphan.distance, |
| 459 orphan.replicate.id)) | 457 orphan.replicate.id)) |
| 460 if output_histogram_file: | 458 if output_histogram_file: |
| 461 tmp_histogram_path = get_temporary_plot_path(plot_format) | 459 tmp_histogram_path = get_temporary_plot_path() |
| 462 frequency_histogram([freq], tmp_histogram_path) | 460 frequency_histogram([freq], tmp_histogram_path) |
| 463 shutil.move(tmp_histogram_path, output_histogram) | 461 shutil.move(tmp_histogram_path, output_histogram) |
| 464 return {'distribution': freq} | 462 return {'distribution': freq} |
