您的位置:首页 > 教程 > Fireworks > R语言绘制带误差线的条形图

R语言绘制带误差线的条形图

2022-06-27 11:54:52 来源:易采站长站 作者:

R语言绘制带误差线的条形图

条形统计图是用一个单位长度表示一定的数量,根据数量的多少画成长短不同的直条.带误差的条形图可以通过误差线来判断显著性。

在这里插入图片描述

继续使用我们的汽车销售数据(公众号回复:汽车销售,可以获得该数据)来演示,先导入数据

library(foreign)
library(ggplot2)
library(tidyverse)  
bc <- read.spss("E:/r/test/tree_car.sav",
                use.value.labels=F, to.data.frame=T)
names(bc)

在这里插入图片描述

我们来看下数据,car就是汽车售价,age是年龄,gender是性别,inccat是收入,这里分成4个等级,ed是教育程度。
假设我们想知道不同教育水平的男女在买汽车的价格上有什么不同,可绘制带误差和可信区间的折线图,关键就是要算出它的标准误se和95%ci.
我们先生成一个计算标准误se和95%ci的自定义函数,这是国外一位大佬设计的函数,我见好用直接搬运过来了。

summarySE <- function(data=NULL, measurevar, groupvars=NULL, na.rm=FALSE,
                      conf.interval=.95, .drop=TRUE) {
  library(plyr)
  
  # New version of length which can handle NA's: if na.rm==T, don't count them
  length2 <- function (x, na.rm=FALSE) {
    if (na.rm) sum(!is.na(x))
    else       length(x)
  }
  
  # This does the summary. For each group's data frame, return a vector with
  # N, mean, and sd
  datac <- ddply(data, groupvars, .drop=.drop,
                 .fun = function(xx, col) {
                   c(N    = length2(xx[[col]], na.rm=na.rm),
                     mean = mean   (xx[[col]], na.rm=na.rm),
                     sd   = sd     (xx[[col]], na.rm=na.rm)
                   )
                 },
                 measurevar
  )
  
  # Rename the "mean" column    
  datac <- rename(datac, c("mean" = measurevar))
  
  datac$se <- datac$sd / sqrt(datac$N)  # Calculate standard error of the mean
  
  # Confidence interval multiplier for standard error
  # Calculate t-statistic for confidence interval: 
  # e.g., if conf.interval is .95, use .975 (above/below), and use df=N-1
  ciMult <- qt(conf.interval/2 + .5, datac$N-1)
  datac$ci <- datac$se * ciMult
  
  return(datac)
}

生成函数后,我们使用自定义函数summarySE生成标准误se和95%ci. Measurevar填入你要衡量比较的指标,这里填入汽车售价,groupvars这里填入性别和教育程度。

carss<- summarySE(bc, measurevar="car", groupvars=c("gender","ed"))

生成了我们需要的做图数据

在这里插入图片描述

画条形图和画折线图不同的是,教育这个指标我们要转换成分类变量

carss$ed <- factor(carss$ed)

转换好以后就可以做图了,先做一个带误差线的

ggplot(carss, aes(x=ed, y=car, fill=gender)) + 
  geom_bar(position=position_dodge(), stat="identity") +
  geom_errorbar(aes(ymin=car-se, ymax=car+se),
                width=.2,                    # Width of the error bars
                position=position_dodge(.9))

在这里插入图片描述

画个带置信区间的

ggplot(carss, aes(x=ed, y=car, fill=gender)) + 
  geom_bar(position=position_dodge(), stat="identity") +
  geom_errorbar(aes(ymin=car-ci, ymax=car+ci),
                width=.2,                    # Width of the error bars
                position=position_dodge(.9))

在这里插入图片描述

进行美化一下,一个可以用于发表的图就做成了

在这里插入图片描述

如何利用置信区间来判断显著性

在这里插入图片描述

在这里插入图片描述

到此这篇关于R语言绘制带误差线的条形图的文章就介绍到这了,更多相关R语言 带误差线条形图内容请搜索易采站长站以前的文章或继续浏览下面的相关文章希望大家以后多多支持易采站长站!

如有侵权,请联系QQ:279390809 电话:15144810328

相关文章

  • Fireworks怎么画五子棋黑白棋盘?

    Fireworks怎么画五子棋黑白棋盘?

    用Fireworks CS5软件能够制作出漂亮的五子棋黑白棋盘,五子棋黑白棋盘的制作,主要通过网格和克隆效果以及组合路径中接合的使用来进行的。下面就以一个实例来进行说明。 1、新建一
    2019-10-16
  • Fireworks怎么将图片做成gif动图?

    Fireworks怎么将图片做成gif动图?

    相信很多人都喜欢GIF图片,搞笑,难忘。有人也想自己编辑,但是苦无软件。今天小编分享一块神器Adobe firework。 1、制作第一张小图,这个比较简单,只需设定画布的大小,在画布上写
    2019-10-16
  • Fireworks制作飞翔的海鸥GIF动画图片教程

    Fireworks制作飞翔的海鸥GIF动画图片教程

    当我们在网上冲浪时,看到形形色色的动画效果,一定会有一种冲动,那就是将这个漂亮的动画存储到自己的硬盘上,那么你知道这些动画是如何制作的的吗? 利用Firework 3制作动画总共
    2019-10-16
  • fireworks切片工具定位功能怎么使用?

    fireworks切片工具定位功能怎么使用?

    fireworks似乎慢慢地淡出了人们的视线,反而ps的地位仍然稳如泰山。但是fireworks是网页制作的三剑客之一,他在制作网页上的体现的能力还是很强的,对于切片工具的位图属性,极其便
    2019-10-16
  • Fireworks怎么为文字添加描边?

    Fireworks怎么为文字添加描边?

    Fireworks可以给文字添加描边效果,并且可以设置描边的粗细,下面我们就来看看详细的教程。 1、打开文档建好画布,然后用文字工具在右侧画板上输入文字 2、然后找到如图所示的位置
    2019-10-16
  • Fireworks绘制甲壳虫小轿车

    Fireworks绘制甲壳虫小轿车

    下面我们先拿出一张汽车原图作参考来画。 原图 效果图 用钢笔工具勾线铺色。 勾出阴影大致形状。 上一页12 下一页 阅读全文
    2019-10-16
  • Fireworks质感金属球仿真绘画教程

    Fireworks质感金属球仿真绘画教程

    最近常常看我们站点的文档评论,发现用 FireWorks 画大众甲壳虫 的评论是很多的,有很多人觉得是用照片逆向而成的假教程,当然大部分是夸赞,并且想得到详细过程和指点,我想说的
    2019-10-16
  • Fireworks快速制作带血的金属刀具教程

    Fireworks快速制作带血的金属刀具教程

    其实方法比较简单,也没有把刀做的很花哨,知道金属可以这样搞,有时间的话,谁都可以打造一件漂亮的家伙。 先来看下最终效果: 上一页12 3 4 下一页 阅读全文
    2019-10-16