求助高人:样本间隔的统计

shiyiming

求助高人:样本间隔的统计
   下列程序是想统计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

data t;
   set aa;
   retain sum 0;
if x = 1 then
   sum = 0;
   else sum + 1;
run;

shiyiming

谢谢指教!! 好象不对?

解释一下：变量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

[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

对于第一次分别出现1和2的时候，你的间隔算法是什么，这个没写清楚。
按照你给出的结果来推理你的第一次处理逻辑，那么实现代码不会超过10行。
so easy