如何为数据集的每个级别执行线性混合模型？

Question

我有一个数据集，变量响应为“ Yield”。在许多“环境”（Env）和“年”中，每个“ CTL”的收益率都是依赖的，每个“ CTL”具有三个“ rep”（重复）。我需要为每个环境执行线性混合模型（我有很多）。为此，我使用了以下LMM：

My data is like this (with many rows more):

我的线性混合模型是跟随模型（对于每个环境）：

fit.f1<- lmer(Yield~ (1|YR)  + (1|Rep) + (1|CLT) + 
                         (1|YR:CLT), data = mydata)

我一直在尝试执行“ for循环”以针对每种环境进行分析：

out <- for(i in 1:length(levels(mydata$env))) {
  fit.f1<- lmer(Yield~ (1|YR)  + (1|RP) + (1|CLT) + 
                         (1|YR:CLT), data = "[i]")

我希望获得一个具有每种环境响应的向量，并具有来自每种环境的混合模型的拟合数据。