align_gg
align_gg() is similar to ggplot in that it
initializes a ggplot data and mapping. Same
with other align_* functions. align_gg()
allowing you to provide data in various formats, including matrices,
data frames, or simple vectors. This data can be also inherited from the
layout.
align_gg() always applies a default mapping for the axis
of the data index in the layout. This mapping is
aes(y = .data$.y) for horizontal stacking (including left
and right heatmap annotation) and aes(x = .data$.x) for
vertical stacking (including top and bottom heatmap annotation). For
more information, refer to the “ggplot2 Specification” section in the
align_gg() documentation.
You can also use the ggalign() function, which is an
alias for align_gg().
ggheatmap(small_mat) +
scale_fill_viridis_c(guide = "none") +
hmanno("t") +
ggalign(data = rowSums) +
geom_point(aes(y = value))
Plot data
align_gg()/ggalign() requires the specific
data format for its operations. If you need to transform or filter data
for individual geoms, you can use the data
argument within each geom. However, if you have multiple
geoms and want a consistent transformation applied across
all, you can utilize the plot_data argument in the
align_gg()/ggalign() function. This allows you
to transform the default data for all subsequent geoms.
ggheatmap(small_mat) +
scale_fill_viridis_c(guide = "none") +
hmanno("t") +
align_kmeans(3L) +
ggalign(plot_data = function(data) subset(data, .panel == 1L)) +
geom_bar(aes(y = value, fill = .row_names), stat = "identity")
Cross panel sumamry
When used in a heatmap layout, and the data is inherited from the
heatmap data, a special column .extra_panel will be added,
which is the panel information for column (left or right annotation) or
row (top or bottom annotation). This is useful if you want to create
summary plot using another axis panel. In such cases, it’s often
necessary to disable the automatic setting of limits
(limits = FALSE in align_gg()).
set.seed(1L)
v <- stats::rnorm(50L)
split <- sample(letters[1:2], 50L, replace = TRUE)
ggheatmap(v) +
scale_fill_viridis_c() +
theme(strip.text = element_text(), strip.background = element_rect()) +
hmanno("r") +
align_group(split) +
hmanno("t", size = 0.5) +
ggalign(limits = FALSE) +
geom_boxplot(aes(.extra_panel, value, fill = .extra_panel),
# here, we use `print()` to show the underlying data
data = function(data) {
print(head(data))
data
}
) +
scale_fill_brewer(palette = "Dark2", name = "branch")
#> .column_index .row_index value .extra_panel .panel .x
#> 1 1 1 -0.6264538 b 1 1
#> 2 2 1 0.1836433 b 1 1
#> 3 3 1 -0.8356286 b 1 1
#> 4 4 1 1.5952808 a 1 1
#> 5 5 1 0.3295078 a 1 1
#> 6 6 1 -0.8204684 a 1 1This approach replicates the functionality of ComplexHeatmap::anno_summary(), but is versatile enough to be used with any heatmap, not just single-column or single-row heatmaps.
ggheatmap(small_mat) +
theme(axis.text.x = element_text(angle = -60, hjust = 0)) +
hmanno("t") +
align_dendro(aes(color = branch), k = 3L) +
scale_color_brewer(palette = "Dark2") +
hmanno("r", size = 0.5) +
ggalign(limits = FALSE) +
geom_boxplot(aes(y = .extra_panel, x = value, fill = factor(.extra_panel))) +
scale_fill_brewer(palette = "Dark2", name = "branch")
Plot titles
ggplot2 only allow add titles in the top or add caption
in the bottom. The ggalign package extends this capability, allowing you
to place titles around any border of the plot using the
patch_titles() function.
ggheatmap(small_mat) +
patch_titles(left = "left patch title", bottom = "bottom patch title") +
theme(axis.text.x = element_text(angle = -60, hjust = 0)) +
hmanno("t") +
align_dendro(aes(color = branch), k = 3L) +
scale_color_brewer(palette = "Dark2") +
patch_titles(top = "top patch title") +
hmanno("r", size = 0.5) +
ggalign(limits = FALSE) +
geom_boxplot(aes(y = .extra_panel, x = value, fill = factor(.extra_panel))) +
scale_fill_brewer(palette = "Dark2", name = "branch") +
patch_titles(right = "right patch title")
