考虑到食草动物的存在/不存在，我正在运行贝叶斯模型来估计植物上的果实数量。我得到每个均值的后验分布。然后，我运行一个单独的模型来估计每个水果的种子数量，考虑到食草动物的存在/不存在，再次使用每个平均值的后验分布。有没有办法结合这些后验分布来估计总种子产量？这可以像乘以后验一样简单吗？类似于计算每株植物的平均果实数和每颗果实的平均种子数，然后将它们相乘以粗略估计种子产量？

或者，是否可以在一个采样器中同时运行两个模型，然后使用生成的数量。例如，使用 JAGS（或 STAN）：

for(n in 1:N){
    fruitNumber[n] ~ dnorm(fruitHat[n], sd_fruit)
    seedMass[n] ~ dnorm(seedHat[n], sd_seed)

    fruitHat[n] <- B0 + B1*Herb[n]
    seedHat[n] <- G0 + G1*Herb[n]
}

fruitHerb <- B0
fruitNoHerb <- B0+B1
seedHerb <- G0
seedNoHerb <- G0 + G1

totalSeedHerb <- fruitHerb * seedHerb
totalSeedNoHerb <- fruitNoHerb * seedNoHerb