我们之前有分享过一系列的R语言画图代码(PCA图、Pathway图、火山图、RDA图、热图),今天再来补充一个盒形图(箱形图)的代码。
以下代码只是示例,不能直接搬来用哦,注意看注释。
--------------代码开始了-------------
setwd("E:/") #改变工作目录
data=read.table("data.txt",header=T) #读取数据,“header=T”第一行为表头
mycolors=c("#A50F15","#DE2D26","#FB6A4A","#FCAE91","#FEE5D9","white","#EFF3FF","#BDD7E7","#6BAED6","#3182BD","#08519C","#74c476",rgb(107,174,213,max=255))
matrix_data= as.matrix(data)
boxplot(matrix_data,
boxwex=0.8, #盒型宽度
main =
"ABC", #图表命名
ylim = c(10,
35), #y轴坐标范围,可以不写,可以替换成x
xlab =
"X_LAB", #x轴名称,可以替换成y
staplewex=0.6, #边缘线(订书钉)宽度
range=1.5, #边缘线扩展范围,默认1.5
border=8, #盒形图颜色,1~8,整数
col= mycolors,
#col =
"bisque", #盒子填充颜色
outline=T, #是否绘制异常值,默认为TRUE
#horizontal=TRUE, #水平,横纵坐标对换
#notch=FALSE/TRUE, #蜂腰
)
#添加标题
title("ABC")
#加图例
legend(0.5,33, #前一个是距离左侧距离,后一个是距离底端距离,现在的位置是左上角
c("a","b","c"), #图例名称
fill =
c("yellow","orange","brown") #图例填充颜色
)
#添加平均数点
points(c(1:12), #“(1:12)”要添加到1~12个盒子
c(read.table("mean.txt")), #读取平均值数据,我不会写直接运算的代码,用了最笨的办法,需要其他同事完善。
pch=15, #点的样式,0~25,图例可以在R中输入“?points”查看
cex=0.5) #点的大小
------------运行结果展示------------
当然要用起来会对代码有细微的调整,根据注释去修改一下对应的代码就可以了。PS:以上的代码缺少一些参数(例如:中位数线的粗细,颜色等)。数据分析师培训