Get coefficient names
Examples
library(muscat)
library(SingleCellExperiment)
data(example_sce)
# create pseudobulk for each sample and cell cluster
pb <- aggregateToPseudoBulk(example_sce,
assay = "counts",
cluster_id = "cluster_id",
sample_id = "sample_id",
verbose = FALSE
)
# voom-style normalization
res.proc <- processAssays(pb, ~group_id)
#> B cells...
#> 0.22 secs
#> CD14+ Monocytes...
#> 0.32 secs
#> CD4 T cells...
#> 0.28 secs
#> CD8 T cells...
#> 0.14 secs
#> FCGR3A+ Monocytes...
#> 0.27 secs
# Differential expression analysis within each assay,
# evaluated on the voom normalized data
res.dl <- dreamlet(res.proc, ~group_id)
#> B cells...
#> 0.2 secs
#> CD14+ Monocytes...
#> 0.27 secs
#> CD4 T cells...
#> 0.22 secs
#> CD8 T cells...
#> 0.14 secs
#> FCGR3A+ Monocytes...
#> 0.27 secs
# show coefficients estimated for each cell type
coefNames(res.dl)
#> [1] "(Intercept)" "group_idstim"