有效地计算β分布置信区间

Question

我正在处理一个非常大的数据集，我需要计算使用每个数据点创建的beta分布的置信区间。我开始使用scipy.stats.beta中的ppf函数，但计算时间太长。有没有更有效的方法来计算β分布的置信区间？ python中的解决方案更可取，但是另一种语言的解决方案也可以与python一起使用。