我正在努力研究带有子采样的广义随机区组设计（GRBD）的lme。 该实验由2个条带组成，每个条带代表处理（4年/ 2年 - 在作物轮作下）。 每个条带由3个区块组成，每个区块4个实验区，2个附加处理（混合作物/非混合作物）中的每一个的2个图。在每个图中，我采用3个样本进行分解，有些样本缺失（NA）。 我想比较4种处理（4年混合; 4年 - 非混合; 2年混合; 2年 - 非混合）性能。
我的初始模型如下：
lme（decomposed_fraction~ treatment，random = ~1 | block / plot / repetition，data = Dt，na.action = na.omit）
但是，我不确定这个模型是否正确反映了GRBD的分析并正确考虑了实验设计。 欢迎任何建议和解释！
谢谢， 排架