R 用调整后的协变量绘制线性模型

我试图绘制物种丰度（X轴）与基因表达（Y轴）的线性散点图。然而，有了这个，我就有了我的基因表达的BMI和年龄数据，我需要在模型中加以说明，这样基因表达数据的残差就可以根据物种丰度来绘制

我有一个大的数据集，里面有我的物种丰富度和基因表达值，还有两个额外的列，包含每个个体的BMI和年龄

我的一小部分数据如下所示：

MyData <- ENSG00000056277.11 = c(1.20347265388767, 
    -0.21867257952087, -1.12334136458201, 2.03112643012426, 0.962974348122562), ENSG00000102078.11 = c(-0.926983693379893, 
    0.736304974879333, 1.08736731090216, 0.255232380993567, -1.0414482826186), ENSG00000134597.9 = c(0.562650927598327, 
    -0.469918965753258, -0.368421398397285, -0.752772477131681), k__Archaea = c(NA, -0.022663988, 0.206137726, NA, -0.054418187), Clostridia = c(0.096692281, 0.21466473, NA, NA, 0.245113763), Firmicutes = c(-0.18336613, -0.158690334, NA, -0.148194045, 
NA), BMI = c(23.49590602, 
    21.78587367, 31.84017809, 21.8426439, 26.36642452), Age = c(66.65479452, 
    66.65205479, 69.8, 69.80273973, 68.16164384)

---

MyData共享的数据不工作并返回错误。
LM <- lm(k__Archaea ~ ENSG00000056277.11)