fit<- rma.mv(yi, vi, mods = ~ arm + crop.type - 1,
             random = ~ arm | paper/study,
             data = dat,
             struct = "HCS"
)