单细胞RNA测序

ggplot点图geom_point并添加分组标注

2020-03-28  本文已影响0人  yingyonghui

自编代码实现Seurat TSNEPlot函数功能

p <- ggplot() + geom_point(data=tsne.cor.filtered, aes(x=tSNE_1,y=tSNE_2,colour=Cluster),size =0.5) +
  scale_colour_manual(values=mycolors[1:length(levels(tsne.cor.filtered$Cluster))])  + 
  guides(color = guide_legend(override.aes = list(size = 3))) +
  theme(legend.title=element_blank(),legend.key = element_blank(),legend.text=element_text(colour= 'black',size=14),
  axis.text= element_text(colour= 'black',size=14),axis.title= element_text(size = 14),axis.line= element_line(colour= 'black'),
  panel.background=element_rect(fill="white"))

#计算标签所在的位置
centers <- tsne.cor.filtered[,c('tSNE_1','tSNE_2','Cluster')] %>%
  dplyr::group_by(Cluster) %>%
  summarize(x=median(x=tSNE_1), y=median(x=tSNE_2))

p <- p  + geom_text(data=centers,aes(x,y,label=Cluster))
p

#将同一组的散点圈起来
ggplot(data=intestine.plot.meta.dat) + 
    geom_point(aes(x=tSNE_1,y=tSNE_2,colour=Clusters),size =0.5) +
    scale_colour_manual(values=mycolors[1:54])  +
    stat_ellipse(aes(x=tSNE_1,y=tSNE_2,fill=Clusters),geom="path", level=0.95, alpha=0.2)+
    guides(color=guide_legend('Clusters',override.aes=list(size=3),ncol=2),fill=F)
上一篇 下一篇

猜你喜欢

热点阅读