蜜蜂群图绘制

输出结果

绘图命令

1	bee_plot_with_color("../../Desktop/seqencing depth.csv","xlab","ylab")

函数代码

bee_plot_with_color=function(In_file,xlab_text,ylab_text){
library(beeswarm)
read.csv(In_file,header=F)->a
#biocLite(c("beeswarm","ggplot2"))

label_tag=a[1,]
label_tag=as.matrix(label_tag)
a=a[-1,]
ncol(a)->num1
mat=c()
for(i in 1:num1){
as.numeric(as.matrix(a[,i]))->temp
temp[!is.na(temp)]->temp
length(temp)->num
rep(i,num)->label_temp
rbind(mat,cbind(temp,label_temp))->mat

}
beeswarm(mat[,1]~mat[,2], pch = 1,col = rainbow(10),labels=label_tag ,xlab=xlab_text , ylab=ylab_text)
#pch对应的不同数值可以调整散点图中点的样式；
boxplot(mat[,1]~mat[,2],add=T,names=label_tag)->S
}

输入文件：

每列对应一类数据，head行，对应类的名称，下面每行是该类的具体数据，最终针对每类数据绘制对应的散点图和4分位图。

demo数据