求高手指点R语言对缺失数据的多重插补程序

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

立即领取

感谢您参与论坛问题回答

经管之家送您两个论坛币！

+2 论坛币

      RT，本人崭新手一枚，往各位大师不吝赐教。

      例：1995-2005年间，有对应0~20岁的死亡率（即矩阵），但现在2000年的死亡率数据完全缺失，用R如何多重插补？

      另存储的数据（txt文件）如果有缺失，R好像读不出来data？那该如何进行下去呢？

      非常非常感谢！！

扫码加我 拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

分享0 收藏0 回帖

关键词：缺失数据 高手指点 多重插补 R语言 求高手 程序

建议是不要补了，从样本区间和缺失数据年份看，正好是中间,可以做前后对比分析嘛。插值之后的2000年数据和1999年的2001年的数据都相关，浪费两年的数据，不值。

欢迎各位推荐一些书~~找了好久没找着合适的，十分感谢~

nuomin 发表于 2012-8-30 19:30
建议是不要补了，从样本区间和缺失数据年份看，正好是中间,可以做前后对比分析嘛。插值之后的2000年数据和1 ...

谢谢您的提醒~~只是接下来需用到完整的时间序列，选择之下也是多重插补（涉及回归分析）最合适了。。

有木有人救我啊~~

你的问题还没有解决吗

推荐《R语言实战》-《R in action》

用mice包
library(mice)
data_imp<-mice(data,m)
fit<-with(train_imp,analysis)

#m是插补数据集
#analysis 是一个表达式对象，设定统计方法

1. install.packages("mice")
   
2. imp <- mice(data)
   
3. fit <- with(imp,formula)
   
4. pooled <- pool(fit)

复制代码

以上是插补，得出的最终结果在pooled中
通过

1. pooled$qbar

复制代码

可以查看回归系数

ヾ過去式°ο° 发表于 2015-12-20 14:04
以上是插补，得出的最终结果在pooled中
通过
可以查看回归系数

你好，有两个问题：with里面的formula是用来干嘛的啊；另外生成了5个插补数据集，最后怎么生成一个最终的插补好的数据集？谢谢赐教啊，实在不懂了

ヾ過去式°ο° 发表于 2015-12-20 14:04
以上是插补，得出的最终结果在pooled中
通过
可以查看回归系数

你好，有两个问题：with里面的formula是用来干嘛的啊；另外生成了5个插补数据集，最后怎么生成一个最终的插补好的数据集？谢谢赐教啊，实在不懂了