如何使用goodfit檢索卡方擬合優度檢驗的結果？

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。 您可能更喜歡在（或可能不是）之後更改列名稱。 如果您想從統計測試或模型輸出創建數據幀，則掃帚套件非常棒。