发模型和程序探讨。。。

标题党呀

程序如下：
data a;                                                                                                                                 
do n=1 to 1000;                                                                                                                        
group=1;                                                                                                                                
do i = 1 to 5;                                                                                                                          
status=round(ranuni(20100628),1);                                                                                                      
surtime = ranexp(20100628) / 0.02;                                                                                                      
if surtime >50 then do;                                                                                                                 
surtime=50;                                                                                                                             
status=0;                                                                                                                              
end;                                                                                                                                    
output;                                                                                                                                 
end;                                                                                                                                    
end;                                                                                                                                    
run;                                                                                                                                    
                                                                                                                                       
                                                                                                                                       
data b;                                                                                                                                
do n=1 to 1000;                                                                                                                        
group=2;                                                                                                                                
do i = 1 to 5;                                                                                                                          
status=round(ranuni(20100629),1);                                                                                                      
surtime = ranexp(20100629) / 0.02;                                                                                                      
if surtime >50 then do;                                                                                                                 
surtime=50;                                                                                                                             
status=0;                                                                                                                              
end;                                                                                                                                    
output;                                                                                                                                 
end;                                                                                                                                    
end;                                                                                                                                    
run;                                                                                                                                    
                                                                                                                                       
                                                                                                                                       
data c;                                                                                                                                
do n=1 to 1000;                                                                                                                        
group=3;                                                                                                                                
do i = 1 to 5;                                                                                                                          
status=round(ranuni(20100630),1);                                                                                                      
surtime = ranexp(20100630) / 0.02;                                                                                                      
if surtime >50 then do;                                                                                                                 
surtime=50;                                                                                                                             
status=0;                                                                                                                              
end;                                                                                                                                    
output;                                                                                                                                 
end;                                                                                                                                    
end;                                                                                                                                    
run;                                                                                                                                    
                                                                                                                                       
data d;                                                                                                                                 
set a b c;                                                                                                                              
run;

接上：
/*拆分为1000个数据集*/                                                                                                                  
%macro test;                                                                                                                           
data                                                                                                                                    
  %do i=1 %to 1000;                                                                                                                     
    n&i                                                                                                                                 
  %end;                                                                                                                                 
  ;                                                                                                                                    
set d;                                                                                                                                 
%do i=1 %to 1000;                                                                                                                       
  if n=&i then output n&i;                                                                                                              
%end;                                                                                                                                   
run;                                                                                                                                    
%mend test;                                                                                                                             
%test;                                                                                                                                 
                                                                                                                                       
/*两种检验方法*/                                                                                                                        
                                                                                                                                       
%macro test1;                                                                                                                           
%do i=1 %to 1000;                                                                                                                       
/*Wilcoxon test*/                                                                                                                       
proc npar1way data=n&i WILCOXON;                                                                                                        
class group;                                                                                                                           
var surtime;                                                                                                                           
output out=wilcoxon&i (keep=P_KW );                                                                                                     
PROC APPEND BASE=wilcoxon force;                                                                                                        
run; quit;                                                                                                                              
                                                                                                                                       
/*Log-Rank test*/                                                                                                                       
ods listing close ;                                                                                                                     
ods output homtests=life&i;                                                                                                            
proc lifetest data=n&i method=lt;                                                                                                      
time surtime * status(0);                                                                                                               
strata group;                                                                                                                           
data life&i (keep= test ProbChiSq);                                                                                                     
set life&i;                                                                                                                             
if test='Log-Rank';                                                                                                                     
PROC APPEND BASE=lifetest force;                                                                                                        
run;quit;                                                                                                                              
ods listing;                                                                                                                           
%end;                                                                                                                                   
%mend test1;                                                                                                                           
%test1;                                                                                                                                 
                                                                                                                                       
/*找出P<=0.05的观测*/                                                                                                                  
data wilcoxon;                                                                                                                          
set wilcoxon;                                                                                                                           
if P_KW<=0.05 ;                                                                                                                        
run;                                                                                                                                    
                                                                                                                                       
                                                                                                                                       
data lifetest;                                                                                                                          
set lifetest;                                                                                                                           
if ProbChiSq<=0.05;                                                                                                                     
run;

3# BraveMadMan
我不是标题党