我正在估计一个障碍模型，它由两个数据生成过程组成：一个用于零/正关系，另一个用于大于零的正数。第一部分通常使用二元 logit 模型，第二部分通常使用零截断泊松或负二项式模型。

在 R 中，psclhurdle()包中的命令共同估计模型，并返回单个 AIC 和对数似然值。我需要加入一些随机效果，因此我使用的是同名包中的命令。这个包没有一个单一的障碍模型命令，所以我必须单独估计它们（首先是二进制，然后是泊松）。glmmADMB()

我的问题：应该如何计算两个单独估计模型的联合对数似然？我认为这是可能的，因为它hurdle似乎就是这样做的。