Question

对于40度自由度，我需要确定75％和99％的单侧尾部的临界t值。

以下是双侧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

Josh的评论很明显。如果您不熟悉我建议使用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()

R中的临界t值

2 个答案: