Функция плотности полиномиальной вероятности
Y = mnpdf(X,PROB)
Y = mnpdf(X,PROB) возвращает pdf для полиномиального распределения с вероятностями PROB, оценивается в каждой строке X. X и PROB m-на-k матриц или 1-на-k векторов, где k - число полиномиальных ячеек или категорий. Каждая строка PROB должен суммироваться с единицей и размерами выборки для каждого наблюдения (строки X) задаются строковыми суммами sum(X,2). Y является вектором m-by-1, и mnpdf вычисляет каждую строку Y использует соответствующие строки входных данных или при необходимости реплицирует их.
Вычисление pdf полиномиального распределения с размером выборки n = 10. Вероятности: p = 1/2 для результата 1, p = 1/3 для результата 2, и p = 1/6 для результата 3.
p = [1/2 1/3 1/6]; n = 10; x1 = 0:n; x2 = 0:n; [X1,X2] = meshgrid(x1,x2); X3 = n-(X1+X2);
Вычислите pdf дистрибутива.
Y = mnpdf([X1(:),X2(:),X3(:)],repmat(p,(n+1)^2,1));
Постройте график pdf на трехмерной фигуре.
Y = reshape(Y,n+1,n+1); bar3(Y) h = gca; h.XTickLabel = [0:n]; h.YTickLabel = [0:n]; xlabel('x_1') ylabel('x_2') zlabel('Probability Mass') title('Trinomial Distribution')
Обратите внимание, что визуализация не отображается x3, которая определяется ограничением x1 + x2 + x3 = n.