数据可视化——鞠艳梅-----BGI9522

最近刚学R语言~画了个简单的条形图 原数据是通过blast比对之后经过处理得到的真菌的名字，我需要 1）先统计每种真菌的频数 2）选取频数>2的画条形图 下面是我的代码 name<-read.table("c:/Users/Administrator/desktop/MGS/table") dataset<-data.frame(line1=name[,1], line2=name[,2]) a<-subset(dataset,line2>2) par(cex=0.7,mai=c(1.5,1,1,1)) bar<-barplot(a$line2,names.arg=a$line1,las=2,main="fungi_freq", xlab="fungi_name", ylab="frequency", col="blue", col.lab="red", ylim=c(0,1.2*max(a[,2]))) library(sca) ratio<-percent(a[,2]/sum(a[,2])) text(bar,a[,2],labels=ratio, cex=0.8,font=2,pos=3)