[quote]引用第6楼TTG于2007-09-10 22:18发表的“”:



summary(coxph(Surv(TIME,STATUS)~AGE,method="breslow",data=colon))

这样算出来，就和书上一样了，不过我想知道这个method用什么参数有什么讲究吗？[/quote]





有关 R 里的 survival analysis，推荐 Therneau 那个关于 survival package 的文档



http://www.mayo.edu/hsr/people/therneau/survival.ps



这个文档提供了几乎所有的技术细节，你甚至可以轻松写个函数，得到与coxph一样的结果，即使功能不如coxph那么全





PS. 以这个文档为基础，后又增加了一些章节，重新编辑成书，Modeling Survival Data: Extending the Cox Model



不过，就Cox的基础知识来说，那个文档已经足够
这个文档好，我最喜欢这样的文档了！
summary(coxph(Surv(TIME,STATUS)~AGE,method="breslow",data=colon))<br />
Call:<br />
coxph(formula = Surv(TIME, STATUS) ~ AGE, data = colon, method = "breslow")<br />
<br />
  n= 30 <br />
     coef exp(coef) se(coef)    z       p<br />
AGE 0.255      1.29   0.0584 4.36 1.3e-05<br />
<br />
    exp(coef) exp(-coef) lower .95 upper .95<br />
AGE      1.29      0.775      1.15      1.45<br />
<br />
Rsquare= 0.598   (max possible= 0.991 )<br />
Likelihood ratio test= 27.4  on 1 df,   p=1.68e-07<br />
Wald test            = 19.0  on 1 df,   p=1.29e-05<br />
Score (logrank) test = 8.06  on 1 df,   p=0.00451<br />


可是我的确实一样啊
[quote]引用第8楼momozilla于2007-09-11 00:55发表的“”:







有关 R 里的 survival analysis，推荐 Therneau 那个关于 survival package 的文档



.......[/quote]

非常感谢，收下好好学习。
[quote]引用第10楼TTG于2007-09-11 05:49发表的“”:

可是我的确实一样啊

.......[/quote]

哪有你这么耍赖的！你居然把顶楼的答案改了！原来的age根本不是这个数。
[quote]引用第12楼rtist于2007-09-11 19:32发表的“”:



哪有你这么耍赖的！你居然把顶楼的答案改了！原来的age根本不是这个数。[/quote]

冤枉啊，我不过是把p值的小数点后加了几位，原来的age到是有点错，我把0.25写成9.25了，可是这个不影响运算啊，那个数据部分我保证没做更改。
[quote]引用第13楼TTG于2007-09-12 03:38发表的“”:



冤枉啊，我不过是把p值的小数点后加了几位，原来的age到是有点错，我把0.25写成9.25了，可是这个不影响运算啊，那个数据部分我保证没做更改。[/quote]

谁冤枉你了？我说过你改数据了么？我清清楚楚地说你把答案给改了！！！

你自己说0.25和9.25相差几十倍？？自己提问都不花点时间把帖子检查一下，你知道别人在你这种typo上要浪费多少时间？？？难道我们的时间都不如你的值钱了？？？
[quote]引用第14楼rtist于2007-09-12 16:54发表的“”:



谁冤枉你了？我说过你改数据了么？我清清楚楚地说你把答案给改了！！！

你自己说0.25和9.25相差几十倍？？自己提问都不花点时间把帖子检查一下，你知道别人在你这种typo上要浪费多少时间？？？难道我们的时间都不如你的值钱了？？？[/quote]

写错了答案我实在非常抱歉，浪费了各位的时间我也深感不安。

不过我修改答案也是为了不再让这个浪费大家的时间，并无有意耍赖的意思。

之所以没有在帖子上再做注明，是因为我认为这和rtist兄给出的答案不一致的原因所在关系不大，对于为什么使用同样的算法而结果不同，这一点我依然不甚明白。

最后我对自己的不慎以致在一贴连续两次出错再次郑重道歉，并恳请各位尤其是rtist兄的原谅。

apology accepted.

I'm angry because I really spent A LOT of time trying to find out why the AGE coefficient was SO different.
Really sorry,and thank you for your help.
5 年 后

回复 第6楼 的 Rtist：请问怎么时依协变量来对样本做PH假设检验？我用的SPSS软件。谢谢！！！我QQ：907172721