返回计算 R 中的 t 统计量

Question

我从独立样本 t 检验中获得了这些信息。

p-value=0.989

sample size for group 1 (n1)= 17

sample size for group 2 (n2)= 18

mean difference= 0.30

我需要这个测试的 t 统计量。我对 p 值的代码进行了反复试验以获得 t 统计量，即 0.014 2*pt(0.014, 33, lower.tail=F) #p=0.989

我需要一个获取 p 值和 df 作为输入并提供 0.014 作为输出的 R 代码。

Answer 1

假设指示函数位于 0 和 6 之间，求其根：

给予：

uniroot(function(x, p) p - 2 * pt(x, 33, lower.tail = FALSE), c(0, 6), p = 0.989)

给予：

$root
[1] 0.01390167

$f.root
[1] 7.863003e-06

$iter
[1] 3

$init.it
[1] NA

$estim.prec
[1] 6.103516e-05

这是另一个例子：

2 * pt(1.901367, 19, lower.tail = FALSE)
## [1] 0.07252959

uniroot(function(x, p) p - 2 * pt(x, 19, lower.tail = FALSE), c(0, 6), p = 0.07252959)

给予：

$根 [1] 1.901365

$f.root [1] -2.672236e-07

$iter [1] 9

$init.it [1] 不适用

$estim.prec [1] 6.103516e-05

Answer 2

不需要uniroot。使用 qt:

> qt(0.989/2, df = 33, lower.tail = FALSE)
[1] 0.01389174

Back calculating t-statistic in R