当我为数字和二分变量运行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时出现错误?,我们在Stack Overflow上找到一个类似的问题:https://stackoverflow.com/questions/61034288/