全部版块 我的主页
论坛 数据科学与人工智能 数据分析与数据科学
6128 3
2010-01-29
我想对这样的一组数据进行补齐插值,想到用proc expand,可是似乎实现不了,希望大侠们帮忙一下,谢谢!

CITY    YEAR    MONTH    DAY    DATA
aaaa    1991      04             01        12345
aaaa    1991      04             02        53472
aaaa    1991      04             03        23456
aaaa    1991      04             05        25436
aaaa    1991      04             06        56743
aaaa    1991      04             08        26754
aaaa    1991      04             09        74532
aaaa    1991      04             10        86435
......
......
bbbb    1991      04             01        12345
bbbb    1991      04             02        53472
bbbb    1991      04             03        23456
bbbb    1991      04             05        25436
bbbb    1991      04             07        56743
bbbb    1991      04             08        26754
bbbb    1991      04             12        74532
bbbb    1991      04             14        86435
......
数据中有几万个城市,每个城市按时间序列有数据,是每天的数据,有很多年,并且每个城市数据开始和结束的年月日都有可能不一样,这些时间是不连续的,其中有很多缺失值,可是数据量太大,不可能一个个去找。我想将这些缺失的数据补齐,并且把data的数据插值补上,另外,我想限定一个条件,只有连续缺失7天的数据用这个方法补,否则就留空(以空值代替)。
我想到要用proc expand,可是不太会写,所以达不到自己想要的结果,希望大侠们指点一下。

上面的数据是为了方便说明。传部分数据,感谢各位帮忙!
附件列表

city.xls

大小:111 KB

 马上下载

二维码

扫码加我 拉你入群

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

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

全部回复
2010-1-29 09:49:31
抛一块砖,希望各位大侠帮忙

这是我自己写的,完全达不到我想要的结果,并且还很有问题,呵呵

data city;
set test;           ;
date = MDY (month,day,year);
format date date9.;
run;

proc sort data=city
by city date;
run;

data range;
  retain from to;
  set city end=lastobs;
  if _n_=1 then from=date;
  if lastobs then do;
      to = date;
      output;
      end;
run;

data dates;
  set range;
  date=from;
  do while (date<= to);
      output;
      date=intnx('day', date,1);
      end;
  keep date;
run;

proc standard;
  merge city dates;
  by city;
run;

proc expand data=city  out=temp  from=day method=join;
  var data;
  by city;
  id date;
run;

我想法是先把日期转换成sas的日期格式,然后把不连续的时间补成连续的,然后用expand差值。
可是出现的问题,在生成新的时间序列表的时候,只生成一个日期序列,并没有按照每个城市生成时间序列。expand后还会生成一列1960年的日期值,而且并没有差值……

我不知道怎么改,希望各位指点。谢谢!
二维码

扫码加我 拉你入群

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

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

2010-1-29 17:41:40
[cry][cry][cry][cry]
二维码

扫码加我 拉你入群

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

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

2010-1-31 16:07:42
似乎没有人对这个问题感兴趣啊……
二维码

扫码加我 拉你入群

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

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

相关推荐
栏目导航
热门文章
推荐文章

说点什么

分享

扫码加好友,拉您进群
各岗位、行业、专业交流群