Задача(приблизительный перевод):
Цитата:
В одном городе записывали на протяжении 100 дней количество аварий. Получили результаты:
Кол-во аварий Частота
0 8
1 28
2 31
3 18
4 9
5 6
Определить эмпирически, какое здесь распределение.
Найти оценку максимального правдоподобия (maximum-likelihood estimation).
Построить статистический тест для проверки, соответствует ли нами предложенное теоретическое распределение этому эмпирическому распределению.
Решить задачу надо через matlab.
С первым просто, по графику это распределение Пуассона.
Со вторым имею код:
Код:
X = [8; 28; 31; 18; 9; 6];
lambdahat = poissfit(X)
phat = mle(X)
lambdahat =
16.6667
phat =
16.6667 9.8601
Можно ли использовать здесь poissfit(), если мы только предполагаем, что это распределение Пуассона? Почему при mle() выдаётся 2 значения? Какой физический смысл 16.6667?
И как можно реализовать последний пункт?