概率编程是否足以解决类似于“二项比例置信区间”的问题？

Question

我有以下设置（摘要）：

• 参数p的伯努利变量z（未知）
• 其值通常取决于变量x（具有很多维）
• 我有很多观测值（x，z）

我想估计z的参数p。

我已经考虑过“二项比例置信区间”之类的方法，但是这些技术似乎并未利用任何有关x的信息：x具有非常复杂的结构（高维），各个维之间存在很多依存关系，我可以用概率编程很好地表达这一点。但是我不确定这是否是正确的工具，您对此有何看法？