【R>>barplot】堆叠柱状图+百分比展示 - 简书 (jianshu.com)
(125条消息) R语言的各种假设检验_hellolijunshy的博客-CSDN博客_r语言假设检验
M39B_XUO}OBDEY~~DSSNIL9.png
################一、从单细胞数据中取出信息
cell <- tumor@meta.data
input <- cell[,c(5,15)]#seurat_clusters,tissue_type
a <- data.frame(table(input$seurat_clusters,input$tissue_type))
##算百分比
library(plyr)
a<- ddply(a,.(Var1),transform,percent=Freq/sum(Freq)*100)
a$label = paste0(sprintf("%.1f", a$percent), "%")
################二、批量算p值
n<-length(unique(a$Var1))
out <- numeric()
#a$percent <- a$percent/100
for ( i in 0:(n-1)){
#i=1
shi<-a[a$Var1==i,]
p <- binom.test(round(shi[2,4]),100,p = 0.5)$p.value
out <- c(out,p)
}
out
pvalue <- out
colnames(a)[1]<-"cluster"
colnames(a)[2]<-"group"
ggplot(a,aes(cluster,percent,fill=group))+
geom_bar(stat="identity",position = position_stack())+
scale_fill_manual(values = c("#008ECA","#DB423E"),label=c("Control","Test"))+#更改ggplot填充颜色
scale_y_continuous(labels = scales::percent_format(scale = 1))+ #百分比y轴
labs(x="Tumor",y="",fill="")+
annotate("text", x ="0", y = 100.5,label = c("***"))+
annotate("text", x ="1", y = 100.5,label = c("***"))+
annotate("text", x ="2", y = 100.5,label = c("***"))+
annotate("text", x ="3", y = 100.5,label = c("***"))+
annotate("text", x ="7", y = 100.5,label = c("***"))+
annotate("text", x ="8", y = 100.5,label = c("***"))+
annotate("text", x ="9", y = 100.5,label = c("***"))+
theme_classic()+
theme(legend.position = "top",
legend.text = element_text(size=12),
axis.text = element_text(size=12),
axis.title = element_text(size=12))
