R语言的画图代码及差异性分析[通俗易懂]

方差分析及可视化

大家好,又见面了,我是你们的朋友全栈君。

R语言结合统计学方法画图实现

例如:


提示:写完文章后,目录可以自动生成,如何生成可参考右边的帮助文档


前言

加载我们常用的包,如ggplot2 , 在这里我只使用R的可视化功能,前期的数据清洗及处理,将数据转化为R可识别的数据格式, 我均是通过Python的pandas包来实现的,强烈推荐使用python清洗数据!!!! 自动化加星号!!!!(显著性!!!!)


提示:以下是本篇文章正文内容,下面案例可供参考

一、ggplot2是什么?

示例:ggplot2 是R语言可视化的一个包,可视化功能非常强大。

二、使用步骤

1.加载我们数据分析及画图所需要的的包

代码如下(示例):

library(xlsx)
library("Hmisc") ## to use "rcorr"
library(psych)
library(lavaan)
library(semTools)
library(haven)
library(ggcorrplot)
library(semPlot)
library(foreign)
library(MASS)
library(tidyr)
library(jtools)
library(dplyr)
library(cluster)
library(fpc)
library(ggplot2)
library(dendextend)
library(purrr)
library(ggpubr) # to use ggplotbar

2.读入数据

代码如下(示例):

setwd("F:/课题/数据/处理好的数据/20/")
raw_dt = read.xlsx2("*****.xlsx", sheetIndex = 1)

画图代码

箱线图

ggplot(dt1,aes(x=name,y=speed,fill=factor(name))) +
  stat_boxplot(geom = 'errorbar')+
  geom_boxplot() +
  stat_compare_means(label.y=5,method = 'anova') +
  labs(title = "你的标题")+
  xlab("字符串")
##代码2
boxplot(speed~group_name, dt1,  ylab = "speed, m/s",xlab= "A=60,B=40,C=20,单位:min",col = "bisque" )



增加显著性标志,加星星

ggplot(dt2,aes(x=Usedtime,y=Trust,fill=factor(Usedtime))) +
  geom_boxplot() +
  stat_boxplot(geom = 'errorbar')+
  stat_compare_means(label.y=5.8,method = 'anova')+
  stat_compare_means(comparisons = list(c(1,2),c(1,3),c(1,4),c(2,3),c(2,4),c(3,4)),label.y=c(4.9,5.2,5.5,5.7,5.9,6),label ='p.signif')+
#list(c(1,2),c(1,3)),label.y=c(4.9,5.2,5.5,5.7,5.9,6),label ='p.signif',hide.ns=T)+
  labs(title = "特斯拉车主")+
  xlab("使用时间")

双因素方差分析

双因素方差分析指的是,我们对于想要研究的某个变量,有两个变量会影响。
我这里研究的是驾驶模式(手动/自动)和驾驶时间(长/中/短)对于驾驶绩效的影响。
双因素分析代码

dt1 = read.csv("C:/Users/Heart/Desktop/接管绩效汇总_接管成功_control.csv")
attach(dt1)
aggregate(acc_x,by=list(group_name,control),FUN=mean,na.rm=TRUE)
aggregate(acc_x,by=list(group_name,control),FUN=sd,na.rm=TRUE)
fit_acc_x <- aov(acc_x~group_name*control)
summary(fit_acc_x)

总结

提示:这里对文章进行总结:

版权声明:本文内容由互联网用户自发贡献,该文观点仅代表作者本人。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如发现本站有涉嫌侵权/违法违规的内容, 请联系我们举报,一经查实,本站将立刻删除。

发布者:全栈程序员-站长,转载请注明出处:https://javaforall.net/128465.html原文链接:https://javaforall.net

(0)
全栈程序员-站长的头像全栈程序员-站长


相关推荐

  • navaicat 15 for mysql 激活码【在线破解激活】[通俗易懂]

    navaicat 15 for mysql 激活码【在线破解激活】,https://javaforall.net/100143.html。详细ieda激活码不妨到全栈程序员必看教程网一起来了解一下吧!

    2022年3月16日
    48
  • eplan win10 能用激活码激活码_通用破解码

    eplan win10 能用激活码激活码_通用破解码,https://javaforall.net/100143.html。详细ieda激活码不妨到全栈程序员必看教程网一起来了解一下吧!

    2022年3月17日
    592
  • Origin绘图快速上手指南「建议收藏」

    Origin绘图快速上手指南「建议收藏」1、创建工程打开origin后,点击菜单栏“文件”,选择“项目另存为”,给项目命名,并存到某个工作路径。2、导入数据然后将excel中的数据(只要数据)选中后复制到Book1中,从第5行开始粘贴。可以在侧面打开“项目管理器”,给表格“Book1”重命名为“曲线数据”。还可以在表格的“长单位”处给每列数据加上标签。3、那么这时可以直接使用Origin的自动绘图功能了。选择A、B、C所有列,然后点击菜单栏的“绘图”,选择一个折线图,双击即可绘图。这样呢就是将两条曲线放到同一张图中了。如果想要自定

    2022年5月31日
    57
  • 水流量霍尔传感器工作原理_51单片机温度传感器

    水流量霍尔传感器工作原理_51单片机温度传感器硬件准备(1)YF-B1流量传感器一个(2)51开发板一个基础知识(1)YF-B1流量传感器只有三根线。即数据线、VCC、GND。数据线输出为占空比为50%的方波。当水流通过水流转子组件时,磁性转子转动并且转速随着流量的变化而成线性变化。霍尔传感器(霍尔元件采样)输出相应的脉冲信号。其流量脉冲特性计算公式为:脉冲f(Hz)=8.1x流量Q(L/min)-3。(2)51单片机的中断…

    2022年9月29日
    3
  • qt qtreewidgetitem_qtwidgets

    qt qtreewidgetitem_qtwidgetsQTreeWidget在添加了数据之后点击扩展箭头展开和收缩的时候会显得卡顿,这个是qt内核的一个bug,正常情况下应该是立即响应的,既然是qt的bug,就只能换种方式实现TreeWidget来解决这个bug,具体办法大致可以有一下几种,具体看需求来决定用哪一种方法更好:如果需求要求实现的TreeWidget不要求是完全的TreeWidget,扩展的层次是有限制的,就是只有2层,3层或者…

    2022年10月1日
    3
  • 0603贴片电阻阻值对照表_贴片电阻快速识别表与识别方法[通俗易懂]

    0603贴片电阻阻值对照表_贴片电阻快速识别表与识别方法[通俗易懂]贴片电阻快速识别表与识别方法,要想快速识别贴片电阻,那就必须知道什么是贴片电阻?贴片电阻就是金属粉和玻璃釉粉混合,采用丝网印刷法印在基板上制成的电阻器。与贴片电容不一样的是,贴片电容一般比电容厚,且电阻两边为银白色镀锡,中间黑色。下面我们来看看具体的两个识别方法。快速识别贴片电阻,方法一:如上图所示,正面主要是中间黑色带电阻代码,方便识别贴片电阻大小,两头是锡面,银白色,方便焊接。背面中间是高纯度…

    2022年4月19日
    472

发表回复

您的邮箱地址不会被公开。 必填项已用 * 标注

关注全栈程序员社区公众号