请教高手如何使得我的程序简化, The shorter the better!

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

立即领取

感谢您参与论坛问题回答

经管之家送您两个论坛币！

+2 论坛币

ask.txt (2.29 KB)

2013-11-4 15:14:47 上传

Original SAS Code

Please help me to simplify the code above, shoter is better!
Thank you!

扫码加我 拉你入群

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

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

分享0 收藏0 回帖

关键词：shorter better Short 请教高手 The shorter better sigma power

1. proc fcmp outlib=work.funcs.test;
   
2. subroutine pow(alpha, sigma, c1, c2, power0, power, n);
   
3.         outargs n, power;
   
4.         do n=1 by 1 until(power>=power0);
   
5.                 Z_alpha = quantile('t', 1-alpha, n);
   
6.                 Z1 = Z_alpha - (c1-c2)*sqrt(n)/sigma;
   
7.                 Z2 = Z_alpha - (c1+c2)*sqrt(n)/sigma;
   
8.                 power = 1 - (probt(Z1,n) + probt(Z2,n));
   
9.         end;
   
10. endsub;
    
11. run;
    
12. options cmplib=(work.funcs);
    
13. data final;
    
14.         alpha=0.05;
    
15.         do sigma=0.29, 0.34;
    
16.                 do lower=0.7, 0.8;
    
17.                         do ratio=0.9, 0.95, 1;
    
18.                                 do power0=0.8, 0.9;
    
19.                                         call pow(alpha, sigma, -log(lower), -log(ratio), power0, power, n);
    
20.                                         output;
    
21.                                 end;
    
22.                         end;
    
23.                 end;            
    
24.         end;
    
25.         stop;
    
26.         label   
    
27.                 lower = 'Low Bound'
    
28.                 ratio = 'True Ratio'
    
29.                 ;
    
30.         format power 9.7 ratio 9.2;
    
31. run;

复制代码