r语言 - 精细和灰色竞争风险调整模型 (FGR) 中的固定系数/偏移量

Question

我想适合一个罚款＆amp;灰色竞争风险调整后的模型，包括偏移。在其他类型的型号中，我习惯于简单地放入>偏移(x(，这将增加系数1的偏移。我尝试使用软件包风险重试的FGR函数进行相同的操作。我没有收到警告消息，但是我随后注意到，对于其他变量，带有和没有偏移的模型的系数完全相同

示例：

#install.packages(riskRegression")
library(riskRegression)
matrix <- matrix(c(3,6,3,2,5,4,7,2,8,2,
               0.8,0.6,0.4,0.25,0.16,0.67,0.48,0.7,0.8,0.78,
               60,55,61,62,70,49,59,63,62,64,
               15,16,18,12,16,13,19,12,15,14,
               0,2,1,0,1,1,0,1,2,0,
               345,118,225,90,250,894,128,81,530,268),
             nrow=10,ncol=6)
df <- data.frame(matrix)
colnames(df) <- c("x","y","z", "a","event","time")
fit <- FGR(Hist(time,event)~ offset(x)+a+y+z, data=df, cause=1)
fit
fit2 <- FGR(Hist(time,event)~ a+y+z, data=df, cause=1)
fit2

如果您运行此脚本，则可以看到A，Y和Z的系数不会更改，而您没有得到无法使用偏移的警告(因此，显然它只是简单地忽略了偏移(x((。

有人知道将x作为偏移(即以1个系数固定为1(包含在FGR中的方法吗？(编辑：或用固定X计算A，Y和Z正确系数的另一种方法？(

Answer 1

您可以将生存软件包用于带有偏移的细灰模型。只需包装您想要使用offset(var)偏移的变量即可。我将下面的模型设置为模型事件1。请参见下面的代码：

library(survival)
matrix <- matrix(c(3,6,3,2,5,4,7,2,8,2,
                   0.8,0.6,0.4,0.25,0.16,0.67,0.48,0.7,0.8,0.78,
                   60,55,61,62,70,49,59,63,62,64,
                   15,16,18,12,16,13,19,12,15,14,
                   0,2,1,0,1,1,0,1,2,0,
                   345,118,225,90,250,894,128,81,530,268),
                 nrow=10,ncol=6)
df <- data.frame(matrix)
colnames(df) <- c("x","y","z", "a","event","time")
coxph(Surv(time,event==1)~ offset(x)+a+y+z, data=df)