在公众号上面看到了多基因比较箱线图的绘制 结果跑出来了 但是没有结果数据 只有图。要怎么才能把结果数据也保存下来呢？
代码如下：
library(reshape2)
library(ggpubr)
inputFile="input.txt" #输入
outFile="barplot.pdf" #输出
#读取输入文件
rt=read.table(inputFile,sep="\t",header=T,check.names=F,row.names=1)
x=colnames(rt)[1]
colnames(rt)[1]="Type"
#把数据转换成gglpot2输入文件
data=melt(rt,id.vars=c("Type"))
colnames(data)=c("Type","Gene","Expression")
head(data)
#绘制boxplot
p=ggboxplot(data, x="Gene", y="Expression", color = "Type",
ylab="Gene expression",
xlab="",
legend.title=x,
palette = c("blue","red"),
width=0.6, add = "none")
p=p+rotate_x_text(60)
p1=p+stat_compare_means(aes(group=Type),
method="wilcox.test",
symnum.args=list(cutpoints = c(0, 0.001, 0.01, 0.05, 1), symbols = c("", "", "", " ")),
label = "p.signif")
#输出图片
pdf(file=outFile, width=6, height=5)
print(p1)
dev.off()