我想使用ggplot2创建一个漂亮的树形图。
这是我正在做的事情的可重现的例子:
Swift 3
这是结果
我想扩展矩形,使其覆盖x轴。如果我改变了,例如,
library(ggplot2)
library(ggdendro)
data(mtcars)
x <- as.matrix(scale(mtcars))
dd.row <- as.dendrogram(hclust(dist(t(x))))
mtcars_dendrogram <- ggdendrogram(dd.row, rotate = TRUE, theme_dendro = FALSE) +
labs(x="", y="Distance") +
ggtitle("Mtcars Dendrogram") +
theme(panel.border = element_rect(colour = "black", fill=NA, size=.5),
axis.text.x=element_text(colour="black", size = 10),
axis.text.y=element_text(colour="black", size = 10),
legend.key=element_rect(fill="white", colour="white"),
legend.position="bottom", legend.direction="horizontal",
legend.title = element_blank(),
panel.grid.major = element_line(colour = "#d3d3d3"),
panel.grid.minor = element_blank(),
panel.border = element_blank(),
panel.background = element_blank(),
plot.title = element_text(size = 14, family = "Tahoma", face = "bold"),
text=element_text(family="Tahoma"))
mtcars_dendrogram <- mtcars_dendrogram +
annotate("rect", xmin = 0.6, xmax = 5.4, ymin = 0, ymax = 6.4, fill="red", colour="red", alpha=0.1) +
annotate("rect", xmin = 5.6, xmax = 7.4, ymin = 0, ymax = 6.4, fill="blue", colour="blue", alpha=0.1) +
annotate("rect", xmin = 7.6, xmax = 11.4, ymin = 0, ymax = 6.4, fill="orange", colour="orange", alpha=0.1) +
geom_hline(yintercept = 6.4, color = "blue", size=1, linetype = "dotted")
mtcars_dendrogram
到
annotate("rect", xmin = 5.6, xmax = 7.4, ymin = 0, ymax = 6.4, fill="blue", colour="blue", alpha=0.1)
然后我明白了
这是我想要获得的(这个结果是用Photoshop改变的)
非常欢迎任何帮助。非常感谢。
答案 0 :(得分:1)
你可以伪造左轴:
mtcars_dendrogram <- mtcars_dendrogram +
annotate("rect", xmin = 0.6, xmax = 5.4, ymin = -1, ymax = 6.4, fill="red", colour="red", alpha=0.1) +
annotate("rect", xmin = 5.6, xmax = 7.4, ymin = -1, ymax = 6.4, fill="blue", colour="blue", alpha=0.1) +
annotate("rect", xmin = 7.6, xmax = 11.4, ymin = -1, ymax = 6.4, fill="orange", colour="orange", alpha=0.1) +
geom_hline(yintercept = 6.4, color = "blue", size=1, linetype = "dotted") +
theme(axis.text.y = element_blank(),
axis.line.y = element_blank(),
axis.ticks.y = element_blank()) +
geom_text(aes(y = 0, x = 1:11,
label = c("carb", "wt", "hp", "cyl", "disp", "qsec", "vs", "mpg", "drat", "am", "gear")),
hjust = "right",
nudge_y = -.1))