我寫了一個Perl腳本,執行許多單樣本t檢驗。我獲得了數以千計的他們的自由度(df)的統計。我需要升級腳本以返回它們的p值(有太多無法在表中手動查找它們)。有什麼樣的公式可以用t-statistic和d.f作爲輸入嗎?我希望有人能夠幫助我,提前多謝!如果我有t統計量和d.f.,我如何計算p值? (在Perl中)?
A.A.
我寫了一個Perl腳本,執行許多單樣本t檢驗。我獲得了數以千計的他們的自由度(df)的統計。我需要升級腳本以返回它們的p值(有太多無法在表中手動查找它們)。有什麼樣的公式可以用t-statistic和d.f作爲輸入嗎?我希望有人能夠幫助我,提前多謝!如果我有t統計量和d.f.,我如何計算p值? (在Perl中)?
A.A.
使用Statistics::Distributions似乎很簡單:
print Statistics::Distributions::tprob($dof,$tstat);
如果您正在進行雙尾檢驗,那麼您的p值= 2 * P(T> t),其中t是您計算的檢驗統計量。所以基本上,您需要一種方法來模擬T-dist以整合(T-dist從t到INFINITY)。這裏是一個演示:http://www.stat.tamu.edu/~west/applets/tdemo.html
我不熟悉Perl及其庫,但希望這可以讓你開始。您可以編寫一個基本的集成器並檢查一些值以確保它足夠準確。
一個搜索MetaCPAN的揭示了以下內容:
1)。沒有足夠的信息(2)。向我們展示您的編碼嘗試... – 2009-05-05 03:27:01
這是一個非常大的腳本。它工作正常。我需要的只是一個公式或一段使用t統計量的腳本,它是d.f. (我已經計算並存儲到內存中)來計算一個p值.. :-) 但是如果有什麼特定的內容想要了解腳本,我會很樂意爲您提供信息。 – Abdel 2009-05-05 03:39:03
沒有更多的信息,沒有辦法告訴你想要什麼/需要什麼。 – 2009-05-05 04:02:32