求助：三组数据分析的方法？

shiyiming

准备按1期，2期，3期，三组数据分析生存率，查了教材，都只有两组数据分析的？
请教一下，三组数据分析的方法，使用cox半参数分析吗？如何编程呢？

shiyiming

何谓一期二期三期？说清楚点，大家可以帮着想想

shiyiming

谢谢！
肿瘤病人的分期，例如1期病人的生存时间分别是10月，12月，13月，13月，15月；2期病人的生存时间分别是7月，8月，10月，12月，12月；3期病人生存时间分别是3月，3月，4月，5月，6月。
如何进行生存分析？
我是新手，对统计不太了解。

shiyiming

to lh605605605

所以此处的期即为病人的病情？3期最严重？我推测你要比较肿瘤病人分期的生存状况。
If not to adjust for other covariates, Kaplan-Meier in proc LIFETEST can be used to plot survival curves and to give the log-rank test for the difference among the survivals. Otherwise, for multivariable model, COX in proc PHREG can be used to compare the differences among groups and give the hazard ratios between groups, after adjusting for other cofounders. Of course, you need to know some concepts, such as censor, event, survival, hazard, and so on. I would recommend you reading a book by Paul Allison. This is an old book (in 1998) but very instructive on survival analysis using SAS.

The sample code should be like:

[code:1w9i7dj3]proc lifetest plots=(survival);
   cancer: time time*censor(0);
                 strata group;*group in (1 2 3);
run;

proc phreg;
   class group;
   model time*censor(0) = group/rl;
run;[/code:1w9i7dj3]

JingJu

shiyiming

谢谢！
我还是不太明白，能把程序写完整吗？

shiyiming

to lh605605605
Give me a sample data

shiyiming

谢谢jingju11!
简单举个例子，癌症病人分期为1，2，3期，统计了生存时间，要进行生存分析。
1期病人的生存时间分别是10月，12月，13月，13月，15月；
2期病人的生存时间分别是7月，8月，10月，12月，12月；
3期病人生存时间分别是3月，3月，4月，5月，6月

shiyiming

to lh
J 给你写出程序你还没弄懂,说明你目前不适合分析这样的数据
要么先啃书,要么请别人帮你分析为好.

shiyiming

说的对。
我的理解是，三组病人每组五个。观察窗口是从确诊为几期癌症到病人死亡。值得提醒的是，1.没有死亡的病人也有价值。如果存在，必须包括在数据之中2.三期病人比较一二两组，生存时间都短。所以很可能估计出无穷小或者是无穷大的hazard几率来。

shiyiming

我的确是菜鸟，所以要请教各位。我说的详细点，请各位帮帮忙！
数据：
8   1 1 1
52  1 1 0
58  1 1 1
63  1 1 1
63  1 1 1
220 1 1 0
365 0 1 0
452 0 1 0
496 0 1 0
528 0 1 0
560 0 1 0
676 0 1 0
13  1 2 1
18  1 2 1
23  1 2 1
70  1 2 0
76  1 2 0
180 1 2 0
195 1 2 0
210 1 2 0
232 1 2 0
300 1 2 0
396 1 2 0
490 0 2 0
540 0 2 0
11  1 3 1
13  1 3 0
13  1 3 1
16  1 3 1
23  1 3 0
33  1 3 1
43  1 3 0
67  1 3 0
87  1 3 0
111 1 3 0
123 1 3 0
156 1 3 0

程序是：data lifereg2;
infile 'd:\lifereg2.dat';
input days censor group renal @@;
proc phreg;
model days*censor(0)=group/rl;
run;
请问这样对吗？
运行结果：                                                           SAS 系统                           


                                           Analysis of Maximum Likelihood Estimates

                             Parameter      Standard                                  Hazard      95% Hazard Ratio
          Variable    DF      Estimate         Error    Chi-Square    Pr > ChiSq       Ratio      Confidence Limits

          group        1       0.92529       0.28834       10.2978        0.0013       2.523       1.434       4.439

请问结果对吗？