Python scipy的rv_discrete，多变量和期望值？

如何解决Python scipy的rv_discrete，多变量和期望值？

我刚刚浏览了scipy的rv_discrete函数的文档。我特别想根据我提供的值创建离散支持多元rv。

从documentation看来，只有一个值列表才有可能。

所以我想出了代码：

a = [[1/6,1/6,1/6],[1/6,1/6]]
#print(np.shape(a))
xks=[[1,2,3],[1,3]];
joint = stats.rv_discrete(name='Joint',values=(xks,a));

这正确吗？

此外，既然rv是多维的，我如何使用“期望”？例如，我需要仅沿维度计算期望值。或者只是边缘化。