修复了 Fine&Gray 竞争风险调整模型 (FGR) 中的 coefficient/Offset
Fixed coefficient/Offset in Fine&Gray competing-risk adjusted model (FGR)
我想拟合 Fine&Gray 竞争风险调整模型,包括抵消。在其他类型的模型中,我习惯于可以简单地输入 >offset(x),这将添加一个系数为 1 的偏移量。
我尝试使用包 riskRegression 中的 FGR 函数来做同样的事情。我没有收到警告消息,但我随后注意到模型中有和没有 offset(x) 的系数与其他变量完全相同
示例:
#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)).
有人知道在 FGR 中包含 x 作为偏移量(即系数固定为 1)的方法吗? (编辑:或者另一种计算固定 x 的 a、y 和 z 的正确系数的方法?)
您可以将生存包用于 Fine-Gray 具有偏移量的模型。只需用 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)
我想拟合 Fine&Gray 竞争风险调整模型,包括抵消。在其他类型的模型中,我习惯于可以简单地输入 >offset(x),这将添加一个系数为 1 的偏移量。 我尝试使用包 riskRegression 中的 FGR 函数来做同样的事情。我没有收到警告消息,但我随后注意到模型中有和没有 offset(x) 的系数与其他变量完全相同
示例:
#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)).
有人知道在 FGR 中包含 x 作为偏移量(即系数固定为 1)的方法吗? (编辑:或者另一种计算固定 x 的 a、y 和 z 的正确系数的方法?)
您可以将生存包用于 Fine-Gray 具有偏移量的模型。只需用 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)