R中的临界t值

Question

我需要确定75%和99%的单侧尾部的临界t值,对于40个自由度.

以下是双侧99%临界t值的代码:

qt(0.01, 40)

但我怎样才能确定单侧临界t值？

Answer 1

您发布的代码为单侧测试提供了临界值(因此,您的问题的答案很简单:

abs(qt(0.25, 40)) # 75% confidence, 1 sided (same as qt(0.75, 40))
abs(qt(0.01, 40)) # 99% confidence, 1 sided (same as qt(0.99, 40))

注意,t分布是对称的.对于双面测试(比如99%置信度),您可以使用临界值

abs(qt(0.01/2, 40)) # 99% confidence, 2 sided

Answer 2

乔希的评论很到位。如果您不是非常熟悉临界值，我建议您使用 qt，阅读手册 ( ?qt) 并查看查找表 ( LINK )。当我第一次从 SPSS 转移到 RI 时，我创建了一个函数，该函数使关键 t 值的查找变得非常容易（我现在永远不会使用它，因为它需要太多时间，并且输出中通常提供的 p 值是一个有争议的问题） ）。这是代码：

critical.t <- function(){
    cat("\n","\bEnter Alpha Level","\n")
    alpha<-scan(n=1,what = double(0),quiet=T)
    cat("\n","\b1 Tailed or 2 Tailed:\nEnter either 1 or 2","\n")
    tt <- scan(n=1,what = double(0),quiet=T)
    cat("\n","\bEnter Number of Observations","\n")
    n <- scan(n=1,what = double(0),quiet=T)
    cat("\n\nCritical Value =",qt(1-(alpha/tt), n-2), "\n")
}

critical.t()