Cox比例風險模型和時間依賴性考克斯R中

模型我有這個survival data描述死亡率三類服務的（沙龍，餐廳和快遞）在十年的時間研究。Cox比例風險模型和時間依賴性考克斯R中

該數據包含三個變量：服務類型（1 =轎車，2 =餐廳，3 =快車），年數（從1到11的整數，其中11表示大於10年）和審查。

我有兩個問題：

1）我已經安裝了Cox比例風險模型，但什麼是檢查比例風險假設的方式。那就是我們假設每個人的風險比和基準風險與時間無關。

2）如何在R中擬合時間相關的Cox模型？

這裏是我的代碼：

#Cox Proportional Hazards 
    cox <- coxph(Surv(Years, Censor) ~ data$`Service Type`) 
    summary(cox)

來源

2016-03-07 Rann

這既是摘但它不是一個編程問題，它是一個統計問題（它可以更好地適用於stats.stackexchange），它看起來也很寬泛，我推薦閱讀[John Fox的生存分析附錄]（https ：//socserv.socsci.mcmaster.ca/jfox/Books/Companion/appendix/Appendix-Cox-Regression.pdf）。它涵蓋了這些主題以及更多內容，並帶有大量示例R代碼。 – Gregor

基線危險不被假定爲與時間無關，數據不支持時間依賴模型，我將回答有關如何檢查比例性的問題在預測中，因爲它是一個班輪。 –

時間依賴性分析需要時間依賴性的預測變量的開始和結束時間。這些數據不支持這種分析。 –

在pkg生存檢查比例的cox.zph功能。請注意，使用「服務類型」爲列名提供的編程麻煩，人們可以通過允許週期來代替空格作爲與函數read.table默認動作很容易避免：

data <- read.table(url("http://www.stat.ufl.edu/~winner/data/bizmort.dat"), col.names=c("Service Type","Years", "Censor") 
# Also note that the censoring indicatior is reversed so will use 1-Censor 
require(survival) 

cox <- coxph(Surv(Years, 1-Censor) ~ factor(Service.Type), data=data) 
# The test fro proportionality: 

> cox.zph(cox) 
         rho chisq  p 
factor(Service.Type)2 0.0306 0.98 0.322 
factor(Service.Type)3 0.0429 1.91 0.167 
GLOBAL     NA 2.33 0.312

來源

2016-03-08 00:31:21

Cox比例風險模型和時間依賴性考克斯R中

回答

相關問題