Задача по ТерВеру в Matlab

[Serbin]

Задача(приблизительный перевод):

Цитата:

В одном городе записывали на протяжении 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?

И как можно реализовать последний пункт?

[Александрович]

В последнем пункте использовать критерий хи-квадрат.