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

Question

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

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

Answer 1

Попробуйте пакет метлы и команду «аккуратным», как это:

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. Возможно, вы захотите изменить имена столбцов после (или, возможно, нет). Брум-пакет замечательный, если вы хотите создать data.frame из статистического теста или вывода модели.