Я пытаюсь вычислить некоторую простую статистику соответствия качества для различных дискретных дистрибутивов, которые я снабдил пакетом fitdistrplus. Однако, когда я пытаюсь запустить функцию gofstat для всех трех функций одновременно, как это описано в документации к пакету, R выдает ошибку, которую я не могу исправить.
Я посмотрел на ошибку, но не могу найти ее специально для такого рода проблем. Функция работает, когда я пишу функцию 3 раза для каждого дистрибутива, так что это связано с моим синтаксисом при написании кода для 3 дистрибутивов одновременно.
#Fit a Poisson distribution
fp <- fitdist(raw_data$summary, "pois")
summary(fp)
#Fit a Negative Binomial distribution
fb <- fitdist(raw_data$summary, "nbinom")
summary(fb)
#Fit a Geometric Distribution
fg <- fitdist(raw_data$summary, "geom")
summary(fg)
#Now we look at the Chi-sq statistic to ascertain which is the best fit
gofstat(list(fp,fb,fg), fitnames=c("Poisson", "Negative Binomial", "Geometric"))
Я получаю сообщение об ошибке при запуске этого кода:
Error in names(Chi2temp$chisqpvalue) <- names(Chi2temp$chisqdf) <- fitnames[i] : attempt to set an attribute on NULL