0
我有ggplot热图,集群中有kmeans(15个集群)。在ggplot热图中添加集群的图例
MotifPresence.kmeans = kmeans(MotifPresence, 15)
class.promoter = MotifPresence.kmeans$cluster
MotifPresence.kmeans.ordered = MotifPresence[order(class.promoter), MPbindist.hclust$order]
x=colnames(MotifPresence)[MPbindist.hclust$order]
y=rownames(MotifPresence)
data <- expand.grid(X=x, Y=y)
data$Z <- as.vector(t((MotifPresence.kmeans.ordered)))
ggplot(data, aes(X, Y, z= Z)) + geom_tile(aes(fill = Z)) +
theme(axis.text.x = element_text(angle = 45, hjust = 1),
axis.text.y = element_blank(), axis.ticks.y=element_blank())+
scale_fill_gradient(name = "Score", low = "white", high = "darkgreen") +
ylab("") +xlab("")
如何在此热图上添加一些标记以显示这些集群之间的边界?大概水平线或每个群集有不同颜色的附加列?