是否可以在PyMC3中使用NUTS，其模型涉及参数的特征分解？

Question

我有一个模型，其中可能性涉及计算矩阵中所有项的总和

P = U exp（tD）U ^ -1

哪里

UDU ^ -1 = Q

和Q是我的参数矩阵。如果我想在PyMC3中使用NUTS，NUTS必须能够计算P中所有元素相对于Q中每个元素的导数。这是否可能使用PyMC3使用的Theano中的象征性区分？

Answer 1

PyMC3使用Theano进行计算和自动编辑。 Theano对张量代数（其中矩阵代数是子集）有很好的支持，所以我认为应该支持你的模型。