在我的硕士论文中，我正在研究林奇综合征患者成年期体重增加、BMI、腰围、腰臀比、腰高比和结直肠肿瘤（结直肠癌和结直肠腺瘤）之间的关系。为了对此进行调查，我将执行一个 cox 比例风险回归模型，并对这些暴露进行分类。但是，我还将通过使用带有 cox 的多元分数多项式来执行连续曝光。但是，我的数据集中有来自同一家庭的人。因此，我需要使用稳健的方差估计。我能够通过根据他们的家庭号码对参与者进行聚类来进行分类分析（我将此文档用于我的脚本：https ://www.rdocumentation.org/packages/survival/versions/3.2-10/topics/考克斯）。
因此，这是我的脚本示例，其中我将 BMI 作为 2 个类别的曝光，协变量基于 DAG：
ph_finalbminew1.2 <- coxph(Surv(ftime_CRT, CRT_afterstart) ~ BMIcatnew + age + totalindex + SumOfkCal + ultraprocessed + fruit_vegetables + fish + sex, cluster = famnr, data = df_mbulbaai)
该脚本有效。但是，当我尝试在 mfp 模型中执行集群项时，它不起作用。
mfp 模型的脚本试验：
df_bmi <-mfp(Surv(ftime_CRT, CRT_afterstart)~ fp(BMIcenter, df=4)+fp(SumOfalcohol, df=4)+fp(SumOfkCal, df=4)+fp(fruit_vegetables, df=4)+fp(totalindex, df=4)+ fp(red_pmeat, df=4)+ fp(SSBs, df=4)+fp(ultraprocessed, df=4)+fp(agecenter, df=4)+fp(fish, df=4)+SMOKER+educat_both+sex, cluster = famnr, family = cox, data=df_mbulbaai, verbose=TRUE)
```r
R says: `unused argument (cluster = famnr)`
df_bmi <-mfp(Surv(ftime_CRT, CRT_afterstart)~ fp(BMIcenter, df=4)+fp(SumOfalcohol, df=4) + fp(SumOfkCal, df=4) + fp(fruit_vegetables, df=4) + fp(totalindex, df=4) + fp(red_pmeat, df=4) + fp(SSBs, df=4)+ fp(ultraprocessed, df=4) + fp(agecenter, df=4) + fp(fish, df=4)+ SMOKER + educat_both + sex, cluster(famnr), family = cox, data=df_mbulbaai, verbose=TRUE)
R 说： Error in cluster(famnr) : object 'famnr' not found
谁能帮我解决这个问题？甚至可以在 R 中做到这一点吗？谢谢你的检查。我很感激！