生存時間解析のサンプルデータ

統計ソフトR

http://minato.sip21c.org/medstat/it2011-12r.html

対馬(2008).SPSSで学ぶ医療系多変量データ解析

SASデータセット

転倒転落を想定した架空データを作ってみました。

Data fall;
input time censor treat;
/*
  time 転倒転落が発生するまでの日数
  censor 打ち切り:0 転倒:1
  treat  治療群:1 コントロール群:0
*/
cards;
9	1	1
13	0	1
20	1	1
31	1	1
40	0	1
48	0	1
55	0	1
55	1	1
75	1	1
120	0	1
3	1	0
3	1	0
7	0	0
7	1	0
12	0	0
15	1	0
27	0	0
33	0	0
43	0	0
45	1	0
;
PROC lifetest data=fall plots=(s, lls);
  time time*censor(1);
  strata treat;
RUN;
proc phreg simple data=fall;
  model time*censor(1) = treat;
  run;