У меня есть набор данных с прибытием автомобилей в минуту.
Я нарисовал гистограмму и подогнал ее к распределению Пуассона со следующими R-кодами.
#Aladdin Arrivals
Datast <- read.csv("Vehiclecount.csv", header = T, sep=";", dec=",")
hist(Datast$Arrival, xlab="Arrivals",
probability = TRUE,col=16, ylim = c(0,0.2), xlim =c(0, 30),
main = "Arrivals from Aladdin Street")
lines(dpois(x=0:25, lambda=13.20), col=2,lwd=3)
legend("topright", c("Probability of Vehicle Arrivals ",
"Poisson Distribution Curve"), fill=c(col=16, col=2))
Приведенный выше код успешно запустился, и я получил подогнанные линии поверх гистограммы.
Но когда я хочу использовать функцию goodfit()
, чтобы узнать, каково значение p, я получил следующую ошибку;
«Ошибка в оптимизации (чи2, диапазон (количество)): «xmin» не меньше, чем «xmax»»
dfs <- dpois(x=1:25, lambda=13.20)
summary(dfs)
goodfit(dfs, type="poisson", method="MinChisq")
Как я могу решить эту проблему? Есть ли другая функция для использования?