Как получить результаты теста на соответствие критерия хи-квадрат, используя Goodfit?

Question

Я пытаюсь написать часть программного обеспечения на R, которая находит наиболее подходящее семейство распределений для набора данных, выполняя тест хи-квадрат для данных (в отношении указанного семейства) и находя наилучшее значение хи-квадрат.

Однако при использовании функции goodfit, по-видимому, единственный способ получить статистику хи-квадрат - запустить функцию и использовать команду summary (gf). Это приводит только к удобочитаемому выводу, и мне нужно что-то, что я могу нарисовать в виде gf$chisqvalue, чтобы я мог сравнить его с результатами других тестов, которые я выполняю. Есть ли способ получить эту статистику как переменную?

Answer 1

Попробуйте пакет метлы и команду "tidy" следующим образом:

library(vcd)
library(broom)

x <- rnbinom(200, size = 3, prob = 0.2)
res <- goodfit(x, type = "nbinomial", method = "MinChisq")
summary(res)

dt_res = tidy(summary(res))

Это позволит получить удобочитаемую информацию и сохранить ее в файле данных. Вы можете предпочесть изменить имена столбцов после (или, возможно, нет). Пакет метлы отлично подходит, если вы хотите создать data.frame из статистического теста или вывода модели.