我一直在用 Scipy 做一些数据分析。到目前为止，我通过连续分布实现了这一点：

• 我可以使用最大似然拟合将概率分布拟合到一组数据点。例如使用stats.chi2.fit(data_points).
• 我可以测试我的样本是否来自特定的分布。例如使用stats.kstest(data_points, 'norm')

我尝试对离散分布（如泊松、二项式和几何）进行相同的分析，但 Scipy 似乎没有为离散分布提供类似的功能。拟合和测试离散分布的标准方法是什么？