library(ggpubr)
setwd("C:\\Users\\ALIENWARE\\")            
rt=read.table("risk.txt",sep="\t",header=T,row.names=1,check.names=F)        
rt=rt[,3:(ncol(rt)-2)]
Type=read.table("clinical.txt",sep="\t",check.names=F,header=T)
row.names(Type)=Type[,1]
sameSample=intersect(row.names(Type),row.names(rt))
rt=rt[sameSample,]
Type=Type[sameSample,]
#准备
data=data.frame()
for(i in colnames(rt)){
  data=rbind(data,cbind(expression=rt[,i],gene=i,Type=as.vector(Type[,2])))
}
write.table(data,file="data.txt",sep="\t",row.names=F,quote=F)
#作图
data=read.table("data.txt",sep="\t",header=T,check.names=F)      
p=ggboxplot(data, x="gene", y="expression", color = "Type",
     ylab="lncRNA expression",
     xlab="",
     palette = rainbow(length(levels(data$Type))) )
p=p+rotate_x_text(45)
pdf(file="boxplot.pdf",width=7,height=5)                          
p+stat_compare_means(aes(group=Type),symnum.args=list(cutpoints = c(0, 0.001, 0.01, 0.05, 1), symbols = c("***", "**", "*", "ns")),label = "p.signif")
dev.off()