SAS中数组的应用

jingju的程序就是这样算的

jingju11 发表于 2012-6-13 21:08
Interesting. How? Jingju

呵呵，搞复杂了。还是倒着计0方便。

jingju11 发表于 2012-6-13 21:08
Interesting. How? Jingju

试了下，如下：不是很严谨。谨以楼主提供的部分数据为例：

1. data raw;
   
2. input id  h0 h0_5 h1 h2 h3 h4 h5 h6;
   
3. datalines;
   
4. 1        6        4        3        1        1        0        0        0
   
5. 2        7        6        5        3        1        0        0        0
   
6. 3        10        10        10        10        .        .        .        .
   
7. 4        7        6        4        3        3        2        2        0
   
8. 5        8        3        0        0        1        1        1        1
   
9. 6        8        6        .        .        .        .        .        .
   
10. 7        8        6        3        2        2        2        2        2
    
11. 8        6        1        0        0        0        0        0        0
    
12. 9        6        3        1        1        0        0        0        0
    
13. 10        8        5        2        0        1        3        4        6
    
14. 11        7        5        3        1        0        0        0        1
    
15. 12        7        5        3        2        0        2        3        4
    
16. 13        8        6        4        4        2        3        2        2
    
17. 14        7        3        0        0        0        0        0        7
    
18. 15        7        6        4        3        1        0        0        0
    
19. 16        8        6        1        0        0        0        0        0
    
20. 17        8        5        3        1        1        1        0        0
    
21. 18        5        4        3        0        0        0        0        0
    
22. 19        5        3        3        1        0        0        0        0
    
23. 20        7        6        5        5        4        4        3        3
    
24. ;
    
25. 
    
26. data want (drop=i j y sum missing );
    
27. array t[*]         h0_5 h1 h2 h3 h4 h5 h6;
    
28.   set raw;
    
29. do i=1 to dim(t) until (sum=0);
    
30.    do; if missing(t[i]) then missing=1;         else missing=0;
    
31.    sum=0;
    
32.    end;
    
33.   do j=i to dim(t) ;
    
34.       sum+t[j];
    
35.   end;
    
36.   y=vname(t[i]);
    
37. if sum^=0 then time="未观察到疼痛消失时间";
    
38. else if missing=0  then time=tranwrd(substr(y,2,length(y)-1),"_",".");
    
39.        else time="疼痛纪录有缺失";
    
40. end;
    
41. run;

复制代码

result:

Obs    h0_5    h1    h2    h3    h4    h5    h6    id    h0    time

           1      4      3     1     1     0     0     0     1     6    4
           2      6      5     3     1     0     0     0     2     7    4
           3     10     10    10     .     .     .     .     3    10    疼痛纪录有缺失
           4      6      4     3     3     2     2     0     4     7    6
           5      3      0     0     1     1     1     1     5     8    未观察到疼痛消失时间
           6      6      .     .     .     .     .     .     6     8    疼痛纪录有缺失
           7      6      3     2     2     2     2     2     7     8    未观察到疼痛消失时间
           8      1      0     0     0     0     0     0     8     6    1
           9      3      1     1     0     0     0     0     9     6    3
          10      5      2     0     1     3     4     6    10     8    未观察到疼痛消失时间
          11      5      3     1     0     0     0     1    11     7    未观察到疼痛消失时间
          12      5      3     2     0     2     3     4    12     7    未观察到疼痛消失时间
          13      6      4     4     2     3     2     2    13     8    未观察到疼痛消失时间
          14      3      0     0     0     0     0     7    14     7    未观察到疼痛消失时间
          15      6      4     3     1     0     0     0    15     7    4
          16      6      1     0     0     0     0     0    16     8    2
          17      5      3     1     1     1     0     0    17     8    5
          18      4      3     0     0     0     0     0    18     5    2
          19      3      3     1     0     0     0     0    19     5    3
          20      6      5     5     4     4     3     3    20     7    未观察到疼痛消失时间

So good.

如果一组数据中同时存在正负数，怎么求这组数的总和呢

学习了

data raw;
input id h0 h0_5 h1-h6;
datalines;
1 6 4 3 1 1 0 0 0
2 7 6 5 3 1 0 0 0
3 10 10 10 10 . . . .
4 7 6 4 3 3 2 2 0
5 8 3 0 0 1 1 1 1
6 8 6 . . . . . .
7 8 6 3 2 2 2 2 2
8 6 1 0 0 0 0 0 0
9 6 3 1 1 0 0 0 0
10 8 5 2 0 1 3 4 6
11 7 5 3 1 0 0 0 1
12 7 5 3 2 0 2 3 4
13 8 6 4 4 2 3 2 2
14 7 3 0 0 0 0 0 7
15 7 6 4 3 1 0 0 0
16 8 6 1 0 0 0 0 0
17 8 5 3 1 1 1 0 0
18 5 4 3 0 0 0 0 0
19 5 3 3 1 0 0 0 0
20 7 6 5 5 4 4 3 3
;
run;
data want1(drop=i );
        array t[*] h0_5 h1 h2 h3 h4 h5 h6;
          do n=1 to total;
  set raw point=n nobs=total;
                do i=1 to dim(t)  until (t[i]=0  or missing=1 or i=7);
                        if missing(t[i]) then missing=1;
                                else do missing=0; sum+t[i]; end;
                        put n= i=  missing= sum= t[i]=;
                end;
        y=vname(t[i]);
                put  y=;
        if missing=1 then time="疼痛纪录有缺失";
      else if t[i]=0 then time=tranwrd(substr(y,2),"_",".");
                                else time="未观察到疼痛消失时间";                       
output;
  end;
  stop;
  run;

panglinqi 发表于 2019-10-11 15:02
data raw;
input id h0 h0_5 h1-h6;
datalines;

要注意两个点
一是：（如下），确保逐行运行，不多不少。
do n=1 to total;
  set raw point=n nobs=total;
output;
  end;
二是：明确该循环的终止条件，以及每个条件对应的time值类别。
do i=1 to dim(t)  until (t=0  or missing=1 or i=7);