如何使用 goodfit 检索卡方拟合优度检验的结果?

How do I retrieve the results of a chi-squared goodness-of-fit test using goodfit?

我正在尝试用 R 编写一个软件,通过对数据(关于所述系列)执行卡方检验并找到最佳分布系列来找到最适合一组数据的分布系列卡方值。

然而,当使用 goodfit 函数时,似乎检索卡方统计量的唯一方法是通过 运行 该函数并使用 summary(gf) 命令。这只会产生人类可读的输出,我需要一些可以 gf$chisqvalue 形式绘制的东西,这样我就可以将它与我 运行 的其他测试的结果进行比较。有什么方法可以将此统计信息作为变量检索吗?

像这样尝试 broom 包和命令 "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 中。 您可能更愿意在(或可能不)之后更改列名。 如果您想从统计测试或模型输出创建 data.frame,那么 broom 包非常有用。