diff --git a/rpkm_name.py b/rpkm_name.py index 9363f108584c9986cebdd692f5456c63da384b5b..6f40d4eacf67b05733f7c3f100147e1b62d42686 100755 --- a/rpkm_name.py +++ b/rpkm_name.py @@ -50,7 +50,7 @@ def rpkm(peak_file,aln_file,exp_name,columns): columns.append(exp_name) ## RPKM = numReads / (geneLength/1000 * totalNumReads/1,000,000 ) peak_counts = peak_file.multi_bam_coverage(bams=[aln_file]) - total_counts = reduce(lambda x, y: x + y, [ int(l.rstrip('\n').split('\t')[2]) for l in pysam.idxstats(aln_file)]) + total_counts = reduce(lambda x, y: x + y, [ int(l.rstrip('\n').split('\t')[2]) for l in pysam.idxstats(aln_file,split_lines=True)]) rpkm = peak_counts.each(normalized_to_length, 4, float(math.pow(10,9))/total_counts).saveas("test.bed") rpkm_df = rpkm.to_dataframe() #os.remove('test.bed')