cox生存曲线

做cox生存曲线 ，SAS 因为我做出来的 是0-1.0 纵轴，我想把0-0.5 部分截取，然后加上横轴，纵轴的标签。这是我的程序，向高手请教怎么修改？？  
  时间变量为 ：fabingtime_0
  删失： fabing_0(0);
  分层 b;

proc lifetest   data=aa
method=km plot=(s)  graphics lineprinter outsurv=surv;
time fabingtime_0* fabing_0(0);
strata b;
symbol1 line=1 color=blue  ;
symbol2 line=2 color=black ;
run;


附件里 有SAS导出的数据 ，请高手？？？

1. proc import datafile='H:\linwork\cox.xls'
   
2.                                 out=cox
   
3.                                 dbms=excel replace;
   
4. run;
   
5. 
   
6. /*proc lifetest data=cox */
   
7. /*method=km plot=(s)  graphics lineprinter outsurv=surv;*/
   
8. /*time fabingtime_0* fabing_0(0);*/
   
9. /*strata b;*/
   
10. /*symbol1 line=1 color=blue  ;*/
    
11. /*symbol2 line=2 color=black ; */
    
12. /*run;*/
    
13. 
    
14. goptions reset=(symbol axis);
    
15. symbol1 c=blue; symbol2 c=orange;
    
16. proc lifetest data=cox outsurv=outsurv plots=(s) noprint;
    
17.       time fabingtime_0* fabing_0(0);
    
18.       strata b;
    
19.    run;
    
20. 
    
21. /* check to see if last time is censored, and if so
    
22.    replace missing survival with that of last known
    
23.    event time */
    
24. data outsurv2;
    
25.         set outsurv;
    
26.         drop _old;
    
27.         _old=lag(survival);
    
28.         if _censor_ eq 1 and survival eq . then survival=_old;
    
29.         if fabingtime_0 ne 0 and survival ne 0 then do;
    
30.                 log_of_days=log(fabingtime_0);
    
31.                 log_negative_log_SDF = log(-log(survival));
    
32.         end;
    
33.         run;
    
34. 
    
35. 
    
36. goptions reset=(symbol axis);
    
37. symbol1 i=steplj c=blue;
    
38. symbol2 i=steplj c=orange;
    
39. symbol3 v=none c=black;
    
40. symbol4 v=none i=none c=black;
    
41. symbol5 v=circle i=none c=black;
    
42. axis1 order=(0.5 to 1 by 0.05) label=(angle=90);
    
43. proc gplot data=outsurv2;
    
44.         plot survival*fabingtime_0 = stratum / vaxis=axis1 ;
    
45.         plot2 survival*fabingtime_0=_censor_ / nolegend noaxis;
    
46.         run;
    
47.         quit;

复制代码

http://www2.sas.com/proceedings/sugi30/262-30.pdf

原图：