关于缺失值插补法的一个例子

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

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

立即领取

感谢您参与论坛问题回答

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

+2 论坛币

今天在sql数据库里面查询的时候，发现某列数据有比较多的缺失值，要求该列数据的缺失值用前面的非缺失值填充。但前一个可能也是缺失的，用sql写了好久，算是搞定了，可惜公司数据库没权限连接R，不然就简单了。如有一列数据 y<-c(10,NA,NA,50,NA)
则实现的编程为
for(i in 1:5) {
if (is.na(y[i])==TRUE) y[i]<-y[i-1]
}



用sas的话，则写如下代码
data e;                                                                                                                                
set a;                                                                                                                                 
retain z;                                                                                                                              
if y^=. then z=y;      

都比坑爹的sql好用多了……
大家如果有向量化的方法实现这个效果，欢迎回帖。

扫码加我 拉你入群

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

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

分享0 收藏0 回帖

关键词：缺失值 SQL数据库 数据库 sql sql数据库

看这个帖子http://cos.name/cn/topic/129915/

这个例子有bug啊   如果y=c(NA,10,NA,50,NA)    就不对了

yywan0913 发表于 2014-11-14 10:55
这个例子有bug啊   如果y=c(NA,10,NA,50,NA)    就不对了

恩，我的数据第一行是有值的，就没考虑第一行没值的情况

李会超 发表于 2014-11-13 17:57
今天在sql数据库里面查询的时候，发现某列数据有比较多的缺失值，要求该列数据的缺失值用前面的非缺失值填充 ...

length.x = 20000
x = runif(length.x, min = 1, max = 100)
x[ sample(2:length.x, 6000) ] = NA

setNA1 = function(x) {
    x[ is.na(x) ] = x[ which(is.na(x)) - 1 ]
    if (any(is.na(x))) return(setNA1(x)) else return(x)
}

setNA2 = function(x) {
    for (i in 2:length(x)) if (is.na(x[ i ])) x[ i ] = x[ i - 1 ]
    return(x)
}

system.time(setNA1(x))
system.time(setNA2(x))

从耗费时间的对比上，方法1要略优于方法2。