我使用 coxph 作为生存模型，使用三个定性变量作为独立变量。我知道它们之间没有交互，所以我只是做了以下模型：
fit <- coxph(Surv(Time, status) ~ factor(v1) + factor(v2) + factor(v3)
我的问题是我需要它们之间产品的置信区间。我可以使用公式手动估计两个乘积中的一个
Var(v1*v2) = v1^2 Var(v2) + v2^2 Var(v1) + 2v1v2COV(v1v2)
但是对于他们三个的乘积我不知道怎么做，我的意思是，我分别有 v1、v2 和 v3 的三个估计器 b1、b2、b3。
我想知道如何估计exp(v1 + v2 + v3)具有三个风险因素的人的风险比的置信区间。
我想知道是否有任何形式可以用 R、包survival或另一个包来做