我编写了一个执行许多单样本t检验的Perl脚本。我得到了数以千计的自由度(df)统计量。我需要升级脚本以返回它们的p值(在表中手动查找它们太多了)。我可以用t-statistic和d.f作为输入来使用某种公式吗?我希望有人可以帮助我,提前多多谢谢!
A.A。
答案 0 :(得分:7)
使用Statistics::Distributions似乎非常简单:
print Statistics::Distributions::tprob($dof,$tstat);
答案 1 :(得分:6)
搜索MetaCPAN会发现以下内容:
答案 2 :(得分:3)
如果您正在进行双尾测试,那么您的p值= 2 * P(T> t),其中t是您计算的测试统计量。基本上,你需要一种方法来模拟T-dist以便整合(T-dist从t到INFINITY)。这是一个演示:http://www.stat.tamu.edu/~west/applets/tdemo.html
我不熟悉Perl及其库,但希望这可以让你开始。您可以编写一个基本的积分器并检查一些值,以确保它足够准确。