以下的程序怎么编程

shiyiming

即使要追求简洁，也要保证以正确为前提。
sXLION的程序有以下欠考虑１。处理缺失值，２。K＝０，３。K＞＝N
没有任何不敬。

shiyiming

经常看SAS-L和SUGI，都是那里边出来的。。
SAS-L上的大师确实太牛~~~
monotonic()跟DB里的rownum一样，但是功能现在还不健全，所以没有列出来。
proc slq还有很多其他的options比较强大。。不过都没有结束测试，还没有上reference..

shiyiming

to deshengsu
你可以把里面的数值做缺失值看看。
另考虑k=0和k>=N，您认为有意义吗？

shiyiming

data result (keep=mean_a_dif mean_b_dif);
     retain mean_a mean_b;
    set ahuige end=last;
    total_a+a;
    if not missing(a) then num_a+1;
    total_b +b;
    if not missing(b) then num_b+1;
    if _N_=&k. then do;
       mean_a=total_A/num_a;
       mean_b=total_b/num_b;
       total_a=0;
       num_a=0;
       total_b=0;
       num_b=0;
    end;
    if last then do;
       mean_a_dif=mean_a-total_a/num_a;
       mean_b_dif=mean_b-total_b/num_b;
       output;
    end;
run;

shiyiming

[quote="deshengsu":3que81hy]即使要追求简洁，也要保证以正确为前提。
sXLION的程序有以下欠考虑１。处理缺失值，２。K＝０，３。K＞＝N
没有任何不敬。[/quote:3que81hy]

不好意思，才明白你的意思。

不将缺失值计算在内，mean=总和/项数（不包含缺失项）                  

如果是2的话：

[code:3que81hy]data ahuige;
  input A B;
cards;
. 3
2 3
5 7
7 9
;
run;

%let k=2;
data ex(keep=ma mb);
   set ahuige nobs=n;
  if &k. ge _N_ then
      do;
          aa1+A;bb1+B;
           if missing(A) then missa1+1;
                   if missing(B) then missb1+1;
        end;
      else do;
           aa2+A;bb2+B;
                    if missing(A) then missa2+1;
                    if missing(B) then missb2+1;
            end;

  if _N_=n then
        do;
           ma=aa1/(&k.-missa1)- aa2/(n-&k.-missb1);
           mb=bb1/(&k.-missa2)- bb2/(n-&k.-missb2);
           output;
        end;
   run;[/code:3que81hy]


ps:不好意思，记错了。

shiyiming

1，将缺失项计算在内，mean=总和/项数（包含缺失项） 这种算法是SAS函数mean（）默认的。

shiyiming

data trial;
y=mean(1,.,3);
;
proc print;run;

shiyiming

[code:26prqt1b]data ahuige;
  input A B;
cards;
1 2
1 2
. 3
2 3
5 7
7 9
;
run;
%let k=2;
proc sql;
select mean(a), mean(b),&k+1 into :ma1, :mb1, :k1 from ahuige(obs=&k.);
select &ma1.- mean(a),&mb1.- mean(b)  into :ma2,:mb2 from ahuige(firstobs = &k1. );
quit;[/code:26prqt1b]

shiyiming

[code:3uuadprv]data ahuige;
  input A B;
cards;
. 3
2 3
5 7
7 9
;
run;
%let k=2;

proc sql ;
select mean(a), mean(b) into :ma1, :mb1 from ahuige where monotonic() le &k.;
select &ma1.-mean(a),&mb1.-mean(b) into :ma,:mb from ahuige where monotonic() gt &k.;
quit;[/code:3uuadprv]

shiyiming

真的非常佩服楼上的各位大师！可以问问你们是怎么找出这么多种巧妙脱俗的方法地呢？你们可以进一步谈谈解决这个问题的思路及感想吗？莫非冥冥之中真有注定？
这里真是一个好地方！下次有问题一定跪求！