File tree Expand file tree Collapse file tree 1 file changed +2
-6
lines changed Expand file tree Collapse file tree 1 file changed +2
-6
lines changed Original file line number Diff line number Diff line change @@ -35,23 +35,19 @@ def preprocess_data(adata, params):
3535
3636 sc .pp .filter_genes (adata , min_cells = 30 )
3737 if params .data_layer in adata .layers .keys ():
38+ logger .info (f'Using data layer: { params .data_layer } ...' )
3839 adata .X = adata .layers [params .data_layer ]
3940 else :
4041 raise ValueError (f'Invalid data layer: { params .data_layer } , please check the input data.' )
4142
4243 if params .data_layer in ['count' , 'counts' ]:
43-
44- adata .X = adata .layers [params .data_layer ].copy ()
45- # HVGs based on count
44+ # HVGs based on count
4645 sc .pp .highly_variable_genes (adata ,flavor = "seurat_v3" ,n_top_genes = params .feat_cell )
4746 # Normalize the data
4847 sc .pp .normalize_total (adata , target_sum = 1e4 )
4948 sc .pp .log1p (adata )
5049
5150 elif params .data_layer in adata .layers .keys ():
52- logger .info (f'Using { params .data_layer } data...' )
53- adata .X = adata .layers [params .data_layer ].copy ()
54- # HVGs based on normalzied data
5551 sc .pp .highly_variable_genes (adata ,flavor = "seurat" ,n_top_genes = params .feat_cell )
5652
5753 return adata
You can’t perform that action at this time.
0 commit comments