MAST differential expression with batch information
Hi all,
I read that MAST can be used with batch information to properly find DE genes, however I can't quite find how this can be done. I have used MAST with Seurat before, but with no batch information. Can someone tell me how this can be done?
• 2,070 views
•
link
updated
by
ATpoint
•
written
by
Rahul Tyagi •
•
written
by
Rahul Tyagi •
0 answers
No answers yet.
Log in to answer this question.
More posts like this
-
MAST for scRNA-seq
written by learner-MD •
-
Plotting the density of CNV with karyoploteR
written by ManuelDB
-
Bad volcano plot
written by firestar
-
Sequencing file conversion
written by bill •
-
Can I use seurat to analyze the data after MNNs batch normalized?
written by jiandan123 •
-
issues with single nuclei sequencing
written by Rahul Tyagi •
-
how to deal with batch effect in TCGA RNA-seq dataset
written by tujuchuanli
-
Batch effect correction in DE analysis of single cell RNA-seq data and visualization
written by Poorya Parvizi •
-
cytoscape and network analysis
written by riyaannmathews •
-
HIGH BCV in EdgeR, any ideas?
written by Biogeek
Hello, MAST seems to take a design formula (like any other model-based framework such as DESeq2, edgeR, etc). Here is an example. This is the code run for the benchmarking paper from Soneson and Robinson (2018).
You would add the batch information as a so-called covariate to the design. This https://rdrr.io/bioc/MAST/man/zlm.html is the function that takes the design, and this you then feed e.g. into https://rdrr.io/bioc/MAST/man/lrTest.html
I am not a MAST user, maybe they have alternative tests beyond the LRT, you would need to consult the manual http://rglab.github.io/MAST/articles/MAST-Intro.html
Does that all make sense to you?