当我为数字和二分变量运行t检验时,没有问题,我可以看到结果。问题是当我运行相同的t检验的ggttest时。发生错误,并说找不到我的变量之一。我不知道为什么会这样。我使用的aml数据集来自程序包引导。在下面您可以看到代码:

https://i.stack.imgur.com/7kuaA.png

library(gginference)
time_group.test16537 = t.test(formula = time~group,
                              data = aml,
                              alternative = "two.sided",
                              paired = FALSE,
                              var.equal = FALSE,
                              conf.level = 0.95)
    time_group.test16537


    ggttest(time_group.test16537,
            colaccept="lightsteelblue1",
            colreject="gray84",
            colstat="navyblue")

最佳答案

问题来自ggttest中的以下代码行:

datnames <- strsplit(t$data.name, splitter)
len1 <- length(eval(parse(text = datnames[[1]][1])))
len2 <- length(eval(parse(text = datnames[[1]][2])))

它试图找到组和时间的长度,但没有看到它来自data.frame。相当糟糕的错误...

对于您的情况,假设每个组中少于30个,并且它绘制t分布,所以:
library(gginference)
library(boot)
gginference:::normt(t.test(time~group,data=aml),
colaccept = "lightsteelblue1",colreject = "grey84",
colstat = "navyblue")

r - 为什么当我运行ggttest时出现错误?-LMLPHP

关于r - 为什么当我运行ggttest时出现错误?,我们在Stack Overflow上找到一个类似的问题:https://stackoverflow.com/questions/61034288/

10-12 19:44