diff --git a/maxatac/analyses/normalize.py b/maxatac/analyses/normalize.py index 1aa0809..76b3353 100644 --- a/maxatac/analyses/normalize.py +++ b/maxatac/analyses/normalize.py @@ -30,7 +30,7 @@ def run_normalization(args): 3) Find the genomic min and max values by looping through each chromosome 5) Loop through each chromosome and minmax normalize the values based on the genomic values. - :param args: signal, output_dir, chromosome_sizes, chromosomes, max_percentile, method, blacklist + :param args: signal, output_dir, chromosome_sizes, chromosomes, max_percentile, method, blacklist, max_zooms :return: A minmax normalized bigwig file """ @@ -83,7 +83,7 @@ def run_normalization(args): with pyBigWig.open(args.signal) as input_bw, pyBigWig.open(output_filename, "w") as output_bw: header = [(x, chromosome_length_dictionary[x]) for x in sorted(args.chromosomes)] - output_bw.addHeader(header) + output_bw.addHeader(header, maxZooms = args.max_zooms) # For every chromosome in header, perform normalization for chrom_name, chrom_length in header: @@ -115,7 +115,7 @@ def run_normalization(args): ends=chrom_length, span=1, step=1, - values=normalized_signal.tolist() + values=normalized_signal.astype(np.float16).tolist() ) # Measure time of averaging