R/plot_marker_express_ALL.R
plot_marker_express_ALL.Rd
This function plots the top n marker genes for a all cell types based off of
the stats
table from get_mean_ratio()
in a multi-page PDF file.
The gene expression is plotted as violin plot with plot_gene_express()
and
adds annotations to each plot.
plot_marker_express_ALL(
sce,
stats,
pdf_fn = NULL,
n_genes = 10,
rank_col = "MeanRatio.rank",
anno_col = "MeanRatio.anno",
gene_col = "gene",
cellType_col = "cellType",
color_pal = NULL,
plot_points = FALSE
)
SummarizedExperiment-class object
A data.frame()
generated by get_mean_ratio()
and/or
findMarkers_1vAll()
.
A character()
of the pdf filename to plot to, if NULL
returns all plots
An integer(1)
of number of markers you'd like to plot
The character(1)
name of column to rank genes by in
stats
.
The character(1)
name of column containing annotation in
stats
.
The character(1)
name of column containing gene name in
stats
should be the same syntax as rownames(sce)
.
The character(1)
name of colData()
column containing
cell type for sce
data. It should match cellType.target
in stats
.
A named character(1)
vector that contains a color pallet
matching the cell_type
values.
A logical(1)
indicating whether to plot points over the
violin, defaults to FALSE
as these often become over plotted and quite
large (especially when saved as PDF).
A PDF file with violin plots for the expression of top marker genes for all cell types.
Other expression plotting functions:
plot_gene_express()
,
plot_marker_express()
,
plot_marker_express_List()
#' ## Fetch sce example data
if (!exists("sce_DLPFC_example")) sce_DLPFC_example <- fetch_deconvo_data("sce_DLPFC_example")
#> 2024-10-31 13:27:17.769337 loading file /github/home/.cache/R/BiocFileCache/66141272054e_sce_DLPFC_example.Rdata%3Frlkey%3Dv3z4u8ru0d2y12zgdl1az07q9%26st%3D1dcfqc1i%26dl%3D1
## load example marker stats
data("marker_test")
# Plot marker gene expression to PDF, one page per cell type in stats
pdf_file <- tempfile("test_marker_expression_ALL", fileext = ".pdf")
plot_marker_express_ALL(
sce_DLPFC_example,
cellType_col = "cellType_broad_hc",
stat = marker_test,
pdf_fn = pdf_file
)
#> [[1]]
#>
#> [[2]]
#>
#> [[3]]
#>
#> [[4]]
#>
#> [[5]]
#>
#> [[6]]
#>
#> [[7]]
#>
#> pdf
#> 2
if (interactive()) browseURL(pdf_file)