r 函数返回多个值_第四讲 R描述性统计分析

00f9b8ca633044187100a95163c8ac5a.gife28f814c8a22e73413887fbda63af508.gif

在“R与生物统计专题”中,我们会从介绍R的基本知识展开到生物统计原理及其在R中的实现。以从浅入深,层层递进的形式在投必得医学公众号更新。

在上一讲中,我们介绍了第三讲 R编程基础-矩阵和数据框 (戳这里即可跳转)到现在为止,大家已经学完了R语言的基础知识,马上就开始进阶地学习一些R语言相关的统计学技能了。

今天的更新,我们会带您学习R的基本统计学技能:描述性统计分析。

1. 将数据导入R

1.1 准备好你的数据命名约定
  • 避免名称带有空格。

        好的列名:patient_age或patient.age。

        列名错误:patient age。

  • 避免使用带有特殊符号的名称:?,$,*,+,#,(,),-,/,},{,|,>,

  • 避免以数字开头的变量名。请改用字母。

    好的列名称:patient_1st_meal或x1st_meal。

    列名错误:1st_male

  • 列名必须是唯一的。不允许重复的名称。

  • R区分大小写。这意味着名称不同于名称或名称。

  • 避免数据中出现空白行。

  • 删除文件中的任何评论。

  • 用NA替换缺少的值(不可用)。

  • 如果你有包含日期的列,请使用四位数格式

    格式良好:20160101。

    格式错误:01/01/16

1.2 将数据保存在外部.txt标签或.csv文件中1.3 如下将数据导入R
# 对于.txt文件my_data 
# 对于.csv文件my_data 

在这里,我们将使用名为iris的内置R数据集。

# 导入R内自带的iris数据集library(datasets)data(iris)
# 将数据存储在变量my_data中my_data 
1.3 检查数据

你可以使用head()和tails()函数检查数据,这将分别显示数据的第一部分和最后一部分。

# 显示前六行内容head(my_data, 6)

输出结果如下

Sepal.Length Sepal.Width Petal.Length Petal.Width Species

1          5.1         3.5          1.4       0.2  setosa

2          4.9         3.0          1.4       0.2  setosa

3          4.7         3.2          1.3       0.2  setosa

4          4.6         3.1          1.5       0.2  setosa

5          5.0         3.6          1.4       0.2  setosa

6          5.4         3.9          1.7       0.4  setosa

2. 常用的描述性统计的R函数

一些用于计算描述性统计量的R函数:

✎ R函数 

6f13dcce22dc0b020f45e5bc3553d839.png

3. 单个组的描述性统计

3.1 集中趋势的度量:均值,中位数,众数

粗略地说,集中趋势衡量的是数据的“平均”或“中间”。最常用的衡量指标包括:

集中趋势

平均值:平均值。它对异常值很敏感。

中位数:中间值。这是一个强有力的替代手段。

众数:最频繁出现的值

在R中

  • 函数mean()和median()可以分别计算平均值和中位数;

# 计算平均值mean(my_data$Sepal.Length)

[1] 5.843333

# 计算中位数median(my_data$Sepal.Length)

[1] 5.8

3.2 可变性的度量

可变性度量给出了数据“分散”的程度。

范围

极值:最小值和最大值

范围:最大值减去最小值

# 计算最小值min(my_data$Sepal.Length)

[1] 4.3

# 计算最大值max(my_data$Sepal.Length)

[1] 7.9

# 范围range(my_data$Sepal.Length)

[1] 4.3 7.9

四分位间距

四分位数将数据均匀分为4部分。四分位数间距(IQR):对应于第一和第三四分位数之间的差异-有时被用作标准偏差的可靠替代方案。

  • R功能:

quantile(x, probs = seq(0, 1, 0.25))
  • x:需要样本分位数的数值向量。

  • probs:在[0,1]之间的概率数值向量。

  • 例:

quantile(my_data$Sepal.Length)

0%  25%  50%  75% 100%

4.3  5.1  5.8  6.4  7.9

# 计算十分位数(0.1,0.2,0.3,…,0.9):quantile(my_data$Sepal.Length, seq(0, 1, 0.1))
# 计算四分位间距:IQR(my_data$Sepal.Length)

[1] 1.3

方差和标准差

方差表示与均值的平均平方差之和。标准差是方差的平方根。它测量数据中数值与平均值的平均偏差。

# 计算方差var(my_data$Sepal.Length)
# 计算标准差sd(my_data$Sepal.Length)
绝对中位数

绝对中位数(Median absolute deviation,MAD):数据中值与中值的偏差,即先计算出数据与它们的中位数之间的残差(偏差),MAD就是这些偏差的绝对值的中位数。

# 计算中位数median(my_data$Sepal.Length)
# 计算绝对中位数mad(my_data$Sepal.Length)
统计描述方式的选择
  • 范围。它不经常使用,因为它对异常值非常敏感。

  • 四分位间距。对于异常值,它非常强大。它多与中位数结合使用。

  • 方差。完全无法解释的,因为它不使用与数据相同的单位。除了用作数学工具外,很少被使用。

  • 标准偏差。方差的平方根。它以与数据相同的单位表示。在均值是集中趋势的分布(多指正态分布)的情况下,通常使用标准偏差。

  • 绝对中位数。对于具有离群值的数据,这是一种估算标准偏差的可靠方法。但是不经常使用。

  • 总而言之,四分位间距和标准差是用于报告数据变异性的两种最常用的度量。

3.3 计算变量和整个数据框的整体摘要summary()函数

函数summary()可用于显示一个变量或整个数据框的多个统计变量概况。

  • 单个变量的概况。

返回六个值:平均值,中位数,第25和,75四分位数,最小值和最大值。

summary(my_data$Sepal.Length)

输出结果如下:

Min. 1st Qu.  Median    Mean 3rd Qu.    Max.4.300   5.100   5.800   5.843   6.400   7.900
数据框概况

在这种情况下,函数summary()将自动应用于每列。结果的格式取决于列中包含的数据类型。例如:

  • 如果列是数字变量,则返回均值,中位数,最小值,最大值和四分位数。

  • 如果该列是一个因素变量(factor),则返回每个组中的观察数。

summary(my_data, digits = 1)

输出结果如下:

Sepal.Length  Sepal.Width  Petal.Length  Petal.Width        SpeciesMin.   :4     Min.   :2    Min.   :1     Min.   :0.1   setosa    :501st Qu.:5     1st Qu.:3    1st Qu.:2     1st Qu.:0.3   versicolor:50Median :6     Median :3    Median :4     Median :1.3   virginica :50Mean   :6     Mean   :3    Mean   :4     Mean   :1.23rd Qu.:6     3rd Qu.:3    3rd Qu.:5     3rd Qu.:1.8Max.   :8     Max.   :4    Max.   :7     Max.   :2.5
3.4 缺失值的情况Tips

当数据包含缺失值时,即使仅缺少一个值,某些R函数也会返回错误或NA。

例如,即使向量中仅丢失一个值,mean()函数也将返回NA。使用参数na.rm = TRUE可以避免这种情况,该参数告诉函数在计算之前删除所有NA。使用均值函数的示例如下:

mean(my_data$Sepal.Length, na.rm = TRUE)

好了,本期讲解就先到这里。

在之后的更新中,我们会进一步为您介绍R的入门,以及常用生物统计方法和R实现。欢迎关注,投必得医学手把手带您走入R和生物统计的世界。

提前打个预告,接下来我们要正式开始学习R语言的统计学技能啦,下一期将会更新“R的描述性统计分析”。喜欢的同学们快快关注起来吧。

65a047d11be7a21743682c7e4859d09d.gif第一讲 R-基本介绍及安装第二讲 R-编程基础-运算、数据类型和向量等基本介绍第三讲 R编程基础-矩阵和数据框

当然啦,R语言的掌握是在长期训练中慢慢积累的。一个人学习太累,不妨加入“R与统计交流群”,和数百位硕博一起学习。

快扫二维码撩客服,

带你进入投必得医学交流群,

让我们共同进步!

↓↓

19af7fccb9434fc33927247039f19acc.png

- END -

长按二维码关注「投必得医学」,更多科研干货在等你!

ebb4a549ecb0fd0d5894297a696edd60.png

麻烦点一下在看再走呗

5ffeb3be29ff1afa7a654d87d4a05e88.gif

本文来自互联网用户投稿,该文观点仅代表作者本人,不代表本站立场。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如若转载,请注明出处:http://www.mzph.cn/news/538964.shtml

如若内容造成侵权/违法违规/事实不符,请联系多彩编程网进行投诉反馈email:809451989@qq.com,一经查实,立即删除!

相关文章

python双线性插值函数_双线性插值法原理 python实现

码字不易,如果此文对你有所帮助,请帮忙点赞,感谢!一. 双线性插值法原理:① 何为线性插值?插值就是在两个数之间插入一个数,线性插值原理图如下:在位置 x 进行线性插值,插…

张恩民 php,php100教程_php100张恩民php视频 56课时_网站编程_寂涯网络学习基地

1、PHP基本语法和数据类型:(1)、PHP基本语法:1、htm 和 php 混编2、一个语句以 ; (分号结束)3、如何定义一个变量,和变量的使用(2)、PHP数据运算类型四种标量类型:boolean(布尔型)integer(整型)float(浮点型,也作“do…

Oracle 自带sqldeveloper导入数据到表

步骤 1、建一张字段名和类型与要导入数据一致的表。 2、如果是excel文件,另存为 csv文件。 因为excel文件遇到空格会报错,excel文件太大工具会读取文件失败。 而csv文件是秒级读取,有空格也会导入不会报错。 ------------- 更多的Java&am…

百度腾讯中兴华为全部入局Linux Foundation深度学习基金会

2019独角兽企业重金招聘Python工程师标准>>> 虽然名字叫做Linux Foundation,但是这家基金会关注的领域并不局限在Linux操作系统领域。事实上,Linux Foundation基金会还为很多其他开源基金和项目提供支持,比如Cloud Foundry、Autom…

python判断点在矩形内_判断平面内矩形和圆是否有交点的python实现

#! /usr/bin/env python# -*- coding: utf-8 -*-filename : 平面坐标内矩形与圆是否相交usage :authors : U{peterguo}copyright: tencentdate : 2012.08.21version : 1.0.0.0Modified History:2012.08.21 创建import osimport sysimport reimport urllib2import timeimport da…

python读取mysql数据库_Python操作MySQL数据库9个实用实例

在Windows平台上安装mysql模块用于Python开发 用python连接mysql的时候,需要用的安装版本,源码版本容易有错误提示。下边是打包了32与64版本。 MySQL-python-1.2.3.win32-py2.7.exe MySQL-python-1.2.3.win-amd64-py2.7.exe 实例 1、取得 MYSQL 的版本 #…

win主机上搭建php网站运行环境,Windows server 2008搭建php运行环境图文详解(php5.3)

这篇文章主要为大家分享下Windows server 2008搭建php运行环境的步骤,需要的朋友可以参考下下载php组件包首先到http://windows.php.net/download/下载你需要的php版本,这里我下载的是php5.3。下面解压php组件 包到磁盘上。安装Microsoft Visual C 2008 …

python 登陆微博 被删除 token_爬取微博信息,使用了cookie仍然无法登录微博

按照网上的模板自己写了类似的代码爬取微博,可是response回来的html是登录界面的html,应该是没有成功登陆微博,但是和网上的代码是基本一样的from bs4 import BeautifulSoupimport pandas as pdimport requestsurl "https://weibo.cn/2…

plsql developer连接oracle--本地不安装oracle

工具 数据库开发程序 PLSQL Developer oracle 客户端 instantclient,点击下载 (下载与oracle相同版本或最新版本,如果 PLSQL Developer是32位,则要下载32位的)安装 1、解压oracle客户端 到盘某个目录。我的目录是C:\or…

ubuntu networking 与 network-manager

刚遇到的坑,因为操作不当导致网络中断,于是手动配置了/etc/network/interfaces , 修复了系统之后发现ubuntu-desktop中的有线链接不见了,百度了一下说是networking服务和network-manger服务会起冲突。 【Server版本】 在UbuntuServer版本中&…

mysql 升序_mysql 的 查找 与 排序

mysql 的 查找 &#xff1a;SELECT name,age FROM mytable; # 在mytable 中所有的name,age 数据SELECT name,age FROM mytable WHERE age>25; #在mytable 中所有的name,age 数据&#xff0c;其中age必须小于25SELECT name,age FROM mytable WHERE age<25 OR age>30; #…

Hadoop3.0的新特性

转载&#xff1a;http://blog.csdn.net/sinat_31726559/article/details/52103229 点击阅读原文 -------------------------------------- 1. Hadoop3.0简介 Hadoop 2.0是基于JDK 1.7开发的&#xff0c;而JDK 1.7在2015年4月已停止更新&#xff0c;这直接迫使Hadoop社区基于…

python torch exp_学习Pytorch过程遇到的坑(持续更新中)

1. 关于单机多卡的处理&#xff1a;在pytorch官网上有一个简单的示例&#xff1a;函数使用为&#xff1a;torch.nn.DataParallel(model, deviceids, outputdevice, dim)关键的在于model、device_ids这两个参数。DATA PARALLELISM​pytorch.org但是官网的例子中没有讲到一个核心…

php 读文件返回字符串,PHP:file_get_contents('php:// input')返回JSON消息的字符串...

我正在尝试在我的PHP应用程序中读取JSON消息&#xff0c;这是我的php代码&#xff1a;$json file_get_contents(php://input);$obj json_decode($json, TRUE);echo $obj->{S3URL};当我这样做时&#xff0c;出现以下错误&#xff1a;Trying to get property of non-object …

Android中ListView数据使用sAdapter.notifyDataSetChanged();方法不刷新的问题

原文链接&#xff1a;http://blog.csdn.net/caihongdao123/article/details/51513410 点击阅读原文 ------------------------- 1.涉及到数据库 当要动态显示更数据库改动&#xff0c;相信大家应该都用过notifyDataSetChanged();. 例如&#xff1a; ...... …

keepalived配置高可用集群

准备工作 分别在主从上安装keepalived和nginxyum install -y keepalivedyum install -y nginx关闭主从上的防火墙和SELinuxsystemctl stop firewalldsetenforce 0 配置主机 查看主机ip [rootlynn-04 ~]# ifconfig ens33: flags4163<UP,BROADCAST,RUNNING,MULTICAST> mtu…

如何快速掌握python包_如何快速掌握一个python模块?

初学者就别想快了。 我自己是这样的。先上网看看一些基础的教程&#xff0c;非常快的过一下&#xff0c;十几分钟&#xff0c;主要是了解这个module能干什么&#xff0c;特别是一些基本的功能&#xff0c;头脑中建立起初步映射。 然后就是用了&#xff0c;不用看了也白搭。我假…

python设计一个函数定义计算并返回n价调和函数_python函数的调和平均值?

我有两个函数&#xff0c;给出精度和召回分数&#xff0c;我需要做一个调和平均函数&#xff0c;定义在同一个库中&#xff0c;使用这两个分数。函数如下所示&#xff1a;功能如下&#xff1a;def precision(ref, hyp):"""Calculates precision.Args:- ref: a l…

jsp超链接到java文件,jsp页面超链接传中文终极解决办法

在做web前端页面的时候&#xff0c;经常碰到传中文问题。网上有许多方案&#xff0c;但仍不能根治&#xff0c;最终要用js或者java的encode相关方法。常规方案有三部&#xff1a;1.改tomcat的server.xml中URIEncodeing为utf-82.页面中编码设置为utf-83.整个项目编码使用utf-8我…

自定义ListView中的分割线

原文&#xff1a;http://blog.csdn.net/zuolongsnail/article/details/7187302点击阅读 --------------------------------------- ListView中每个Item项之间都有分割线&#xff0c;设置Android:footerDividersEnabled表示是否显示分割线&#xff0c;此属性默认为true。 1.不显…