熊荣川

六盘水师范学院生物信息学实验室

xiongrongchuan@126.com

http://blog.sciencenet.cn/u/Bearjazz

R语言中通常绘制的箱体图的命令如下

boxplot( SVL ~  clade ,  data = mat.m ,  col =  colist , main = " SVL " , notch= TRUE , xaxt= 'n' )

其中SVL为数据表格中，某一列的列标题，这一列就是我们绘制箱体图（boxplot）的数据列，clade这是将数据分类的分类数据列的列标题。问题就在这里，通常我们从表格中提取列标题时，通常为字符串，简单地说就是有双引号的，所以不利于“批量”的进行箱体图的绘制，所以我们这里的思路就是批量的生成命令然后在贴回R中进行批量绘制。

setwd("D:/ziliao/优先阅读/R语言 学习/R语言boxplot箱体图命令生成表箱体图")
mat.csv = read.csv("boxplot命令表.csv")
L = length(colnames(mat)[11:37])
mat.csv[1:L,"V2"] = colnames(mat)[11:37]
mat.csv[1:L,"V4"] = "clade"
mat.csv[1:L,"V7"] = "mat.m"
mat.csv[1:L,"V10"] = "colist"
mat.csv[1:L,"V14"] = colnames(mat)[11:37]
write.table(mat.csv[1:L,],file = "boxplot命令.txt",row.names = F,col.names = F)  
mat.txt = readLines("boxplot命令.txt")
mat.txt = gsub(""","",mat.txt)
mat.txt = gsub("douhao",",",mat.txt)
mat.txt = gsub("ZuoYH",""",mat.txt)
mat.txt = gsub("YouYH",""",mat.txt)
write(mat.txt, file = "boxplot命令2.txt")
 
注：有需要“boxplot命令表.csv”的朋友可以来信索取。