楼主: redaring
5900 3

请教用 proc expand 进行插值的问题 [推广有奖]

  • 0关注
  • 0粉丝

硕士生

21%

还不是VIP/贵宾

-

威望
0
论坛币
48 个
通用积分
0.0600
学术水平
0 点
热心指数
1 点
信用等级
0 点
经验
4910 点
帖子
102
精华
0
在线时间
123 小时
注册时间
2009-11-22
最后登录
2018-2-15

相似文件 换一批

+2 论坛币
k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群
赵安豆老师微信:zhaoandou666

经管之家联合CDA

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

感谢您参与论坛问题回答

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

+2 论坛币
我想对这样的一组数据进行补齐插值,想到用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,可是不太会写,所以达不到自己想要的结果,希望大侠们指点一下。

上面的数据是为了方便说明。传部分数据,感谢各位帮忙!
二维码

扫码加我 拉你入群

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

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

关键词:expand EXP ROC month Data

city.xls

111 KB

沙发
redaring 发表于 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年的日期值,而且并没有差值……

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

使用道具

藤椅
redaring 发表于 2010-1-29 17:41:40 |只看作者 |坛友微信交流群
[cry][cry][cry][cry]

使用道具

板凳
redaring 发表于 2010-1-31 16:07:42 |只看作者 |坛友微信交流群
似乎没有人对这个问题感兴趣啊……

使用道具

您需要登录后才可以回帖 登录 | 我要注册

本版微信群
加好友,备注cda
拉您进交流群

京ICP备16021002-2号 京B2-20170662号 京公网安备 11010802022788号 论坛法律顾问:王进律师 知识产权保护声明   免责及隐私声明

GMT+8, 2024-11-6 09:59