我有以下问题： 我正在分析不同基因型（A，B和C）之间种子囊的数量 每个基因型我有4个重复，在每个重复中，我有8个植物。以下是数据示例： 基因型Replicate_ID Plant_ID Seed_capsules A 1 1 6 A 1 2 10 A 1 3 15 B 2 1 100 B 2 2 40 B 2 3 63 C 3 1 80 C 3 2 90 C 3 3 100
我在数据上使用了glmer，但我不确定我的随机是replicate_ID还是plant_ID或两者。这是我到目前为止尝试过的一个例子： freplicate_ID＆lt; - factor（newdata $ replicate_ID） fplant_ID＆lt; - factor（newdata $ plant_ID） m5＆lt; - glmer（seed_capsules~pop_ID +（1 | freplicate_ID / fplant_ID），family = poisson，data = data）
此外，我如何获得这种模型的诊断图？如何比较基因型？我需要在模型上运行lsmeans吗？像这样例如： lsmeans（m5，pairwise~pop_ID，data = newdata）。 我提前感谢你的帮助。 最好， 安娜