我有一个问题，我想如何从线性混合效应模型中获取残差（例如，使用R中nlme的lme）以便在lme模型删除了无用的特定变量后绘制数据。在下面的代码中，我已经显示了我正在使用的模型。我关心的变量是grp，而我试图找出的变量（因为它们是令人讨厌的变量）在固定效果和网站随机效果中都是性别和年龄之类的。我想本质上相对于grp变量绘制y_var，但是在模型考虑了性别，年龄和地点之后。谁能建议我如何在R中实现呢？
lme_model = lme(fixed = y_var ~ grp + sex + age, random = ~1|site, data = my_df)