1.代码

args = commandArgs(T)
if (length(args) !=4){
    print("Rscript this R <infile> <element1> <element2> <element3> <outfile>")
    q()
}

options(stringsAsFactors=F)
library(reshape2)
library(ggplot2)

data<-read.table(args[1],header=T,sep="\t")
data1<-data[order(data[args[3]]),]
id<- c(1:dim(data1[args[2]])[1])
data2<-cbind(id,data1)
data3<-data2[,-2]
long_data <- melt(data3, id="id")
pdf(args[4],height=4,width=16)

ggplot(long_data,aes(id, value)) + geom_point( alpha=0.25, shape=16,size = 1,stat = "identity",position = "identity") + xlab(args[2]) + ylab('Shannon') +geom_line(aes(id, value,group=variable,color=variable),alpha=0.55)+scale_x_continuous( breaks=c(1:dim(data1[args[2]])[1]) ,labels=as.vector(unlist(data1[args[2]])) , expand = c(0,0))+theme(axis.text.x = element_text( color="black", size=3.5, angle=45,hjust=1),axis.ticks = element_blank(),panel.grid.major.x=element_blank(), panel.grid.minor.x=element_blank())
dev.off()

2.使用示例

image.png

Rscript plot.r danyangbenbijiao.xls Sample H.sample_P 单个样本熵值比较.pdf

3.上游数据示例

image.png

4.图

image.png