「統計学関連なんでもあり」の過去ログ---026

1213.　比例ハザード解析/R/coxph　detour　2003/12/03 (水) 15:41
└1225.　Re: 比例ハザード解析/R/coxph　detour　2003/12/04 (木) 18:40

1213.　比例ハザード解析/R/coxph　detour　　2003/12/03 (水) 15:41

虫を低温において何日で死亡するか比例ハザード解析を行ってみました。Rのsurvivalのcoxphを使いました。死亡日 (time)以外に雌雄など4要因あります。個体数は総計5500くらいです。sf<-survfit(time),data=部分集合)。 plot(sf$time,log(sf$surv))やplot(log(sf$time),log(-log(sf$surv)))，つまり，log やlog(-log())をみると直線からずれているのもあります。このような場合に比例ハザード解析は有効なのでしょうか。http://www.h5.dion.ne.jp/~ge3j-ari/stat/surv.html
の「3）Cox回帰分析では」によると指数分布が前提とされているそうです。どの程度指数からはずれているとよくないのでしょうか。

　　　　　[このページのトップへ]

1225.　Re: 比例ハザード解析/R/coxph　detour　　2003/12/04 (木) 18:40

http://www.h5.dion.ne.jp/~ge3j-ari/stat/surv.html
> の「3）Cox回帰分析では」によると指数分布が前提とされているそうです。どの程度指数からはずれているとよくないのでしょうか。

と書きましたが，「指数分布が前提」というのは私の誤読かもしれません。Rのsurvivalのcoxphのhelpをもう少し良く読んでみると，比例ハザードの検定 cox.zphというのがみつかりました。
http://cran.r-project.org/doc/contrib/Fox-Companion/appendix-cox-regression.pdf
の解説を部分的に読んでみると，cox.zphで表示されるpが1に近ければ，その要因は比例ハザード性，0に近ければ非比例ハザード性，ということらしいです。

　　　　　[このページのトップへ]