gpt4 book ai didi

r - 从pairwise.t.test 中提取t 值

转载 作者:行者123 更新时间:2023-12-04 10:33:55 26 4
gpt4 key购买 nike

我想知道是否可以从 pairwise.t.test 函数中提取 t 值,因为它只报告 p 值。
在运行重复测量方差分析后,我使用 pairwise.t.test() 进行多重比较,该方差报告报告了显着的主效应。

非常感谢您提前!

最佳答案

在尝试解决类似问题时,可以尝试一些简单的事情(当然,不能保证它们在任何给定的情况下都能工作)。首先要尝试的是查看文档( ?pairwise.t.test )并查看 Value 下列出的内容.在这种情况下,它只说:

Object of class "pairwise.htest"



尝试的第二件事是获取一个示例对象并将其分配给一个变量,然后您可以运行 str(obj) .以下使用文档中的示例:
attach(airquality)
Month <- factor(Month, labels = month.abb[5:9])
obj <- pairwise.t.test(Ozone, Month)
str(obj)
List of 4
$ method : chr "t tests with pooled SD"
$ data.name : chr "Ozone and Month"
$ p.value : num [1:4, 1:4] 1 0.000264 0.000195 1 NA ...
..- attr(*, "dimnames")=List of 2
.. ..$ : chr [1:4] "Jun" "Jul" "Aug" "Sep"
.. ..$ : chr [1:4] "May" "Jun" "Jul" "Aug"
$ p.adjust.method: chr "holm"
- attr(*, "class")= chr "pairwise.htest"

不幸的是,这并没有显示我们想要看到的(例如,类似 $ t.stat 的东西)。

您的最后一个选择是查看代码。您可以通过在命令提示符下键入不带括号的函数调用来获取它:
> pairwise.t.test
function (x, g, p.adjust.method = p.adjust.methods, pool.sd = !paired,
paired = FALSE, alternative = c("two.sided", "less", "greater"),
...)
{
<code omitted>
if (pool.sd) {
<code omitted>
}
}
else {
<code omitted>
compare.levels <- function(i, j) {
xi <- x[as.integer(g) == i]
xj <- x[as.integer(g) == j]
t.test(xi, xj, paired = paired, alternative = alternative,
...)$p.value
}
}
PVAL <- pairwise.table(compare.levels, levels(g), p.adjust.method)
ans <- list(method = METHOD, data.name = DNAME, p.value = PVAL,
p.adjust.method = p.adjust.method)
class(ans) <- "pairwise.htest"
ans
}

关键部分是定义的函数 compare.levels ,这仅保留来自基础 t 检验的 p 值。因此,您问题的直接答案是否定的,您无法提取 t 统计量。

关于r - 从pairwise.t.test 中提取t 值,我们在Stack Overflow上找到一个类似的问题: https://stackoverflow.com/questions/28119426/

26 4 0
Copyright 2021 - 2024 cfsdn All Rights Reserved 蜀ICP备2022000587号
广告合作:1813099741@qq.com 6ren.com