我使用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 或与另一个包