Список функций Statistics Toolbox

  В оглавление \ К следующему разделу \ К предыдущему разделу

Синтаксис:

[f,x] = ecdf(y)
[f,x,flo,fup] = ecdf(y)
[...] = ecdf(y,'param1',value1,'param2',value2,...)

Описание:

[f,x] = ecdf(y) - расчет значений эмпирической функции распределения на основе оценки Каплана-Мейера, где y – вектор исходных данных, f – вектор значений эмпирической функции распределения рассчитанной для упорядоченного ряда исходных данных х.

[[f,x,flo,fup] = ecdf(y) - кроме значений f и х позволяет определить нижнюю flo и верхнюю fup границы доверительных интервалов для значений эмпирической функции распределения. Расчет значений границ доверительных интервалов проводится по формуле Гринвуда.

[...] = ecdf(y,'param1',value1,'param2',value2,...) - дополнительные параметры 'param1', value1, 'param2', value2, ... задаются в виде строки и соответствующего ей вектора значений. Дополнительные параметры позволяют задать вид цензурированности наблюдений, частоту значений, уровень значимости и тип выходного результата. Возможные значения строковой переменной 'param' и функции вектора value приведены в следующей таблице:

Примеры:

Пример построения теоретической и эмпирической функции распределения с границами доверительных интервалов для выборки из 50 элементов, распределенных по экспоненциальному закону.

Генерация двух выборок y, d на 50 элементов распределенных по экспоненциальному закону с параметрами распределения равными 10 и 20 соответственно.
>> y = exprnd(10,50,1);
>> d = exprnd(20,50,1);

Исследуемая выборка t определяется как вектор минимальных значений при поэлементном сравнении элементов векторов y и d.
>> t = min(y,d);

Определение условий цензурирования значений исследуемой выборки t.
>> censored = (y>d);

Расчет значений эмпирической функции распределения и границ доверительных интервалов.
>> [f,x,flo,fup] = ecdf(t,'censoring',censored);

Ступенчатые графики эмпирической функции распределения и границ доверительных интервалов.
>> stairs(x,f);
>> hold on;
>> stairs(x,flo,'r:');
>> stairs(x,fup,'r:');

Расчет значений теоретической функции распределения экспоненциального закона с параметром распределения равным 10 и построение общего графика для названных функций.
>> xx = 0:.1:max(t);
>> yy = 1-exp(-xx/10);
>> plot(xx,yy,'g-')
>> grid on
>> hold off;

  В оглавление \ К следующему разделу \ К предыдущему разделу