Kaplan-Meierプロット・Forestプロット作成の応用： グラフ …

1 Kaplan-Meier Forest 1 2 2 3 1 2 3 Advancement in both Kaplan-Meier and forest plots: quantitative results output inside or outside the graph area Ryuji Uozumi1 , Saori Yoshida2 , Takayuki Hirai2 , and Chikuma Hamada3 1 Department of Biomedical Statistics and Bioinformatics, Kyoto University Graduate School of Medicine 2 Biostatistics Team, Pharmaceutical Development Division, Nippon Kayaku Co.,Ltd 3 Department of Information and Computer Technology, Faculty of Engineering, Tokyo University of Science Kaplan-Meier Forest ODS GRAPHICS Kaplan-Meier Forest Kaplan-Meier LIFETEST SGPLOT XAXISTABLE Forest SGPLOT Forest Kaplan-Meier Forest SGPLOT XAXISTABLE YAXISTABLE X2 AXIS Y2 AXIS 1 SAS ODS (Output Delivery System)

2 GRAPHICS ODS GRAPHICS ON LIFETEST Kaplan-Meier SG (Statistical Graphics) SGPLOT SGPANEL SAS SAS (2013) GPLOT SGPLOT SGPLOT (2014) SGPLOT (2012) SG Annotation SG GTL (Graph Template Language) ODS PDF SAS HTML CDISC ADaM (Clinical Data Interchange Standards Consortium Analysis Data Model) ( , 2011, 2012) ODS PowerPoint SAS PowerPoint ( , 2015) Kaplan-Meier 20 ( , 1995) SAS Kaplan-Meier (2010) LIFETEST GPLOT Kaplan-Meier (2011) SGPLOT SGRENDER Kaplan-Meier LIFETEST SGPLOT ( , 2015) (2010) (2011) Kaplan-Meier (number of subjects at risk)

3 LIFETEST ODS GRAPHICS Kaplan-Meier LIFETEST ATRISK Kaplan-Meier Kaplan-Meier Kaplan-Meier Forest Forest SAS Kaplan-Meier Forest LIFETEST SGPLOT ODS "STYLE = JOURNAL" 2 Kaplan-Meier Kaplan-Meier Kaplan-Meier VALung VALung (Veterans Administration Lung cancer trial) Kalbfleisch and Prentice (2002) (n = 137) SAS/STAT phreg (2016) 3 ( ) ( : Time) 1 ( : Cell) 4 (Cell = 'adeno' 'small' 'large' 'squamous' ) Kaplan-Meier LIFETEST PLOTS = S ATRISK 1 (2016) SAS proc lifetest data=VALung plots=s(atrisk atrisktick); time Time*Censor(1); strata Cell; run; 1 : LIFETEST Kaplan-Meier 1 SAS ods graphics.

4 Ods output Survivalplot=Survivalplot; proc lifetest data=VALung atrisk plots=s(atrisk); time Time*Censor(1); strata Cell; run; *------------------------------------*; *-- ( A ) --*; *------------------------------------*; proc sgplot data=Survivalplot2 noautolegend; step x=Time y=Survival / group=StratumNum name='Survival'; scatter x=Time y=Censored / group=StratumNum markerattrs=(symbol=plus); scatter x=tAtRisk y=StratumNum / markerchar=atrisk y2axis; keylegend 'Survival' / location=outside noborder position=bottom; yaxis offsetmin= values=(0 to 1 by ) label='Survival'; y2axis offsetmin= offsetmax= min=1 max=4 label=' ' display=(noticks novalues); xaxis offsetmin= values=(0 to &TimeMax by &TimeInterval) label='Time'; format StratumNum Cellf.; run; 2 : SGPLOT Kaplan-Meier 1 1 VALung 0 0 0 0 0 0 PLOTS = S 2 SGPLOT (2016)

5 2 Y2 AXIS YAXIS OFFSETMIN = YAXIS XAXIS Y2 AXIS X2 AXIS OFFSETMIN = OFFSETMAX = 2 YAXIS OFFSETMIN = 20% 1 20% Y2 AXIS OFFSETMIN = OFFSETMAX = 85% ODS OUTPUT 0 0 SGPLOT A 1 2 Kaplan-Meier Kaplan-Meier Kaplan-Meier SGRENDER TEMPLATE LAYOUT LATTICE Kaplan-Meier ( , 2011) SGPLOT Kaplan-Meier 3 SGPLOT XAXISTABLE Kaplan-Meier XAXISTABLE AtRisk tAtRisk X = StratumNum CLASS = Kaplan-Meier XAXISTABLE YAXISTABLE XAXISTABLE LOCATION = INSIDE 3 SGPLOT Kaplan-Meier SAS/STAT LIFETEST OUTSIDE Kaplan-Meier ( 4) OUTSIDE(k)

6 ODS GRAPHICS k 4 4 14% SAS proc sgplot data=Survivalplot2 noautolegend; step x=Time y=Survival / group=StratumNum name='Survival'; scatter x=Time y=Censored / group=StratumNum markerattrs=(symbol=plus); xaxistable AtRisk / x=tAtRisk class=StratumNum; keylegend 'Survival' / location=outside noborder position=bottom; yaxis values=(0 to 1 by ) label='Survival'; xaxis values=(0 to &TimeMax by &TimeInterval) label='Time'; format StratumNum Cellf.; run; 3 : SGPLOT Kaplan-Meier 2 SAS proc lifetest data=VALung plots=s(atrisk atrisktick outside); time Time*Censor(1); strata Cell; run; 4 : LIFETEST Kaplan-Meier 2 3 Forest Forest Forest Kaplan-Meier Borghaei et al.

7 (2015) 95% Motzer et al. (2015) Turner et al. (2015) p SAS Forest (2013) GPLOT Annotation Forest Forest p VALung (Therapy2 = 1 2 ) (AgeC = 1 < 65 2 >= 65 ) (CellC = 1 2 3 4 ) (Prior = 0 1 ) 95% VALung Therapy2 AgeC (Age IF ) CellC 1 : Forest Subgroup [ ] Intpval p [ ] Type N Grp1 / (%) [ ] Grp2 / (%) [ ] HazardRatio HRLowerCL 95% HRUpperCL 95% HR (95% ) [ ] Indent Forestval Ref SGPLOT Forest phreg 1 phreg SGPLOT B 1 C 1 SGPLOT Forest 5 5 BUBBLE Forest ( , 2002)

8 XAXIS TYPE = LOG Overall Forestval PowerPoint 95% ODS PowerPoint ( , 2015) PowerPoint PowerPoint SAS SAS proc sgplot data=forestdata noautolegend; refline ref / axis=y lineattrs=(thickness=30 color=cxf0f0f0); scatter y=forestval x=HazardRatio / errorbarattrs=(thickness=1 color=blue) markerattrs=(size=0) xerrorupper=HRUpperCL xerrorlower=HRLowerCL; bubble y=forestval x=HazardRatio size=N / bradiusmax= bradiusmin= fillattrs=(color=blue) lineattrs=(color=blue); refline 1 / axis=x lineattrs=(pattern=shortdash) transparency= ; xaxis type=log minor min= max=5 display=(nolabel) valueattrs=(size=7); yaxis display=none; run; 5 : SGPLOT Forest 1 Forest Y2 AXIS 5 5 SGPLOT 6 SGPLOT Forest DATTRMAP = ATTRMAP D 5 proc sgplot data=forestdata noautolegend dattrmap=attrmap nowall noborder nocycleattrs; styleattrs axisextent=data; refline ref / axis=y lineattrs=(thickness=13 color=cxf0f0f0); scatter y=forestval x=HazardRatio / errorbarattrs=(thickness=1 color=blue) xerrorupper=HRUpperCL xerrorlower=HRLowerCL markerattrs=(size=0).

9 /* highlow y=forestval low=HRLowerCL high=HRUpperCL / lineattrs=(color=blue); */ bubble y=forestval x=HazardRatio size=N / bradiusmax= bradiusmin= fillattrs=(color=blue) lineattrs=(color=blue); scatter y=forestval x=HazardRatio / markerattrs=(size=0) x2axis; refline 1 / axis=x lineattrs=(pattern=shortdash) transparency= ; text x=xl y=forestval text=text / position=center contributeoffsets=none; yaxistable subgroup / location=inside position=left textgroup=valtype labelattrs=(size=8) labeljustify=left titlejustify=left textgroupid=text indentweight=indent; yaxistable grp1 grp2 / location=inside position=left labelattrs=(size=8) valueattrs=(size=7) valuehalign=center valuejustify=right; yaxistable HR intpval / location=inside position=right labelattrs=(size=8) valueattrs=(size=7) valuehalign=center; yaxis reverse display=none offsetmin= colorbands=odd colorbandsattrs=(transparency=1); xaxis type=log values=( 2 5) minor display=(nolabel) valueattrs=(size=7); x2axis label='Hazard Ratio (95% CI)' display=(noline noticks novalues) labelattrs=(size=8); run; 6 : SGPLOT Forest 1 ( ) 2 : Forest STYLEATTRS X ( TS1M3 ) REFLINE = 1 - SCATTER (HIGHLOW ) Hazard Ratio (95% CI) (X2 AXIS ) - BUBBLE ( ) - TEXT (<-- Test Better Standard Better -->)