SAS中文论坛

标题: 求助:三组数据分析的方法? [打印本页]

作者: shiyiming    时间: 2010-2-25 16:28
标题: 求助:三组数据分析的方法?
准备按1期,2期,3期,三组数据分析生存率,查了教材,都只有两组数据分析的?
请教一下,三组数据分析的方法,使用cox半参数分析吗?如何编程呢?
作者: shiyiming    时间: 2010-2-26 02:25
标题: Re: 求助:三组数据分析的方法?
何谓一期二期三期?说清楚点,大家可以帮着想想
作者: shiyiming    时间: 2010-2-26 12:42
标题: Re: 求助:三组数据分析的方法?
谢谢!
肿瘤病人的分期,例如1期病人的生存时间分别是10月,12月,13月,13月,15月;2期病人的生存时间分别是7月,8月,10月,12月,12月;3期病人生存时间分别是3月,3月,4月,5月,6月。
如何进行生存分析?
我是新手,对统计不太了解。
作者: shiyiming    时间: 2010-2-26 23:00
标题: Re: 求助:三组数据分析的方法?
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    时间: 2010-3-1 23:00
标题: Re: 求助:三组数据分析的方法?
谢谢!
我还是不太明白,能把程序写完整吗?
作者: shiyiming    时间: 2010-3-1 23:55
标题: Re: 求助:三组数据分析的方法?
to lh605605605
Give me a sample data
作者: shiyiming    时间: 2010-3-8 16:33
标题: Re: 求助:三组数据分析的方法?
谢谢jingju11!
简单举个例子,癌症病人分期为1,2,3期,统计了生存时间,要进行生存分析。
1期病人的生存时间分别是10月,12月,13月,13月,15月;
2期病人的生存时间分别是7月,8月,10月,12月,12月;
3期病人生存时间分别是3月,3月,4月,5月,6月
作者: shiyiming    时间: 2010-3-8 20:42
标题: Re: 求助:三组数据分析的方法?
to lh
J 给你写出程序你还没弄懂,说明你目前不适合分析这样的数据
要么先啃书,要么请别人帮你分析为好.
作者: shiyiming    时间: 2010-3-9 22:37
标题: Re: 求助:三组数据分析的方法?
说的对。
我的理解是,三组病人每组五个。观察窗口是从确诊为几期癌症到病人死亡。值得提醒的是,1.没有死亡的病人也有价值。如果存在,必须包括在数据之中2.三期病人比较一二两组,生存时间都短。所以很可能估计出无穷小或者是无穷大的hazard几率来。
作者: shiyiming    时间: 2010-3-23 20:33
标题: Re: 求助:三组数据分析的方法?
我的确是菜鸟,所以要请教各位。我说的详细点,请各位帮帮忙!
数据:
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

请问结果对吗?
作者: shiyiming    时间: 2010-3-24 01:27
标题: Re: 求助:三组数据分析的方法?
首先祝贺你终于成为了统计学家之分类一:生存分析学家。请不要汗颜,忆往昔,当我从字典里查到survival是什么意思的时候,我就毅然决然地在我的简历里郑重地加上一条:本人擅长统计分析,包括生存分析。。。有个锐利指出,你连测度论都不懂,还说懂生存理论?我想,党和国家领导人不懂养猪,不也关心猪农的疾苦吗?咱们闲话少说,言归正传。
你的分析大致不差,几点意见如下:
1.如何处理group变量。通常如你的情况我们把它放于class之中:CLASS GROUP;原因是1.这里不同的组之差别未必相等2.组之间应该再没有小组了吧,说白了group不能等于1.1, 1.5等等。
2.那个变量renal。无端的舍掉这个变量理由不充分。至少让它试一试,看看起不起作用。我感觉它有作用,还有大作用。
3.KM曲线。生存曲线。有一点小问题,即组一和组二相交。那个不好。也就是他们不定型。一会儿好的一会而不好的。原因也许有二:1.人太少2.组一内有个可疑值 days = 8.这未必错,但是值得检查一下,看看他是否车祸故去的,是否活了80天。
4总而言之,proportional hazard的基本假设略微有点践踏,但是说实话,我认为无大碍。
作者: shiyiming    时间: 2010-3-24 15:30
标题: Re: 求助:三组数据分析的方法?
谢谢指点!关于结果有点不明白。
Variable DF [size=150:1tafj8fm]Estimate[/size:1tafj8fm] Error Chi-Square Pr > ChiSq Ratio Confidence Limits

group 1 [size=150:1tafj8fm]0.92529 [/size:1tafj8fm]0.28834 10.2978 0.0013 2.523 1.434 4.439

Estimate是0.925,是说group1相对其他分期(group2,group3)病人的生存时间吗?应该如何理解?对于group2和group3如何分析?

我是菜鸟,可能问题比较幼稚!!!

谢谢[size=150:1tafj8fm]jingju11[/size:1tafj8fm][size=150:1tafj8fm]!!![/size:1tafj8fm]
[size=150:1tafj8fm]能麻烦帮忙把程序按你讲的改改吗?我的水平太差。[/size:1tafj8fm]
作者: shiyiming    时间: 2010-3-24 23:55
标题: Re: 求助:三组数据分析的方法?
proc phreg;
class group;
model days*censor(0)=group renal/rl;
run;
作者: shiyiming    时间: 2010-3-29 11:59
标题: Re: 求助:三组数据分析的方法?
[color=#800000:3vy0aqnp][size=200:3vy0aqnp]谢谢!!!!!!!![/size:3vy0aqnp][/color:3vy0aqnp]




欢迎光临 SAS中文论坛 (http://www.mysas.net/forum/) Powered by Discuz! X3.2