matlab 画柱状图和误差棒,小白1分钟绘图:带误差棒的柱状图-程序员宅基地

技术标签: matlab 画柱状图和误差棒  

目标图片效果

e81506c31aac

想要实现的图片效果

其实上图说白了就是均值加减标准差,上面的线条是用来标示组间比较的。

制图清单

1.工具: R+Rstudio

2.数据准备:需要如下样式的数据。

第一列是分组,第二列是均值,第三列是标准差。

e81506c31aac

数据格式

3.设置R软件工作目录:

为了方便找到结果的图,可以把工作目录设置为桌面。操作如下:

e81506c31aac

工作目录设置

4.跑代码:打开Rstudio导入数据运行代码就行了~

将数据文件放在桌面,导入数据,注意数据格式为csv,在括号中按Tab键进行导入。

rm(list=ls())#清除掉全局环境中的文件

data

加载需要用到的R软件包

##########加载需要用的R包

if (!require(ggplot2)) {install.packages("ggplot2")}

if (!require(ggsignif)) {install.packages("ggsignif")}

if (!require(ggsci)) {install.packages("ggsci")}

if (!require(export)) {install.packages("export")}

加载数据,构建绘图函数。(直接运行就行,看不懂也没关系。)

#######载入数据

attach(data)

#构建作图函数

myplot

ggplot(data,aes(eval(parse(text=names(data)[1])),eval(parse(text=names(data)[2])),

fill=factor(eval(parse(text=names(data)[2])))))+

geom_bar(position = "dodge",stat = "identity")+

geom_errorbar(aes(ymin=eval(parse(text=names(data)[2]))-eval(parse(text=names(data)[3])),

ymax=eval(parse(text=names(data)[2]))+eval(parse(text=names(data)[3])),

color=factor(eval(parse(text=names(data)[2])))),position = "dodge",width=0.3,size=2)+

scale_y_continuous(expand = c(0,0))+

labs(x=names(data)[1],y=names(data)[2])+

theme_classic(base_size = 22)+

theme(axis.text.x = element_text(angle = -45,hjust = 0.2),

legend.position = "none")+

scale_fill_lancet()+ #####这里用了lancet杂志的配色方案

scale_color_lancet()

}

出图~

###绘图

myplot(data)

e81506c31aac

本文示例数据出的图

修改Y轴长度,不需要修改这一步可以跳过。

###将画出来的图赋值方便修改

p

##如果你想改Y轴的范围,可以修改下面“limits=”中的数值

p+scale_y_continuous(limits = c(0,30),expand = c(0,0))+

scale_fill_lancet()+

scale_color_lancet()

接下来就是导出到ppt里,在ppt里就可以直接加上组间对比的线条。

导出的图所几乎所有对象都是可以修改的,可以自己DIY标题,颜色等等~

###导出成ppt自己修改细节

graph2ppt(file="myplot.pptx",width=7,height=5)

e81506c31aac

ppt中的图

扩展

当然,不用PPT也行,R本身就足够强大,但是需要多了解一点图形语法。

示例图中线条的增加,可以用ggsignif包来实现,当然用ggplot2自带的annotate函数也可以。

###根据本图实际情况的示例代码

p+scale_y_continuous(limits = c(0,40),expand = c(0,0))+

annotate("errorbarh",xmin=1,xmax = 2,y=25,size=1)+

annotate("errorbarh",xmin=1,xmax = 6,y=36,size=1)+

annotate("errorbarh",xmin = 3,xmax=4,y=29,size=1)+

annotate("errorbarh",xmin = 1,xmax=5,y=32,size=1)+

annotate("text",label="*",x=1.5,y=26,fontface="bold",size=5)+

annotate("text",label="*",x=3.5,y=33,size=5,fontface="bold")+

annotate("text",label="ns",x=3.5,y=37.5,size=5,fontface="bold")+

scale_fill_lancet()+

scale_color_lancet()

跑代码出图。

e81506c31aac

效果图

补充一个实际使用中更常用的方法,借助ggpubr包

library(ggpubr)

#设定比较组

my_comparisons

1:nrow(t(combn(levels(as.factor(data$group)),2))))

#添加对比p值

p+ stat_compare_means(comparisons = my_comparisons, label.y = c(29, 35, 40))+

stat_compare_means(label.y = 45) # Add 全局nova p-value

p值保留小数点问题

#以保留两位小数点为例

p + stat_compare_means(

aes(label = paste0("p = ", round(as.numeric(..p.format..,2))))

)

小结

带误差棒柱状图算是实验类论文中很常用的图形。

实现起来也不难,本文绘图的主函数,在笔者五年前的破烂小电脑里运行出图,耗时如下,快得飞起。

e81506c31aac

绘图耗时计算

加上载入数据,运行其他代码,还有喝口水的时间,大概需要一分钟吧,可以尝试拷贝代码在自己的电脑运行起来吼~(─‿‿─)

版权声明:本文为博主原创文章,遵循 CC 4.0 BY-SA 版权协议,转载请附上原文出处链接和本声明。
本文链接:https://blog.csdn.net/weixin_39759600/article/details/115850219

智能推荐

java获取body_java获取http请求的Header和Body的简单方法-程序员宅基地

文章浏览阅读3.6k次。在http请求中,有Header和Body之分,读取header使用request.getHeader("...");读取Body使用request.getReader(),但getReader获取的是BufferedReader,需要把它转换成字符串,下面是转换的方法。public class TestController {@RequestMapping("/a")protected void..._java获取body的值

mysql中union all使用注意点_mysql union all 中使用 含order by子查询 注意事项-程序员宅基地

文章浏览阅读320次。今天在J8看到一道题:表A里有id money两个字段 当money大于100的时候倒序排 小于100的时候按照id排 写出查询语句baidu一下,可以用union查询。找了自己的一张表试了下(select * from (select * from `t_user` where head>3 order by autoid desc) t1)union(select * from (sel..._mysql union all 中使用 oder by

java sheet_JAVA读取Excel中sheet页内容-程序员宅基地

文章浏览阅读1.2k次。packagecom.bestpay.bpbp.bill.init;importjxl.Cell;importjxl.Sheet;importjxl.Workbook;importjxl.read.biff.BiffException;importjava.io.File;importjava.io.IOException;publicclassReadExcelDataForT..._java 读取excel sheet页为表名 sheet页中的字段为字段值 类型 和 注释 ,根据以上

R语言安装glmnet包_glmnet包安装-程序员宅基地

文章浏览阅读8.4k次,点赞9次,收藏8次。R语言安装glmnet包安装glmnet包安装glmnet包正常这样安装时报错install.packages("glmnet")改为如下代码可正常安装。install.packages("glmnet", dependencies=TRUE)具体原因等写完作业有时间再查。_glmnet包安装

mk突变点检测_气候突变检测-程序员宅基地

文章浏览阅读3k次。1、气候突变是什么意思科学家们所说的气候变化,通常指的是“气候渐变”。也就是我们经常说的,“地球在逐渐变暖”,或者是“气候在慢慢变暖”。但现在已有证据显示地球上的气候系统有时更像一个会在瞬间变化的“开关”,即“突变”,而不是渐变的“指针盘”。2、小气候的变化监测在森林生态系统中,进行裸岩、草地、灌丛三种目前龙何分布最广的植被类型的进行对比性小气候监测,进行不知同的立体农业覆盖5种模式进行小气候对比..._mk的dna检测

mysql查看relaylog参数_MySQL relay log 详细参数解释-程序员宅基地

文章浏览阅读1.6k次。前言:MySQL进行主主复制或主从复制的时候会在home目录下面产生相应的relay log,本文档总结这些相关参数的定义及解释.1、什么是relay logThe relay log, like the binary log, consists of a set of numbered files containing events that describe database changes,..._mysql 查看relay log

随便推点

openstack虚拟机创建过程中生成的xml文件位置_如何查openstack创建的虚拟机的xml文件-程序员宅基地

文章浏览阅读3.4k次。mktaka版本的openstack在启动该虚拟机的nova节点上,在/var/lib/nova/instances/xxx目录下有四个文件,console.log disk disk.info libvirt.xml其中libvirt.xml文件里写的就是虚拟机启动信息文件_如何查openstack创建的虚拟机的xml文件

使用Opengl绘制字体_opengl写字的函数-程序员宅基地

文章浏览阅读3.1k次。使用Opengl绘制字体通过使用windos自带的字体,可以在opengl中方便的绘制多种字体。bitmapFont类下面介绍 bitmapFont 字体类及其成员函数:首先是构造函数。第一个构造函数:与winapi函数CreateFont(...)函数参数基本相同。bitmapFont( HDC &hdc,_opengl写字的函数

oracle rac vip不通,Oracle RAC 中vip网关引起错误的解决-程序员宅基地

文章浏览阅读794次。在虚拟机上刚装完Oracle RAC,重新启动系统后。重新启动了oracle所有的服务,一开始所有的服务都正常。但是过了一会,节点2的vip又再次飘移到了节点1上,并且节点2的监听服务也OFFLINE了。反复启动了几次都是同样的效果的。[root@node2 bin]# ./crs_stat -tName Type Target State H..._oracle rac db.vip failed

【DW组队学习—动手学数据分析】第一章:第三节探索性数据分析-课程学习-程序员宅基地

文章浏览阅读239次。【DW组队学习—动手学数据分析】第一章:第三节探索性数据分析-课程学习

ios 系统状态栏样式修改_iOS设置状态栏样式-程序员宅基地

文章浏览阅读341次。标签:状态栏样式,statusBarStyle状态栏的样式系统的默认样式是黑色的,但可以修改,方法有2种方法1:// 返回状态栏的样式- (UIStatusBarStyle)preferredStatusBarStyle{return UIStatusBarStyleLightContent;}// 控制状态栏的现实与隐藏- (BOOL)prefersStatusBarHidden{return..._ios 设置当前页面 状态栏样式

Verilog基础入门-程序员宅基地

文章浏览阅读2.7w次,点赞109次,收藏705次。Verilog基础入门_verilog