I have the following data:
df1 <- structure(list(Wastewater_Treatment_Plant = c("D_S1_L001_R1_001",
"D_S1_L001_R1_001", "D_S1_L001_R1_001", "D_S1_L001_R1_001", "D_S1_L001_R1_001",
"D_S1_L001_R1_001", "D_S1_L001_R1_001", "D_S1_L001_R1_001", "D_S1_L001_R1_001",
"D_S1_L001_R1_001", "D_S1_L001_R1_001", "D_S1_L001_R1_001", "D_S1_L001_R1_001",
"D_S1_L001_R1_001", "D_S1_L001_R1_001", "D_S1_L001_R1_001", "D_S1_L001_R1_001",
"D_S1_L001_R1_001"), Domain = c("Archaea", "Archaea", "Archaea",
"Archaea", "Archaea", "Bacteria", "Bacteria", "Bacteria", "Bacteria",
"Bacteria", "Bacteria", "Bacteria", "Eukaryota", "Eukaryota",
"Eukaryota", "Eukaryota", "Other Sequences", "Viruses"), Phylum = c("Crenarchaeota",
"Euryarchaeota", "Korarchaeota", "Nanoarchaeota", "Thaumarchaeota",
"Acidobacteria", "Actinobacteria", "Aquificae", "Bacteroidetes",
"Candidatus Poribacteria", "Chlamydiae", "Chlorobi", "Streptophyta",
"Xanthophyceae", "unclassified (derived from Eukaryota)", "unclassified (derived from Fungi)",
"unclassified (derived from other sequences)", "unclassified (derived from Viruses)"
), Alignment_Length = c(3573, 34060, 257, 22, 525, 85973, 670251,
8825, 1040376, 1273, 5962, 41954, 13026, 5, 15129, 48, 528, 3451
), Relative_Alignment_Length = c(0.00185587444253646, 0.0176913192031323,
0.000133489989289636, 1.14271586162334e-05, 0.000272693557887388,
0.044655777623338, 0.348139294985867, 0.00458384885401182, 0.540388253296476,
0.000661216950839325, 0.00309675998499926, 0.0217915914811571,
0.00676591673341166, 2.59708150368941e-06, 0.00785824921386343,
2.49319824354184e-05, 0.000274251806789602, 0.00179250565384643
), ymax = c(0.00185587444253646, 0.0195471936456687, 0.0196806836349584,
0.0196921107935746, 0.019964804351462, 0.0646205819748, 0.412759876960667,
0.417343725814679, 0.957731979111154, 0.958393196061993, 0.961489956046993,
0.98328154752815, 0.990047464261561, 0.990050061343065, 0.997908310556929,
0.997933242539364, 0.998207494346154, 1), ymin = c(0, 0.00185587444253646,
0.0195471936456687, 0.0196806836349584, 0.0196921107935746, 0.019964804351462,
0.0646205819748, 0.412759876960667, 0.417343725814679, 0.957731979111154,
0.958393196061993, 0.961489956046993, 0.98328154752815, 0.990047464261561,
0.990050061343065, 0.997908310556929, 0.997933242539364, 0.998207494346154
)), row.names = c(NA, -18L), class = c("tbl_df", "tbl", "data.frame"
))
Domain_Colors <- c("orange", "blue", "grey", "green", "purple")
Phyla_Colors <- c("#FEEDDE", "#FDBE85", "#FD8D3C", "#E6550D", "#A63603", "#440154FF",
"#443A83FF", "#31688EFF", "#21908CFF", "#35B779FF", "#8FD744FF",
"#FDE725FF", "#4D4D4D", "#969696", "#C3C3C3", "#E6E6E6", "green1",
"purple1")
I've created the following chart:
if (!require (dplyr)) {
install.packages("dplyr")
}
library (dplyr)
if (!require (ggnewscale)) {
install.packages("ggnewscale")
}
library (ggnewscale)
if (!require (ggplot2)) {
install.packages("ggplot2")
}
library (ggplot2)
if (!require (tidyr)) {
install.packages("tidyr")
}
library (tidyr)
ggplot(df1) +
geom_rect(aes(fill = Domain, ymax = ymax, ymin = ymin, xmax = 2, xmin = 0)) +
scale_fill_manual(aesthetics = "fill", values = Domain_Colors, breaks = unique(df1$Domain), name = "Domain") +
new_scale_fill() +
geom_rect(aes(fill = Phylum, ymax = ymax, ymin = ymin, xmax = 4, xmin = 2)) +
scale_fill_manual(aesthetics = "fill", values = Phyla_Colors, breaks = unique(df1$Phylum), name = "Phylum") +
coord_polar(theta = "y") +
theme_void()
Is there a way to rearrange the legend such that the phyla that belong to a particular domain are listed under each domain? In other words, I'd like to put the domain (with its associated color) at the head of each group of phyla that belong to it.
Additionally, it would be nice to keep the legend titles 'Domain' and 'Phylum' in the legend, corresponding to the different domains and phyla.
When all is said and done, the legend should look like this:
Thank you!