Hello everyone,
I am using scVI to correct batch effects. However, I have encountered a problem and I need help from all of you.
My layout for my data:
adata AnnData object with n_obs × n_vars = 33486 × 352365 obs: 'features', '_scvi_batch', '_scvi_labels' var: 'orig.ident', 'nCount_RNA', 'nFeature_RNA', 'organ', 'study', 'condition', 'percent.mt'
I hope to clean the batch effect based on the study which is obs.
My code is based on this:
scvi.model.SCVI.setup_anndata(adata, layer = "counts", categorical_covariate_keys=["Sample"], continuous_covariate_keys=['pct_counts_mt', 'total_counts', 'pct_counts_ribo'])
This code do the batch effect based on obs. I am wondering how could I clean batch effect based on var.
0 answers
No answers yet.
Log in to answer this question.
-
Tutorial: Single-cell RNA-seq: Annotation: Celltype auto annotation with MetaTiME
written by Julia Ma •
-
Tutorial: Single-cell RNA-seq: Preprocessing: Data integration and batch correction-3
written by Julia Ma •
-
Tutorial: Single-cell RNA-seq: Preprocessing: Data integration and batch correction-2
written by Julia Ma •
-
Tutorial: Single-cell RNA-seq: Preprocessing: Data integration and batch correction
written by Julia Ma •
-
Tutorial: Single-cell RNA-seq: Preprocessing: Preprocessing the data of scRNA-seq with omicverse-1
written by Julia Ma •
-
Tutorial: Bulk RNA-seq: Batch correction in Bulk RNA-seq or microarray data
written by Julia Ma •
-
How to import filtered_tf_bc_matrix and add it to my anndata object?
written by bioinfo •
-
Questions Arising from Using Scanpy for Data Input and scVI for Batch Effect Removal
written by Andy •
-
from seurat to scanpy data conversion and re-scaling
written by Omics data mining
-
How to obtain all gene details from each of the clusters formed in UMAP
written by shome •
I used this common:
However, it appears that the length of my obs and var are not same. Thus this line of code helped me past the study from var to obs, but I missed some information. I am wondering if I can switch obs to var.