SAS实现生存分析之ProbPlot

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

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

立即领取

感谢您参与论坛问题回答

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

+2 论坛币

       你也许疑惑在某些SAS提供的生存分析过程里,比如说PROC LIFEREG, PROC RELIABILITY (QCpackage), 如何绘制和阅读概率分布图.基本上讲,该图通过比较实际概率分布和理论概率分布,进而得出拟和的相似程度,以判断模型的优劣.所谓的理论分布是你给出的参数模型在某些应变量条件上的取值.具体我将通过我的一片小文进行阐述.

In the following examplethe data is from Example 50.5Probability Plotting—Right Censoring.

      SAS提供不同的方法来估计实际的分布函数,但是方法和结果大同小异. 以下展示的方法为Expected Rank Plotting Position(mi).而对mi进行变换将得到PPLOT的Y轴.X的取值可以是时间或者log(时间).而与其对比的理论分布直线来源于模型分布本身.对于伪不分布,Y =log(1/(1-cdf)).

      cdf是分布参数的函树.其中更加详细的计算我将在我blog里给出.如果计算合适,得出的结果应该和SAS给出的吻合一致.

1. datat;
   
2.       set samp.fan;
   
3.       time=lifetime;
   
4.       censored=censor;
   
5. proc sort ; by time censored;
   
6. run;
   
7. /*    ExpectedRank Plotting Position*/
   
8. Data PPOS_exprank;
   
9.       set tnobs=ri end=Eof;
   
10.       by time;
    
11.       retain r 1;
    
12.       if censored=0 then do;
    
13.             c=ri/(ri +1);
    
14.             r =r*c;
    
15.             ai_exprank=1-r;
    
16.       end;
    
17.       ir =ri;
    
18.       ri ++(-1);
    
19.       if _n_=1 then do;
    
20.             put "Table 1.1 Expected Rank PlottingPositionCalculations";
    
21.             put 72*"-";
    
22.             put @1 "Ordered Obs" @15 "ReverseRank" @30 "ri/(ri+1)" @44 "Ri" @53 "ai=1-Ri" @65 "mi=G^-1(ai)";
    
23.             put 72*"-";
    
24.       end;
    
25.       timec=cats(time,ifc(censored =1,'+',''));
    
26.       rc=ifn(censored =1, .,r);
    
27.       mi=log(1/(1-ai_exprank));
    
28.       put @1 timec $8.-R @15 ir best8. @30 c fract9.@44 rc f9.4-L @53 ai_exprank f9.4-L @65 mi f9.4-L;
    
29.       if Eof then do;
    
30.             put 72*"-";
    
31.             put "+ Censored Times";
    
32.       end;
    
33. run;

复制代码

京剧

扫码加我 拉你入群

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

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

分享0 收藏4 回帖

关键词：ProbPlot 生存分析 plot bpl Rob 生存分析 ProbPlot sas

感谢分享，支持楼主

好贴，谢谢分享。