Create a plot showing expression profiles for individual genes across developmental stages or cell types, with various visualization options.
Arguments
- phyex_set
A PhyloExpressionSet object (BulkPhyloExpressionSet or ScPhyloExpressionSet)
- genes
Character vector of gene IDs to plot. If NULL, top expressing genes are selected
- show_set_mean
Logical indicating whether to show the mean expression across all genes (default: FALSE)
- show_reps
Logical indicating whether to show individual replicates (bulk) or cells (single-cell) (default: FALSE)
- transformation
Character string specifying expression transformation: "log" (log1p), "std_log" (standardized log1p), or "none" (default: "log")
- colour_by
Character string specifying coloring scheme: "strata" (by phylostratum), "stage" (by developmental stage/cell type), or "manual" (default: "manual")
- colours
Optional vector of colors for manual coloring (default: NULL)
- max_genes
Maximum number of genes to plot when genes=NULL (default: 100)
- show_labels
Logical indicating whether to show gene labels (default: TRUE)
- label_size
Font size of gene id labels if shown (default: 0.5).
- show_legend
Logical indicating whether to show legend (default: TRUE)
- facet_by_strata
Logical indicating whether to facet by phylostratum (default: FALSE)
Details
This function creates detailed visualizations of individual gene expression patterns across development (bulk data) or cell types (single-cell data). Genes can be colored by phylostratum or developmental stage, and various transformations can be applied to highlight different aspects of the data.
Examples
# Plot specific genes for bulk data
# p1 <- plot_gene_profiles(bulk_phyex_set, genes = c("gene1", "gene2", "gene3"))
# Plot for single-cell data with faceting by strata
# p2 <- plot_gene_profiles(sc_phyex_set, facet_by_strata = TRUE)
# p2 <- plot_gene_profiles(phyex_set, facet_by_strata = TRUE)