Объект распределения вероятностей ядра
Объект KernelDistribution состоит из параметров, образцового описания и выборочных данных для непараметрического сглаживающего ядро распределения.
Распределение ядра является непараметрической оценкой функции плотности вероятности (PDF) случайной переменной.
Распределение ядра использует следующие опции.
| Опция | Описание | Возможные значения |
|---|---|---|
Kernel | Тип функции ядра | normal, box, triangle, epanechnikov |
BandWidth | Параметр сглаживания ядра | BandWidth > 0 |
Существует несколько способов создать объект распределения вероятностей KernelDistribution.
Соответствуйте распределению к данным с помощью fitdist.
В интерактивном режиме соответствуйте распределению к данным с помощью приложения Distribution Fitter.
cdf | Кумулятивная функция распределения |
icdf | Обратная кумулятивная функция распределения |
iqr | Межквартильный размах |
mean | Среднее значение распределения вероятностей |
median | Медиана распределения вероятностей |
negloglik | Отрицательный loglikelihood распределения вероятностей |
pdf | Функция плотности вероятности |
random | Случайные числа |
std | Стандартное отклонение распределения вероятностей |
truncate | Усеченный объект распределения вероятностей |
var | Отклонение распределения вероятностей |
Загрузите выборочные данные. Визуализируйте терпеливые данные о весе с помощью гистограммы.
load hospital
histogram(hospital.Weight)
Гистограмма показывает, что данные имеют два режима, один для пациенток и один для штекерных пациентов.
Создайте объект распределения вероятностей путем подбора кривой распределению ядра к терпеливым данным о весе.
pd_kernel = fitdist(hospital.Weight,'Kernel')pd_kernel =
KernelDistribution
Kernel = normal
Bandwidth = 14.3792
Support = unbounded
Для сравнения создайте другой объект распределения вероятностей путем подбора кривой нормальному распределению к терпеливым данным о весе.
pd_normal = fitdist(hospital.Weight,'Normal')pd_normal =
NormalDistribution
Normal distribution
mu = 154 [148.728, 159.272]
sigma = 26.5714 [23.3299, 30.8674]
Задайте x значения и вычислите PDF каждого распределения.
x = 50:1:250; pdf_kernel = pdf(pd_kernel,x); pdf_normal = pdf(pd_normal,x);
Постройте PDF каждого распределения.
plot(x,pdf_kernel,'Color','b','LineWidth',2); hold on; plot(x,pdf_normal,'Color','r','LineStyle',':','LineWidth',2); legend('Kernel Distribution','Normal Distribution','Location','SouthEast'); hold off;
При подборе кривой распределению ядра вместо одномодового распределения такой, когда нормальное показывает отдельные режимы для розеточных и штекерных пациентов.