Сгенерируйте набор выборочных данных и подгоняйте кусочное распределение с хвостом Парето к данным при помощи paretotails. Найдите количество сегментов в подобранном распределении с помощью функции объекта nsegments.

Сгенерируйте набор выборочных данных, содержащий 10% выбросов в правом хвосте.

rng('default');  % For reproducibility
right_tail = exprnd(5,100,1);
center = randn(900,1);
x = [center;right_tail];

Создайте paretotails объект путем подгонки кусочно-линейного распределения в x. Задайте контуры хвостов с помощью нижней и верхней кумулятивных вероятностей хвоста. Проходите в 0 и 0,9 так, чтобы подгоняемый объект не содержал нижний хвостовой сегмент и состоял из эмпирического распределения для нижнего 90% набора данных и обобщенного распределения Парето (GPD) для верхних 10% набора данных.

pd = paretotails(x,0,0.9)

pd = 
Piecewise distribution with 2 segments
   -Inf < x < 1.73931  (0 < p < 0.9): interpolated empirical cdf
    1.73931 < x < Inf  (0.9 < p < 1): upper tail, GPD(0.643752,1.62246)

Верните количество сегментов в pd при помощи nsegments функция.

n = nsegments(pd)

n = 2

Вы также можете получить количество сегментов при помощи NumSegments свойство. Доступ к NumSegments свойство при помощи записи через точку.

pd.NumSegments

ans = 2