我正在使用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)
的置信区间,这是具有三个风险因素的人的危险比。
我想知道是否有任何形式可以使用包survival
或其他包