SAS中文论坛

标题: 求助高人:样本间隔的统计 [打印本页]

作者: shiyiming 时间: 2009-3-22 00:07
标题: 求助高人:样本间隔的统计
求助高人:样本间隔的统计
下列程序是想统计1出现时,距上次出现间隔了多少次,怎么不对,请高手帮助?
data aa;
input x @@;
cards;
1 2 1 1 2 2 1 1 1 1 2 2 2 2 1 2 2 1 2 2 2 1
;
run;
data aa1;
set aa;
if x=1  then x_1+1; else x_1=.;
if x^=1 then x_2+1; else x_2=.;
if x=1 and x^=lag(x) then x_1jg=lag(x_2);
if x=1 and x =lag(x) then x_1jg=0;
if x=2 and x^=lag(x) then x_2jg=lag(x_1);
if x=2 and x =lag(x) then x_2jg=0;
run;

我想得到下列形式:
x  x_1jg  x_2jg
1    0       .
2    .       1
1    1       .
1    0       .
2    .       2
2    .       0
1    2       .
1    0       .
1    0       .
1    0       .
2    .       4
2    .       0
2    .       0
2    .       0
1    4       .
2    .       1
2    .       0
1    2    .
2    .       1
2    .       0
2    .       0
1    3       .

谢谢!!

作者: shiyiming 时间: 2009-3-23 16:18
标题: Re: 求助高人:样本间隔的统计
data t;
set aa;
retain sum 0;
if x = 1 then
sum = 0;
else sum + 1;
run;

作者: shiyiming 时间: 2009-3-23 21:11
标题: Re: 求助高人:样本间隔的统计
谢谢指教!! 好象不对?

解释一下：变量x是样本，变量 x_1 是样本1连续出现的个数，变量 x_2 是样本2连续出现的个数
变量 x_1jg 是样本1这次出现距上次出现的间隔期数，变量 x_2jg 是样本2这次出现距上次出现的间隔期数。

请教？样本1这次出现距上次出现的间隔期数x_1jg ，不就是等于样本2连续出现的个数。
用 if x=1 and x^=lag(x) then x_1jg=lag(x_2); 这一句，想得到样本1这次出现距上次出现的间隔期数，这句为什么不对呢？
谢谢！！！

我想得到下列形式:
x x_1 x_2 x_1jg x_2jg
1 1 . 0 .
2 . 1 . 1
1 1 . 1 .
1 2 . 0 .
2 . 1 . 2
2 . 1 . 0
1 1 . 2 .
1 2 . 0 .
1 3 . 0 .
1 4 . 0 .
2 . 1 . 4 /*这里的4表示：样本2的这次出现，距上次出现间隔了4期*/
2 . 2 . 0 /*这里的0表示：样本2连续出现，距上次出现间隔了0期*/
2 . 3 . 0
2 . 4 . 0
1 1 . 4 .
2 . 1 . 1
2 . 2 . 0
1 1 . 2 .
2 . 1 . 1
2 . 2 . 0
2 . 3 . 0
1 1 . 3 .

作者: shiyiming 时间: 2009-3-24 12:23
标题: Re: 求助高人:样本间隔的统计
[code:29vm67q9]
data a;
input x @@;
cards;
1 2 1 1 2 2 1 1 1 1 2 2 2 2 1 2 2 1 2 2 2 1
;
data a;
set a(keep=x);
if x=1 then x_1+1;
  else x_1 =.;
if x=2 then x_2+1;
  else x_2 = .;
  if x=1 and _n_=1 then x_1jg=0;
  if x=2 and _n_=1 then x_2jg=0;
if x=1 and lag(x)=1 then x_1jg=0;
  else if x=1 and lag(x)=2 then do;
            x_1jg=u;
            u=0;
         end;
else if x=2 then u+1;
if x=2 and lag(x)=2 then x_2jg=0;
  else if x=2 and lag(x)=1 then do;
            x_2jg=t;
            t=0;
         end;
else if x=1 then t+1;
drop u t;
run;
[/code:29vm67q9]
大哥~~~，就别要求这么多了，我也就想着最简单的方法，那就是一样的逻辑~~~
要么我就换成array来处理了~~~~
好像要简单一些~~~
只要解决问题就行了~~~:D:D..(如果数据量不是那么大~~），
况且if else if 的效率还是比较高的。。

作者: shiyiming 时间: 2009-3-24 16:38
标题: Re: 求助高人:样本间隔的统计
对于第一次分别出现1和2的时候，你的间隔算法是什么，这个没写清楚。
按照你给出的结果来推理你的第一次处理逻辑，那么实现代码不会超过10行。
so easy

欢迎光临 SAS中文论坛 (https://mysas.net/forum/)