ggmice | 用这只可爱的小老鼠来填补你的缺失值吧!~
创始人
2024-04-06 14:59:47
0

1写在前面

我们在处理数据的时候常常会遇到存在缺失值(NA)的情况,如何处理就仁者见仁智者见智了。🤒
简单粗暴的方法可能就是行删除法(listwise)或者个案删除法(case-wise)了,这种方法在缺失值比较少的情况下比较适用,但在NA较多的情况下可能就会丢失过多信息导致无法继续分析。😘
本期我们介绍一下mice包ggmice包这两只可爱的小老鼠,全名Multivariate Imputation by Chained Equations, mice,即链式方程多重填补。📍

一张图总结基本原理,嘿嘿。👇

alt

2用到的包

rm(list = ls())
library(tidyverse)
library(mice)
library(ggmice)

3示例数据

这里我们使用一下示例数据airquality,再在其中添加一些缺失值

dat <- airquality
dat[4:9,3] <- rep(NA,6)
dat[1:4,4] <- NA
alt

4数据概览

4.1 缺失值查看

这里面我们一共有6个变量,其中4个存在缺失值。🫠

summary(dat)
alt

4.2 缺失值可视化

Note! ggmice提供了一种NA值的可视化方法,一目了然,nice! 🤒

plot_pattern(dat,
square = F,
rotate = F)
alt

4.3 influx-outflux plot

这里和大家简单介绍一下这个influx-outflux plot,总的来说评估了缺失数据其他变量联系程度。😂
一般来说,在建模时,influx以及outflux越大越好。😗

plot_flux(dat,
label = F,
caption = F)
alt

5可视化一下吧

5.1 连续变量

这里我们对连续变量缺失值进行一下可视化,可以看到红色的为缺失值。😘

ggmice(dat, aes(Ozone, Solar.R))+
geom_point()
alt

5.2 分类变量

接着我们对分类变量缺失值进行一下可视化,红色的为缺失值。😃

ggmice(dat, aes(Month, Solar.R)) +
geom_point()
alt

5.3 分面展示

ggmice(dat, aes(Month, Solar.R)) +
geom_point() +
facet_wrap(~ Month == 5,
# labeller = label_both
)
alt

6mice包填补缺失值

6.1 填补缺失值

在这里我们生成几个填补缺失值后的数据,m默认是5,为了减小计算量,这里我设置成3。🤗
Note! 可选method包括:👇

pmm,
logreg,
polyreg,
polr

imp <- mice(dat, m = 3, method = "pmm")
alt

6.2 连续变量缺失值填补后可视化

我们再看一下填补缺失值后的散点图吧,红色的为缺失值填补后。😚

ggmice(imp, aes(Ozone, Solar.R))+
geom_point()
alt

6.3 分类变量缺失值填补后可视化

ggmice(imp, aes(Month, Solar.R)) +
geom_point()
alt

6.4 分面展示

ggmice(dat, aes(Month, Solar.R)) +
geom_point() +
facet_wrap(~ Month == 5,
# labeller = label_both
)
alt

7填补数据集的可视化

7.1 dotplot

我们之前设置了m = 3,这里我们看一下3个数据集的NA填补情况。🥰

ggmice(imp, aes(x = .imp, y = Temp)) +
geom_jitter(height = 0, width = 0.25) +
labs(x = "Imputation number")
alt

7.2 boxplot

ggmice(imp, aes(x = .imp, y = Temp)) + 
geom_jitter(height = 0, width = 0.25) +
geom_boxplot(width = 0.5, size = 1, alpha = 0.75, outlier.shape = NA) +
labs(x = "Imputation number")
alt

8算法收敛

看来default = 5是有原因的,哈哈哈哈哈哈!🤣

plot_trace(imp,
# "Temp"
)
alt

甜筒
最后祝大家早日不卷!~

点个在看吧各位~ ✐.ɴɪᴄᴇ ᴅᴀʏ 〰

📍 往期精彩

📍 🤩 ComplexHeatmap | 颜狗写的高颜值热图代码!
📍 🤥 ComplexHeatmap | 你的热图注释还挤在一起看不清吗!?
📍 🤨 Google | 谷歌翻译崩了我们怎么办!?(附完美解决方案)
📍 🤩 scRNA-seq | 吐血整理的单细胞入门教程
📍 🤣 NetworkD3 | 让我们一起画个动态的桑基图吧~
📍 🤩 RColorBrewer | 再多的配色也能轻松搞定!~
📍 🧐 rms | 批量完成你的线性回归
📍 🤩 CMplot | 完美复刻Nature上的曼哈顿图
📍 🤠 Network | 高颜值动态网络可视化工具
📍 🤗 boxjitter | 完美复刻Nature上的高颜值统计图
📍 🤫 linkET | 完美解决ggcor安装失败方案(附教程)
📍 ......

本文由 mdnice 多平台发布

相关内容

热门资讯

AWSECS:访问外部网络时出... 如果您在AWS ECS中部署了应用程序,并且该应用程序需要访问外部网络,但是无法正常访问,可能是因为...
AWSElasticBeans... 在Dockerfile中手动配置nginx反向代理。例如,在Dockerfile中添加以下代码:FR...
AWR报告解读 WORKLOAD REPOSITORY PDB report (PDB snapshots) AW...
AWS管理控制台菜单和权限 要在AWS管理控制台中创建菜单和权限,您可以使用AWS Identity and Access Ma...
北信源内网安全管理卸载 北信源内网安全管理是一款网络安全管理软件,主要用于保护内网安全。在日常使用过程中,卸载该软件是一种常...
​ToDesk 远程工具安装及... 目录 前言 ToDesk 优势 ToDesk 下载安装 ToDesk 功能展示 文件传输 设备链接 ...
Azure构建流程(Power... 这可能是由于配置错误导致的问题。请检查构建流程任务中的“发布构建制品”步骤,确保正确配置了“Arti...
群晖外网访问终极解决方法:IP... 写在前面的话 受够了群晖的quickconnet的小水管了,急需一个新的解决方法&#x...
AWSECS:哪种网络模式具有... 使用AWS ECS中的awsvpc网络模式来获得最佳性能。awsvpc网络模式允许ECS任务直接在V...
不能访问光猫的的管理页面 光猫是现代家庭宽带网络的重要组成部分,它可以提供高速稳定的网络连接。但是,有时候我们会遇到不能访问光...