R语言缺失值探索的强大R包：naniar

简介

缺失值在数据中无处不在，需要在分析的初始阶段仔细探索和处理。在本次示例中，会详细介绍naniar包探索缺失值的方法和理念，它和ggplot2和tidy系列使用方法非常相似，上手并不困难。

有时，解释缺失值出现的原因可能很简单，比如，可能是由于记录不全，各种意外等，但实现这一解释的过程可能并不简单，而且可能需要比开发探索性数据分析和模型所需的更多时间。

本次学习主要探讨3个问题：

1. 开始探索缺失值
2. 探索缺失值的机制
3. 模型化缺失值

如何开始探索缺失值

当你面对新的数据时，可能首先会使用各种汇总函数查看数据的基本情况，比如：

• summary()
• str()
• skimr::skim
• dplyr::glimpse()
• ...

但是当数据有缺失值时，就会影响接下来的分析。所以首先还要查看数据的缺失情况。

R包visdat可以展示缺失值数据，主要有2个函数：

• vis_dat()
• vis_miss()

vis_dat()

此函数可以可视化整体数据，还会给出数据类型等信息，当然也包括缺失值。

vis_miss()

vis_miss则主要专注于缺失值的可视化。

探索缺失值的关系

通过vis_miss可以知道哪些变量有缺失值，但是我们还想知道这些缺失值之间有什么关系，这时候就用到naniar包了。

例如：

这幅图会直接把缺失值删掉，并不能知道缺失值的情况。

通过使用naniar包，可以达到以下效果：

缺失值和非缺失值一目了然！

既然是ggplot2一样的tidy系列，那也肯定是支持其他特性的，比如分面：

支持主题：

可视化变量中的缺失值

可以看到每一列有多少缺失值。

当然也是支持更换主题：

支持ggplot2的各种特性，比如labs:

支持分面的语法和上面稍有不同：

还有很多其他可视化方法：

一个更加复杂的数据：

使用NA替换缺失值

在一个数据中有很多缺失值用NA来表示可能会更加方便，比如N/A、N A，Not Available，-999等。

naniar中提供了replace_with_na函数把这些缺失值替换为NA。主要有：

• replace_with_na
• replace_with_na_all
• replace_with_na_at
• replace_with_na_if

和dplyr中的replace_na()用法完全一样，不过一个是把NA替换成其他值，一个是把其他值替换成NA。

整洁的缺失数据：shadow matrix

as_shadow()函数直接以数据框的形式返回是否是缺失值，输入什么样子输出就是什么样子！

把数据框形式的缺失值数据和原数据整合到一起：bind_shadow():

也可以使用nabular()函数达到同样的效果：

这两种方法得到的东西是一样的：

通过这2个函数，就可以对缺失值做一些简单的统计了：

还可以通过画图展示缺失值和非缺失值的数据分布：

可视化插补后的缺失值

使用simpltation包进行缺失值插补，并可视化插补后的数据：

这样就不知道哪些是插补的哪些不是了，可以通过以下方法解决：

缺失值的汇总函数

主要通过n_miss()和n_complete()统计缺失值和非缺失值：

当然还提供其他格式的函数，比如百分比形式，小数形式，表格形式，针对某一列的缺失值汇总等：

汇总表格形式的，以下结果说明：有111行没有缺失值，占数据的72%，有40行只有1个缺失值，占数据的26%，有2行含2个缺失值，只占1%。

还有针对列的：

同样的针对列的缺失值汇总：

除此之外，还提供miss_var_run和miss_var_span函数，还可以和group_by连用探索缺失值！

模型化缺失值

对缺失值建立模型！如果不学习这个R包，我是真的想不到还可以这样搞缺失值！

首先看下缺失值占比：

然后我们可以使用决策树之类的模型来预测哪些变量及其值对于预测缺失比例是重要的：

!!!震惊啊！

这个包还有很多其他的用法，本次教程只是浅尝辄止，不得不说，包的作者真是太有才了，缺失值也能玩出花来！

本文转自医学与生信笔记。

如果需要进群，请在公众号菜单栏处添加作者微信，邀请入群，请记得附上备注，谢谢。

点个在看，SCI马上发表。