如何调整coxph（R）中的变量？

Question

我很困惑如何调整coxph函数中的变量。我知道可以使用strata（）进行分层，但是如何调整变量？

线性模型中，可以通过以下

（这里有一个我从这个link看到例子）

fit.diamOnMachine <- lm(diameter˜machine) 
diam.adjusted <- residuals(fit.diamOnMachine) 
fit.diamadjmach <- lm(strength ˜ diam.adjusted + machine) 

在coxph，确实一个做同样的事情调整变量？例如，我有年龄，性别，BMI，事件发生时间，结果和1000个人的蛋白质表达值列表。

我想知道如何调整年龄，性别和体重指数在以下coxph型号：

coxph(Surv(Time_till_event, outcome) ~ protein_expression_values, data = data) 

我应该补充一点，我需要调整到这样的公式中的变量：

coxph(Surv(Time_till_event, outcome) ~ protein_expression_values + age + gender + BMI, data = data) 

我在互联网上经历了很多材料，但我似乎无法找到合适的答案。我从来没有参加线性回归课，只是开始在线阅读/学习。如果您有推荐的材料，我也非常感谢。谢谢。

Answer 1

我假设你的核心目标是生存时间模型作为蛋白质表达值的功能，但你也想控制这可能混淆这种关系（即，年龄，性别，BMI）变量。如果是这样，只需在上面建议的公式中包含所有这些术语（即protein_expression_values + age + gender + BMI）。从该模型输出的蛋白质表达系数将表示对于给定年龄，性别和BMI，表达对存活的影响。

作为一个声明，我从来没有见过任何人做，你在你的问题联系起来的调整，所以它可能是我误解你的目标。