单基因数据在癌症中的表达
上面得到的result.txt文件，读入R中，进行可视化。
代码如下（Rstudio）：
setwd("W:\\singleExp")
#install.packages("ggpubr")
library(ggpubr)
singleExp <- read.table("result.txt",sep="\t",
header=T,check.names=F)
colnames(singleExp)[1]="ID"
geneName <- colnames(singleExp)[2]
# 可视化
ggboxplot(singleExp,
x="tumorType",
y= geneName,
ylab=paste0(geneName," Expression Level"),
xlab="",
fill = "tissueType",
palette = c("#0072B5FF","#BC3C29FF") )+
font("caption", size = 10, color = "orange")+
font("xlab", size = 12, color = "blue")+
font("ylab", size = 14, color = "#993333")+
font("legend.title", color="black", face = "bold",size = 16)+
font("legend.text", color = "black",size = 12)+
rotate_x_text(60)+
stat_compare_means(aes(group=tissueType),
method="wilcox.test",
symnum.args=list(cutpoints = c(0, 0.001, 0.01, 0.05, 1),
symbols = c("***", "**", "*", " ")),
label = "p.signif")
得到的图片如下：
基因表达
