是否有R函數可用於多處理Cox回歸？

Question

例如，如果我想知道吸煙對肺癌的影響，數據如下

##   ID  age sex basemort  time2 assessment smokingstatus mortality
## 1  1 45.1   1    6.659     14          1          0         0
## 2  1 45.1   1    6.659     16          2          1         0
## 3  1 45.1   1    6.659     23          3          1         1
## 4  2 44.1   1    0.943      8          1          0         0
## 5  2 44.1   1    0.943     20          2          2         1
## 6  2 44.1   1    0.943     22          3          3         1

列吸煙狀態代表處理（暴露），這是一個序數變量，其中1 =禁止吸煙，2 =偶爾吸煙，3 =吸煙年數> 10＆<20，4 =吸煙年數> = 20。 是否有R函數可以像這樣進行多元治療進行cox回歸？ 感謝您的任何建議！

Answer 1

您可以看一下軟件包survival和survminer ，有一個函數可以進行cox回歸。

例如，如果您要根據吸煙狀況來衡量隨時間變化的死亡率，則可以執行以下操作：

library(survival)
coxph(Surv(time2, mortality) ~smokingstatus, data =NameofYourDataset)