Can I integrate data with iCellR? about icellr

[rezakj]

Yes, it is very easy. First read each patient data like this:

P1 <- load10x("filtered_gene_bc_matrices/hg19/")
P2 <- load10x("filtered_gene_bc_matrices/hg19/")
P3 <- load10x("filtered_gene_bc_matrices/hg19/")
...
P50 <- load10x("filtered_gene_bc_matrices/hg19/")

Then use the data.aggregation function to make them into one file while also labeling each patient/condition.

Example for 3 patients:
my.data <- data.aggregation(samples = c("P1", "p2","P3"), condition.names = c("P1", "p2","P3"))

Then make an object:

my.obj <- make.obj(my.data)
my.obj

CPCA is for batch alignment, which you run instead of PCA step if you wish to do batch alignment.

[Midsummer723]

Thank you for the quick reply. The CPCA can work for the batch effect,right? If it is, how can I use it to my data . Like the "my.obj <- iba(my.obj,dims = 1:30, k = 10,ba.method = "CPCA", method = "gene.model", gene.list = [email protected])" change the gene.model to orig.ident (sorry , I use the seurat)? OR if not, how can I integrate the data to aviod the batch effect.

[rezakj]

Yes, you can replace [email protected] with orig.ident or any gene list. It is technically just a number of genes you want to use for this. It could be a list of top dispersed genes or top 1000-2000 expressed genes or combination of both or a set of genes you are interested in (eg. c(gene1,gene2,...)).

[Midsummer723]

Thank you very much!