关于生存数据的方法探讨，重大发现

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

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

立即领取

感谢您参与论坛问题回答

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

+2 论坛币

我用SAS做了一下生存时间数据的模拟，用两种方法，一是Log-rank,一是基于秩次的kruskal wallis h test ,结果前者的type I error 显著大于后者，这是为什么？？难道是后一种方法更好？？

扫码加我 拉你入群

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

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

分享0 收藏0 回帖

关键词：Kruskal Wallis Error Wall type 数据 探讨 生存

发模型和程序探讨。。。

标题党呀

程序如下：
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
我不是标题党

知道原因了，生存分析因针对远期和近期差别，有多种方法，选择合适的方法分析是提高生存分析效率的关键！

7# previal
我按照你的思路复制检测。的确如你所言，二者的p不尽相同。我的问题是，为什么你认为P不相同，一种结果便优于另外一个呢？

1. 
   
2. /*create data*/
   
3. data aa;
   
4. do n = 1 to 1000;
   
5.   do group = 1 to 3;
   
6.    do _n_ = 1 to 5;
   
7.     status=round(ranuni(20100628));   
   
8.     surtime = min(50, ranexp(20100628) / 0.02);  
   
9.     if surtime =50 then status = 0;
   
10.     output;
    
11.    end;
    
12.   end;
    
13. end;
    
14. run;  
    
15. /*get two tests*/
    
16. proc npar1way data=aa WILCOXON noprint;
    
17. by n;
    
18. class group;                                                                                                                           
    
19. var surtime;                                                                                                                           
    
20. output out=wilcoxon1 (keep=P_KW );  
    
21. run;
    
22. ods output homtests=life1(keep= test ProbChiSq where = (Test = 'Log-Rank'));                                                                                                            
    
23. proc lifetest data=aa method=lt;  
    
24. by n;
    
25. time surtime * status(0);                                                                                                               
    
26. strata group;   
    
27. run;
    
28. ods output clear;         
    
29. data all;
    
30. set life1(rename =(ProbChisq = p)) wilcoxon1(rename =(p_kw = p) in = w) ;
    
31. if w then test = 'Wilcoxon';
    
32. run;
    
33. /*count significant p values*/
    
34. proc format;
    
35. value Pfmt low-0.05 ='p <=0.05' 0.05<-high = 'p>0.05';
    
36. run;
    
37. proc freq data = all;
    
38. tables test*p/chisq;
    
39. format p Pfmt.;
    
40. run;

复制代码

8# jingju11
是计算的一类错误，即1000次中P<=0.05的个数，不过我现在又进行了修改，考虑了生存分析各种方法的结合，选择最优的结果

9# previal
前辈的程序好简洁，我还要继续努力！