我正在尝试作为练习,因为我仍然是相当新的python和编程,以制作一个脚本,它采用一个样本数据池,并使用从表中的值使用,以产生比stdev更准确的偏差。
示例: 10个样本,我希望列中的t表值为0.975。 在这种情况下,它的10-1,第9行在0.975列中具有t值2.262
有没有简单的方法只需输入python或任何额外的库,我只想从行X列0.975中获取T值并在进一步计算中使用它?
或者我是否必须找到CVS文件并尝试以这种方式导入值? 我试图通过scipy.stats函数,但老实说,对于刚刚开始并且之前没有使用p值或类似统计信息的人来说,它有点压倒性。
答案 0 :(得分:5)
您可以使用scipy.stats.t
的ppf
方法(即分位数函数):
In [129]: from scipy.stats import t
In [130]: alpha = 0.025
In [131]: t.ppf(1 - alpha, df=9)
Out[131]: 2.2621571627409915
t.ppf()
最终会调用scipy.special.stdtrit
,因此您也可以使用该功能,避免经历t.ppf()
的轻微开销:
In [141]: from scipy.special import stdtrit
In [142]: alpha = 0.025
In [143]: stdtrit(9, 1 - alpha)
Out[143]: 2.2621571627409915
(如果您想知道该功能如何以名称stdtrit
结束,则 S tudent T D 是 TR ibution函数我与 T 相反。清楚,对吗?)
答案 1 :(得分:0)
我在Python中找不到一个,但是我发现rpy2在Python中提供了许多R功能,还有它们的t.test功能。