Question

我对使用python计算学生的置信区间感兴趣。

我在Mathematica中使用StudentTCI（）函数，现在需要在python http://reference.wolfram.com/mathematica/HypothesisTesting/ref/StudentTCI.html中编写相同的函数

我不太确定如何自己构建这个函数，但在我开始之前，这个函数是否在python的某个地方？喜欢numpy？（我没有使用numpy，我的顾问建议如果可能的话不使用numpy。）

解决此问题的最简单方法是什么？我可以将源代码从numpy中的StudentTCI（）（如果存在）复制到我的代码中作为函数定义吗？

编辑：我将需要使用python代码构建Student TCI（如果可能的话）。安装scipy已经变成了死胡同。我遇到了其他人都遇到的同样的问题，如果设置这么长时间，我就无法为我分发的代码要求Scipy。

任何人都知道如何在scipy版本中查看算法的源代码？我想我会把它重构为python定义。

Answer 1

我猜您可以使用scipy.stats.t并使用interval方法：

In [1]: from scipy.stats import t
In [2]: t.interval(0.95, 10, loc=1, scale=2)  # 95% confidence interval
Out[2]: (-3.4562777039298762, 5.4562777039298762)
In [3]: t.interval(0.99, 10, loc=1, scale=2)  # 99% confidence interval
Out[3]: (-5.338545334351676, 7.338545334351676)

当然，如果你愿意，你可以自己动手。让它看起来像Mathematica：

from scipy.stats import t


def StudentTCI(loc, scale, df, alpha=0.95):
    return t.interval(alpha, df, loc, scale)

print StudentTCI(1, 2, 10)
print StudentTCI(1, 2, 10, 0.99)

结果：

(-3.4562777039298762, 5.4562777039298762)
(-5.338545334351676, 7.338545334351676)

学生在python中的置信区间

1 个答案: