随笔5：R语言插补面板数据缺失值

原创

已于 2022-05-05 09:57:27 修改 · 3.4k 阅读

CC 4.0 BY-SA版权

文章标签：

于 2021-12-07 22:40:03 首次发布

【2022.5.5说明】增加了一行命令清除内存中的临时数据集

【linput3.0.6】百度网盘链接：https://pan.baidu.com/s/1SeLBKp933RKFu_J01O6eIA
提取码：xna4

Github地址见后文

【2022.5.4说明】

1.增加了测试数据集，用library命令加载包后，输入data("test",package = "linput") 即可加载包内置数据test到内存

2.现在补齐的向量值会被添加到原始数据集中作为新的一列，提高了科研效率。

【linput3.0.5 R包】百度网盘链接：链接：https://pan.baidu.com/s/1Ik8jn1k8R2y-nKlcg9EDKw
提取码：yln2

Github下载地址见后文。

=====================================================================

【2022.5.2说明】近来较多朋友反映包的使用出了一些问题，经过和北京林业大学经济管理学院的梁晨同学多次讨论和测试，推出本R包的3.0.4版。

【linput3.0.4包的说明】修复了评论区中出现的BUG。

【linput3.0.4包作者信息】1.杜云晗.西南财经大学经济学院政治经济学 2.梁晨.北京林业大学经济管理学院经济统计学

【其他说明】linput包的m参数设置为0即为只用内插法，为1即为用了内插法后用外推法。

外推法基于第一个有数值的年份和最后一个有数值的年份进行年均增长率测算然后外推，如果第一个有数值的年份大于等于最后一个有数值的年份，如500(1987)和300（1997）；500(1987)和500(1997)等情形则外推失败。

【使用过程中有什么问题，请及时将测试数据发到1270950448@qq.com，方便我测试BUG】这是本人自编的一个简单的用于处理面板数据缺失值的R语言包，原理是对存在缺失值的面板数据的按顺序年份出现的第一个非缺失值，以及倒数的第一

19 条评论

2301_77007719 2023.04.06
能请教一下，为什么使用之后个别城市的数据会从2013-2021，变成2019 2020 2021 2013 2014...2018，然后中间会空出几行
- 2301_77007719回复2301_77007719 2023.04.07
  好了，找到问题了：数据有问题，您的包非常好用

qq_43891750 2023.01.02
大佬，用外插法的时候一直报这个错，但内插的时候是正常的： Error in if (growth != 0 && test_data2$year[i + step_long] < fir_temp_year & : missing value where TRUE/FALSE needed
- 女神の戦士回复qq_43891750 2023.01.18
  请发数据和测试代码到我邮箱1270950448@qq.com

weixin_43753134 2022.10.17
a<-c("四川省","四川省","四川省","四川省","四川省","四川省","湖北省","湖北省", "湖北省","湖北省","湖北省","湖北省") b<-c(rep(1999:2004,2)) c<-c(NA, 18,19,NA,20,NA,20,21,NA,22,NA,NA) df<-as.data.frame(cbind(a,b,c)) df$b<-as.numeric(df$b) df$c<-as.numeric(df$c) linput(df,df$a,df$b,df$c,m=1) #利用内插法插补缺失值还是出现报错： Note:This package #Linear_inputed_Duyunhan# need the package #tidyverse#'s supporting,now the package 'tidyverse' has already been installed. Error in Summary.factor(c(2L, 2L, 2L, 2L, 2L, 2L, 1L, 1L, 1L, 1L, 1L, : ‘max’ not meaningful for factors 请问是为什么呢？

alwaysabstract 2022.06.15
Note:This package #Linear_inputed_Duyunhan# need the package #tidyverse#'s supporting,now the package 'tidyverse' has already been installed. Error in 1:max_year : result would be too long a vector In addition: Warning messages: 1: In linput(newpaneldata, newpaneldata$label, newpaneldata$DATE, newpaneldata$GDP, : NAs introduced by coercion 2: In linput(newpaneldata, newpaneldata$label, newpaneldata$DATE, newpaneldata$GDP, : NAs introduced by coercion 3: In max(test_data2$year, na.rm = T) : no non-missing arguments to max; returning -Inf 4: In min(test_data2$year, na.rm = T) : no non-missing arguments to min; returning Inf 您好，请问运行代码时一直出现这个错误，是为啥呢，求大佬指教
- 女神の戦士回复alwaysabstract 2022.09.08
  有可能，根据你的报错结果，某些临时变量出现负值了
- alwaysabstract回复女神の戦士 2022.06.26
  我怀疑应该是我的数据空值太多了
- 女神の戦士回复alwaysabstract 2022.06.15
  可以发样本数据到我邮箱1270950448@qq.com,我看看哪里出了BUG
- 女神の戦士回复alwaysabstract 2022.06.15
  把数据和你的代码发到1270950448@qq.com我测试下

qq_45374983 2022.02.28
我的也出现了一楼的那种情况，而且除了省份其他都是数值型[[face]emoji:010.png[/face]
- 女神の戦士回复qq_45374983 2022.04.11
  可以发样本数据到我邮箱1270950448@qq.com,我看看哪里出了BUG，我自己的数据一直在使用这个包但没出过问题

法圆 2021.12.24
> linput(mydataframe,mydataframe$province,mydataframe$year,mydataframe$nox,m=1) Note:This package #Linear_inputed_Duyunhan# need the package #tidyverse#'s supporting,now the package 'tidyverse' has already been installed. Error in Summary.factor(c(2L, 2L, 2L, 2L, 2L, 2L, 27L, 27L, 27L, 27L, : ‘max’ not meaningful for factors 您好，我对数据的处理情况与您的相似，但是我在最后两行代码的部分一直出现这样的错误，麻烦可以帮我解答一下吗？
- LQ666A回复法圆 2022.04.15
  解决了吗？我也是这个问题，
- 女神の戦士回复法圆 2021.12.27
  应该是你的数据结构中某个变量数值类型不对，我这个函数封装的时候，里面涉及到最大最小值的判断，所以你的目标变量必须是数值型变量，比如，year，nox，都必须是数值型，用as.numeric转化下吧。实在不行发数据到1270950448@qq.com我看