r语言boxcox异方差_科学网—R语言boxplot箱体图命令生成表 - 熊荣川的博文

熊荣川

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

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”的朋友可以来信索取。

转载本文请联系原作者获取授权，同时请注明本文来自熊荣川科学网博客。

链接地址：http://blog.sciencenet.cn/blog-508298-756449.html

上一篇：Sequencher软件进行极速DNA序列拼接

下一篇：R语言终止密码子统计