2013-03-01 65 views
2

我比較了兩組觀測的長列表,並且我有每個觀測的t統計量和自由度,但實際上我沒有原始數據。儘管如此,我想獲得與t統計量相關的p值。我知道這是可能的,但我不知道如何讓R報告一個P值。有什麼建議?從t統計得到一個p值

回答

6

task view on Distributions應該是你的第一個停靠港來回答這樣

引述

對於大多數經典的分佈問題,基礎R提供概率 分佈函數(P),密度函數(d),分位數函數 (q)和隨機數生成(r)。除了這個基本功能外,許多CRAN軟件包還提供了其他有用的 分佈。特別是,多元分佈以及 copula可在提供的軟件包中獲得。

您是第t分配後,從概率分佈函數的值,所以你的功能後pt

pt(q, df) 

其中q是你quantile,在這種情況下,你觀察到的測試統計,並且df自由度

我不確定你的意思是the degrees of freedom for *each* observation,我假設你的意思是你知道每個gro的觀察次數向上。

所以,如果你假設組差異中都是平等的,那麼自由度會n1+n2 - 2,否則你將需要佔不等方差和獲得自由的近似度(Welch-Satterthwaite modification

+0

是的!這正是我需要的!謝謝!!! – LauraS 2013-03-01 06:30:02