From ad9d3e1669333b1789eb35ae51941a5c23de2c90 Mon Sep 17 00:00:00 2001
From: LauraU123 <laura.urbanska@stud.unibas.ch>
Date: Thu, 10 Nov 2022 11:09:44 +0100
Subject: [PATCH] updated poisson sampling script
---
scripts/poisson_sampling.py | 6 +++---
1 file changed, 3 insertions(+), 3 deletions(-)
diff --git a/scripts/poisson_sampling.py b/scripts/poisson_sampling.py
index 737cc82..92537e8 100644
--- a/scripts/poisson_sampling.py
+++ b/scripts/poisson_sampling.py
@@ -16,13 +16,13 @@ output: csv file with gene id and count
'''
-def transcript_sampling(total_transcript_number, csv_file, output_csv):
- df = pd.read_csv(csv_file)
+def transcript_sampling(total_transcript_number, tsv_file, output_csv):
+ df = pd.read_csv(tsv_file, sep='\t', lineterminator='\n', names=["id", "level"]))
levels = []
for expression_level in df['level']:
- poisson_sampled = np.random.poisson(total_transcript_number/expression_level)
+ poisson_sampled = np.random.poisson(total_transcript_number*expression_level)
levels.append(poisson_sampled)
transcript_numbers = pd.DataFrame({'id': df['id'],'count': levels})
--
GitLab