在R中固定变量的特定水平上的治疗效果

Question

我有4个变量（变量，比率，DAP和总生物量），其中变量和比率是我的治疗方法，而DAP是固定变量。而总生物量是一个响应变量。 我有10级DAP，4级变种和2级Nrates。

现在，我想看看在每个DAP上，组合治疗（即相互作用）对我的反应变量的影响（因此，对于每个组合治疗，DAP保持固定）。关于如何在R中执行此操作的任何建议？欣赏它！