Dmitry Kobak
Dmitry Kobak
@jlause @LuckyMD Okay, so I am on board with `sc.pp.highly_variable_genes(adata, flavor='pearson_residuals')` and with `sc.pp.normalize_pearson_residuals(adata)` that normalizes all genes and does not change the dimensions (and does not allow it even...
I would use something like ``` adata_for_pca = adata[:,adata.var.highly_variable] sc.pp.normalize_pca_pearson_residuals(adata_for_pca) sc.pp.pca(adata_for_pca, n_comps=n_pc) adata['pca'] = adata_for_pca['pca'] ``` in either bundle.
> Thus, I assume most people that use this recipe would also perform the HVG selection before I think we should make `highly_variable_genes(flavor='pearson_residuals')` as well as `normalize_pearson_residuals()` check if the...