【请教】怎么用SAS计算每N年的一个指标（比如连续五年的平均值、连续三年的标准差等） - 经管之家

› 论坛 › 数据科学与人工智能 › 数据分析与数据科学 › SAS专版

【请教】怎么用SAS计算每N年的一个指标（比如连续五年的平均值、连续三年的标准差等）

13022

29

收藏 2014-11-09

请教大侠们，怎么用SAS计算每N年的一个指标（比如连续五年的平均值、连续三年的标准差等）
具体而言：
如果要计算这个指标，又怎么做？
将公司利润EARNING每三年加总一次。（数据详见附件）

附件列表

大小:185.5 KB

本附件包括：

test1.csv

二维码

扫码加我拉你入群

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

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

全部回复

2014-11-9 21:31:00

用sum+earning得到data a;
merge a a(firstobs=6  rename=(sum=sum_1));
ma5=(sum_1-sum)/5;
即用6年总利润-第一年利润然后运算  标准差类似

二维码

扫码加我拉你入群

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

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

2014-11-9 21:32:21

我理解的不一定对

二维码

扫码加我拉你入群

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

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

2014-11-9 21:37:25

我理解的不一定对

二维码

扫码加我拉你入群

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

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

2014-11-9 21:39:30

朱蕾Julie 发表于 2014-11-9 21:32
我理解的不一定对

感谢你的热心回答！我试了试，还是有点没有摸清您的思路，现在我把一个很简洁的数据集（test）放在了帖子的附件上，能否请您给一个代码呢？

二维码

扫码加我拉你入群

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

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

斤斤计较你

2014-11-9 22:38:39

参考一下：

data class;
  set sashelp.class;
  array var_lag[6];
  var_lag[1] = lag(age);
  do i=1 to 5;
  var_lag[i+1] = lag(var_lag[i]);
  end;
  mean = mean(of var_lag[*]);
  std = std(of var_lag[*]);
run;

二维码

扫码加我拉你入群

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

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

点击查看更多内容…

2014-11-10 07:46:07

斤斤计较你发表于 2014-11-9 22:38
参考一下：

data class;

非常好的举例，给我很大帮助！

二维码

扫码加我拉你入群

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

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

2014-11-10 08:08:23

proc import out=test
datafile="D:\test1.csv"
dbms=csv replace;
getnames=yes;
run;
proc sort data=test out=test1;
by year;
run;
data a;
set test1;
by year;
if last.year ;
sum+earning;
sum_1=sum-earning;
run;
data b;
merge a a(firstobs=3 drop=sum_1 earning rename=(sum=sum_2)) ;
ma=(sum_2-sum_1);
run;

二维码

扫码加我拉你入群

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

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

2014-11-10 09:03:15

朱蕾Julie 发表于 2014-11-10 08:08
proc import out=test
datafile="D:\test1.csv"
dbms=csv replace;

朱蕾的程序很好，但是如果能考虑stkcd就更好了。因为要针对每一个公司来做，比如stkcd=000001时，1-5年的平均值、2-6年的平均值。

二维码

扫码加我拉你入群

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

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

2014-11-10 09:11:20

lizhewenbei 发表于 2014-11-10 09:03
朱蕾的程序很好，但是如果能考虑stkcd就更好了。因为要针对每一个公司来做，比如stkcd=000001时，1-5年的 ...

对每个stkcd 就把
proc sort data=test out=test1;
by year;
run;
中的year 改为 stkcd

二维码

扫码加我拉你入群

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

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

2014-11-10 09:18:39

proc import out=test
datafile="d:\test1.csv"
dbms=csv replace;
getnames=yes;
run;
data a;
set test;
by stkcd year;
if last.year ;
sum+earning;
sum_1=sum-earning;
run;
data b;
merge a a(firstobs=3 drop=sum_1 earning rename=(sum=sum_2)) ;
ma=(sum_2-sum_1);
run;

二维码

扫码加我拉你入群

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

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

2014-11-10 09:43:57

朱蕾Julie 发表于 2014-11-10 09:11
对每个stkcd 就把
proc sort data=test out=test1;
by year;

按照您原来的程序（按照year排序的），好像只能得到stkcd=900957的相关数据（求和、移动平均）。我的意思是，能否把其他公司的计算结果也显示出来？感谢您！

二维码

扫码加我拉你入群

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

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

2014-11-10 09:43:58

朱蕾Julie 发表于 2014-11-10 09:11
对每个stkcd 就把
proc sort data=test out=test1;
by year;

按照您原来的程序（按照year排序的），好像只能得到stkcd=900957的相关数据（求和、移动平均）。我的意思是，能否把其他公司的计算结果也显示出来？感谢您！

二维码

扫码加我拉你入群

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

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

2014-11-10 09:53:40

那个stkcd=900957求的是900957个stkcd的总和
后面那个程序我写错了不好意思啊

二维码

扫码加我拉你入群

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

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

2014-11-10 10:07:07

朱蕾Julie 发表于 2014-11-10 09:53
那个stkcd=900957求的是900957个stkcd的总和
后面那个程序我写错了不好意思啊

proc import out=test
datafile="d:\test1.csv"
dbms=csv replace;
getnames=yes;
run;
data a;
set test;
by stkcd year;
if last.year ;
sum+earning;
sum_1=sum-earning;
if first.stkcd then do sum=earning ;sum_1=0;
end;
run;
data b;
do obsnum=3 to last;
set a point=obsnum ;
by stkcd year;
if last.year;
end;
run;
data c;
merge a b(drop=sum_1 earning rename=(sum=sum_2)) ;
ma=(sum_2-sum_1);
run;

二维码

扫码加我拉你入群

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

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

2014-11-10 10:10:27

现在行不？

二维码

扫码加我拉你入群

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

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

2014-11-10 10:13:51

朱蕾Julie 发表于 2014-11-10 09:11
对每个stkcd 就把
proc sort data=test out=test1;
by year;

按照您原来的程序（按照year排序的），好像只能得到stkcd=900957的相关数据（求和、移动平均）。我的意思是，能否把其他公司的计算结果也显示出来？感谢您！

二维码

扫码加我拉你入群

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

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

2014-11-10 10:13:57

我试过没错了

二维码

扫码加我拉你入群

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

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

2014-11-10 10:15:07

朱蕾Julie 发表于 2014-11-10 10:13
我试过没错了

多谢您，由于系统错误，消息发送了很多次，真抱歉！

二维码

扫码加我拉你入群

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

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

2014-11-10 10:17:00

lizhewenbei 发表于 2014-11-10 10:15
多谢您，由于系统错误，消息发送了很多次，真抱歉！

proc import out=test
datafile="d:\test1.csv"
dbms=csv replace;
getnames=yes;
run;
data a;
set test;
by stkcd year;
if last.year ;
sum+earning;
sum_1=sum-earning;
if first.stkcd then do sum=earning ;sum_1=0;
end;
run;
data b;
do obsnum=3 to last;
set a point=obsnum ;
by stkcd year;
if last.year;
end;
run;
data c;
merge a b(drop=sum_1 earning rename=(sum=sum_2)) ;
ma=(sum_2-sum_1);
run;
没事
这是新程序我试过没错了

二维码

扫码加我拉你入群

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

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

2014-11-10 10:21:43

朱蕾Julie 发表于 2014-11-10 10:17
proc import out=test
datafile="d:\test1.csv"
dbms=csv replace;

万谢。还有一个日志报错，这是怎么回事呢？
854  data b;
855  do obsnum=3 to last;
856  set a point=obsnum ;
857  by stkcd year;
858  if last.year;
859  end;
860  run;

NOTE: 变量 last 未初始化。
ERROR: POINT= 选项与 BY 语句不兼容。
NOTE: 由于出错，SAS 系统停止处理该步。
WARNING: 数据集 WORK.B 可能不完整。该步停止时，共有 0 个观测和 6 个变量。
WARNING: 数据集 WORK.B 由于该步已停止，而没有被替换。
NOTE: “DATA 语句”所用时间（总处理时间）:
   实际时间       0.00 秒
   CPU 时间       0.01 秒

二维码

扫码加我拉你入群

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

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

2014-11-10 10:29:01

不好意思是我发错了等一下

二维码

扫码加我拉你入群

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

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

2014-11-10 10:40:13

proc import out=test
datafile="d:\test1.csv"
dbms=csv replace;
getnames=yes;
run;
data a;
set test;
by stkcd year;
if last.year ;
sum+earning;
sum_1=sum-earning;
if first.stkcd then do sum=earning ;sum_1=0;
end;
run;
data b;
set a (firstobs=3 );
by stkcd year;
if last.year;
run;
data c;
merge a b(drop=sum_1 earning rename=(sum=sum_2)) ;
ma=(sum_2-sum_1);
run;

二维码

扫码加我拉你入群

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

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

2014-11-10 10:50:23

朱蕾Julie 发表于 2014-11-10 10:40
proc import out=test
datafile="d:\test1.csv"
dbms=csv replace;

感谢朱蕾，第一家公司（stkcd=1）的ma计算非常正确，但是从第二家公司（stkcd=2）的ma计算就有些问题了，还请您指点一下

二维码

扫码加我拉你入群

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

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

2014-11-10 10:52:01

lizhewenbei 发表于 2014-11-10 10:50
感谢朱蕾，第一家公司（stkcd=1）的ma计算非常正确，但是从第二家公司（stkcd=2）的ma计算就有些问题了， ...

第二家的也是从2003年那行开始看的

二维码

扫码加我拉你入群

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

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

2014-11-10 10:54:55

朱蕾Julie 发表于 2014-11-10 10:52
第二家的也是从2003年那行开始看的

看明白了！也就是2001和2002年的可能无法计算，要有一个公认的开始年份（2003）。这个程序已经很棒了！感谢你一直陪着！

二维码

扫码加我拉你入群

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

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

2014-11-10 10:58:43

不用谢我弄得很挫耽误你时间了

二维码

扫码加我拉你入群

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

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

2014-11-10 11:04:17

朱蕾Julie 发表于 2014-11-10 10:58
不用谢我弄得很挫耽误你时间了

相当厉害了！比我强很多！

有机会当面感谢，我就在人大。

二维码

扫码加我拉你入群

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

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

2016-2-5 10:57:30

非常感谢~

二维码

扫码加我拉你入群

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

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

2018-4-2 15:56:56

斤斤计较你发表于 2014-11-9 22:38
参考一下：

data class;

有用~~

二维码

扫码加我拉你入群

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

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

相关推荐

栏目导航

热门文章

推荐文章

扫码加好友，拉您进群

各岗位、行业、专业交流群