Exercise

Fit zero-inflated regression using MAST

To identify genes that are differentially expressed (DE) between cell type 9 and the other cell types, we want to fit a zero-inflated regression using the function zlm() from the package MAST.

Instructions

100 XP

Load MAST.
Print sca, which was created from sce.
Fit a zero-inflated regression to sca regressing variables celltype and cngeneson to adjust for the cell detection rate. Assign it to zir_model.
Summarize model features from the object zlm. Assign it to the variable summary_lrt, calculating a likelihood ratio test (LRT) on each gene for celltype9 using the argument doLRT.