Почему при запуске ggttest возникает ошибка?

Question

Когда я запускаю t-критерий для числа c и дихотомической переменной, проблем нет, и я вижу результаты. Проблема в том, когда я запускаю ggttest того же t-теста. Возникает ошибка и говорит, что одна из моих переменных не найдена. Я не знаю, почему это происходит. Набор данных 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")

Answer 1

Проблема возникает с этими строками кода в 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")

Answer 2

t.test не сохраняет ваши данные в выходных данных, поэтому вы не сможете извлечь данные из списка выходных данных t.test.

Единственный способ использовать формулу:

library(gginference)
t_test <- t.test(questionnaire$pulse ~ questionnaire$gender)
ggttest(t_test)

Оригинальный ответ здесь: Как извлечь набор данных из объекта "htest" при использовании формулы в r