如何調整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，表達對存活的影響。

作爲一個聲明，我從來沒有見過任何人做，你在你的問題聯繫起來的調整，所以它可能是我誤解你的目標。