Skip to content
Snippets Groups Projects
Commit 0abbc46e authored by Pierre Neveu's avatar Pierre Neveu
Browse files

Update scection.py

parent dc80e121
Branches master
No related tags found
No related merge requests found
......@@ -18,7 +18,8 @@ def scection_step(data, indices_to_consider, min_expression,expression_range, nm
nb_samples=len(indices_to_consider)
data_to_consider=data[:,indices_to_consider]
# keep only genes that have in one sample and a fold change greater than expression_range across samples
# keep only genes that have a minimum expression greater than min_expression in at least one sample
# and a fold change greater than expression_range across samples
genes_indices=((amax(data_to_consider,axis=1)>=log(min_expression)/log(2))
& (amax(data_to_consider,axis=1)-amin(data_to_consider,axis=1)>=log(expression_range)/log(2))).nonzero()[0]
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment