pymc3（适用于ARIMAX模型）

Question

我有一个观测序列作为三个潜在随机序列之和

其中F和G是解释性变量。观察到F，G和O。 pymc3是否可以在贝叶斯框架（假设a，b，beta具有正常先验）下对此建模？

我知道Z可以定义为pymc3.distributions.timeseries.GaussianRandomWalk，但是如何定义X和Y？