我正在尝试使用brms软件包来研究3种情况之间四种蝴蝶物种的寄生虫发生率变化。除一种外，所有其他物种均已在每种条件下进行了测试，因此出现了我的问题。物种4未在条件3下进行测试。
我的完整模型是：
mod15=brm(#ind_parasitized|trials(#ind_parasitized+#ind_nonparasitized)~0+condition*butterfly_species,data=dat3,iter=5000,warmup=2000,core=4,family=binomial)
此模型返回未能收敛的错误。但是，从模型摘要来看，Rhat，Eff。样本，并且链的混合是好的，但是很明显，因为区域3中缺少构造类4，模型无法忽略/排除。
四处阅读，我已经看到了如何处理缺失的值，但是我不知道该如何处理缺失的约束。我以为可以使用brm_multiple来解决这个问题，但是这种方法不起作用，因为在这种情况下，我尝试测试不再存在的constrast。
任何建议将不胜感激！