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

Question

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

但是，在使用 goodfit 函数时，似乎检索卡方统计量的唯一方法是运行该函数并使用 summary(gf) 命令。这只会产生人类可读的输出，我需要一些可以以 gf$chisqvalue 的形式绘制的东西，以便我可以将它与我正在运行的其他测试的结果进行比较。有没有办法将此统计信息作为变量检索？

Answer 1

尝试使用 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 包非常棒。