File tree Expand file tree Collapse file tree
Expand file tree Collapse file tree Original file line number Diff line number Diff line change @@ -40,8 +40,6 @@ def setup_cache(self) -> None:
4040
4141 def setup (self , dataset , layer ) -> None :
4242 self .adata = ad .read_h5ad (f"{ dataset } _{ layer } .h5ad" )
43- if "X_pca" not in self .adata .obsm :
44- sc .pp .pca (self .adata )
4543
4644 def time_pca (self , * _ ) -> None :
4745 sc .pp .pca (self .adata , svd_solver = "arpack" )
@@ -69,12 +67,6 @@ def time_regress_out(self, *_) -> None:
6967 def peakmem_regress_out (self , * _ ) -> None :
7068 sc .pp .regress_out (self .adata , ["total_counts" , "pct_counts_mt" ])
7169
72- def time_neighbors (self , * _ ) -> None :
73- sc .pp .neighbors (self .adata , n_neighbors = 10 , n_pcs = 40 )
74-
75- def peakmem_neighbors (self , * _ ) -> None :
76- sc .pp .neighbors (self .adata , n_neighbors = 10 , n_pcs = 40 )
77-
7870 def time_scale (self , * _ ) -> None :
7971 sc .pp .scale (self .adata , max_value = 10 )
8072
Original file line number Diff line number Diff line change @@ -220,7 +220,6 @@ def combat( # noqa: PLR0915
220220 "within-batch variance. Filter these batches before running combat."
221221 )
222222 raise ValueError (msg )
223- n_array = float (sum (n_batches ))
224223
225224 # standardize across genes using a pooled variance estimator
226225 logg .info ("Standardizing Data across genes.\n " )
@@ -273,7 +272,7 @@ def combat( # noqa: PLR0915
273272
274273 # we now apply the parametric adjustment to the standardized data from above
275274 # loop over all batches in the data
276- bayesdata_arr = bayesdata .values
275+ bayesdata_arr = bayesdata .to_numpy ( copy = True )
277276 batch_design_arr = batch_design .values
278277 for j , batch_idxs in enumerate (batch_info .values ()):
279278 # we basically subtract the additive batch effect, rescale by the ratio
You can’t perform that action at this time.
0 commit comments