我是 R 的新手，我不熟悉 R 中线性混合模型的功能。
我的数据是关于植物生长的，我将害虫作为固定因素（2 个水平），将土壤侵蚀作为一个固定因素（2 个水平），我还想包括它们的相互作用。
试验在三个可能对害虫和土壤侵蚀速率有影响的不同湿度的地点（随机因素）进行。根据植物在叶片形状上的微小差异（嵌套在场地内），植物位于每个场地。
我目前已经用 lme、aov 和 lmer 构建了一些模型，但我不确定哪个是正确的，或者对模型有更好的建议吗？
以下是我提出的模型：
Model1 <- lme(fixed = growth ~ pests*soil, random = ~1|leafshapes, data = plants)`
我不知道如何在模型中包含“站点”。
测试固定因素（害虫和土壤侵蚀）和随机因素（地点）之间的相互作用：
Model2 <- aov(growth ~ pests*soil + Error(sites/leafshapes), data = plants)
