Jump to content

Распределение Вейбулла

(Перенаправлено из диаграммы Вейбулла )
Вейбулл (2-параметрический)
Функция плотности вероятности
Функция распределения вероятностей
Кумулятивная функция распределения
Кумулятивная функция распределения
Параметры шкала
форма
Поддерживать
PDF
CDF
Квантиль
Иметь в виду
медиана
Режим
Дисперсия
асимметрия
Избыточный эксцесс (см. текст)
Энтропия
МГФ
CF
Расхождение Кульбака – Лейблера см. ниже

В теории вероятностей и статистике / распределение Вейбулла ˈ w b ʊ l / является непрерывным распределением вероятностей . Он моделирует широкий диапазон случайных величин, в основном таких как время до отказа или время между событиями. Примерами являются максимальное количество осадков за один день и время, которое пользователь проводит на веб-странице.

Распределение названо в честь шведского математика Валодди Вейбулла , подробно описавшего его в 1939 году. [1] хотя он был впервые идентифицирован Рене Морисом Фреше и впервые применен Розином и Раммлером (1933) для описания распределения частиц по размерам .

Определение

[ редактировать ]

Стандартная параметризация

[ редактировать ]

Функция плотности вероятности Вейбулла случайной величины равна [2] [3]

где k > 0 — параметр формы , а λ > 0 — параметр масштаба распределения. Его дополнительная кумулятивная функция распределения представляет собой растянутую экспоненциальную функцию . Распределение Вейбулла связано с рядом других вероятностных распределений; в частности, он интерполирует между экспоненциальным распределением ( k = 1) и распределением Рэлея ( k = 2 и [4] ).

Если величина X представляет собой «время до отказа», распределение Вейбулла дает распределение, для которого интенсивность отказов пропорциональна степени времени. Параметр формы k равен этой степени плюс единица, поэтому этот параметр можно интерпретировать непосредственно следующим образом: [5]

  • Значение указывает на то, что интенсивность отказов со временем уменьшается (как в случае эффекта Линди , который, однако, соответствует распределениям Парето [6] а не распределения Вейбулла). Это происходит, если наблюдается значительная «детская смертность» или дефектные изделия выходят из строя рано, а частота отказов снижается с течением времени по мере того, как дефектные изделия отсеиваются из популяции. В контексте диффузии инноваций это означает негативную молву: функция риска представляет собой монотонно убывающую функцию доли усыновителей;
  • Значение указывает на то, что интенсивность отказов постоянна во времени. Это может означать, что случайные внешние события вызывают смертность или неудачу. Распределение Вейбулла сводится к экспоненциальному распределению;
  • Значение указывает на то, что интенсивность отказов увеличивается со временем. Это происходит, если существует процесс «старения» или детали, которые с большей вероятностью выходят из строя с течением времени. В контексте распространения инноваций это означает положительную молву: функция риска представляет собой монотонно возрастающую функцию доли последователей. Функция сначала выпуклая, затем вогнутая с точкой перегиба при .

В области материаловедения параметр формы распределения сил известен как модуль Вейбулла . В контексте распространения инноваций распределение Вейбулла представляет собой «чистую» модель имитации/отвержения.

Альтернативные параметризации

[ редактировать ]

Первая альтернатива

[ редактировать ]

Приложения в медицинской статистике и эконометрике часто используют другую параметризацию. [7] [8] Параметр формы k такой же, как указано выше, а параметр масштаба равен . В этом случае для x ≥ 0 функция плотности вероятности равна

кумулятивная функция распределения

функция квантиля

функция опасности

и среднее значение

Вторая альтернатива

[ редактировать ]

Также можно найти вторую альтернативную параметризацию. [9] [10] Параметр формы k такой же, как и в стандартном случае, а параметр масштаба λ заменяется параметром скорости β = 1/ λ . Тогда для x ≥ 0 функция плотности вероятности равна

кумулятивная функция распределения

функция квантиля

а функция опасности равна

Во всех трех параметризациях опасность уменьшается при k < 1, увеличивается при k > 1 и остается постоянной при k = 1, и в этом случае распределение Вейбулла сводится к экспоненциальному распределению.

Характеристики

[ редактировать ]

Функция плотности

[ редактировать ]

Вид функции плотности распределения Вейбулла резко меняется с изменением значения k . При 0 < k < 1 функция плотности стремится к ∞ при приближении x к нулю сверху и строго убывает. При k = 1 функция плотности стремится к 1/ λ при приближении x к нулю сверху и строго убывает. При k > 1 функция плотности стремится к нулю по мере приближения x к нулю сверху, возрастает до своей моды и убывает после нее. Функция плотности имеет бесконечный отрицательный наклон при x = 0, если 0 < k < 1, бесконечный положительный наклон при x = 0, если 1 < k < 2, и нулевой наклон при x = 0, если k > 2. Для k = 1 плотность имеет конечный отрицательный наклон при x = 0. Для k = 2 плотность имеет конечный положительный наклон при x = 0. Когда k стремится к бесконечности, распределение Вейбулла сходится к дельта-распределению Дирака с центром в точке x = λ. При этом асимметрия и коэффициент вариации зависят только от параметра формы. Обобщением распределения Вейбулла является гиперболатическое распределение типа III .

Кумулятивная функция распределения

[ редактировать ]

Кумулятивная функция распределения для распределения Вейбулла равна

для x ≥ 0 и F ( x ; k ; λ) = 0 для x < 0.

Если x = λ, то F ( x ; k ; λ) = 1 − e −1 ≈ 0,632 для всех значений k . И наоборот: при F ( x ; k ; λ ) = 0,632 значение x λ .

Функция квантиля (обратного кумулятивного распределения) для распределения Вейбулла:

для 0 ≤ p < 1.

Интенсивность отказов h (или функция опасности) определяется выражением

Среднее время наработки на отказ составляет

логарифма Производящая момент функция случайной распределенной Вейбулла величины определяется выражением [11]

где Γ гамма-функция . Аналогично, характеристическая функция журнала X определяется выражением

В частности, n-й момент X определяется исходный выражением

и Среднее значение дисперсия случайной Вейбулла величины можно выразить как

и

Асимметрия определяется

где , что также можно записать как

где среднее значение обозначается µ , а стандартное отклонение обозначается σ .

Избыточный эксцесс определяется выражением

где . Избыток эксцесса также можно записать как:

Функция генерации момента

[ редактировать ]

доступно множество выражений Для производящей функции момента самого X . В качестве степенного ряда , поскольку исходные моменты уже известны, мы имеем

Альтернативно можно попытаться разобраться непосредственно с интегралом

Если параметр k предполагается рациональным числом, выраженным как k = p / q , где p и q — целые числа, то этот интеграл можно вычислить аналитически. [12] Заменив t на − t , можно найти

где G G-функция Мейера .

Характеристическая функция также была получена Muraleedharan et al. (2007) . Характеристическая функция и функция, генерирующая момент трехпараметрического распределения Вейбулла, также были получены Соаресом (2014) прямого подхода.

Минимумы

[ редактировать ]

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

То есть, также будет распределено Вейбуллом с параметром масштаба и с параметром формы .

Трюки с репараметризацией

[ редактировать ]

Исправьте некоторые . Позволять неотрицательны и не все равны нулю, и пусть быть независимыми образцами , затем [13]

  • .

Энтропия Шеннона

[ редактировать ]

Информационная энтропия определяется выражением

где постоянная Эйлера–Машерони . Распределение Вейбулла — это максимальное распределение энтропии для неотрицательной действительной случайной величины с фиксированным ожидаемым значением x . к равный λ к и фиксированное ожидаемое значение ln( x к ), равный ln( λ к ) −  .

Расхождение Кульбака – Лейблера

[ редактировать ]

Расхождение Кульбака – Лейблера между двумя распределениями Вейбулла определяется выражением [14]

Оценка параметров

[ редактировать ]

Обычный метод наименьших квадратов с использованием графика Вейбулла

[ редактировать ]
сюжет Вейбулла

Соответствие распределения Вейбулла данным можно оценить визуально с помощью графика Вейбулла. [15] График Вейбулла представляет собой график эмпирической кумулятивной функции распределения. данных по специальным осям в виде графика Q–Q . Оси против . Причиной такой замены переменных является то, что кумулятивную функцию распределения можно линеаризовать:

которая, как можно видеть, имеет стандартную форму прямой линии. Следовательно, если данные получены из распределения Вейбулла, то на графике Вейбулла ожидается прямая линия.

Существуют различные подходы к получению эмпирической функции распределения по данным: один метод заключается в получении вертикальной координаты для каждой точки с помощью где - это ранг точки данных и количество точек данных. [16] [17]

Линейную регрессию также можно использовать для численной оценки согласия и оценки параметров распределения Вейбулла. Градиент напрямую информирует о параметре формы. и параметр масштаба также можно сделать вывод.

Метод моментов

[ редактировать ]

Коэффициент вариации распределения Вейбулла зависит только от параметра формы: [18]

Уравнивание объемов выборки к , моментная оценка параметра формы можно прочитать либо из справочной таблицы, либо из графика против . Более точная оценка можно найти с помощью алгоритма поиска корня для решения

Затем оценку момента масштабного параметра можно найти с использованием первого уравнения момента как

Максимальная вероятность

[ редактировать ]

Оценка максимального правдоподобия для задан параметр является [18]

Оценка максимального правдоподобия для является решением для k следующего уравнения [19]

Это уравнение определяет только неявно, обычно нужно найти решение для числовыми средствами.

Когда являются крупнейшие наблюдаемые выборки из набора данных объемом более выборок, то оценка максимального правдоподобия для задан параметр является [19]

Кроме того, при этом условии оценка максимального правдоподобия для является [ нужна ссылка ]

Опять же, поскольку это неявная функция, обычно нужно найти числовыми средствами.

Приложения

[ редактировать ]

Используется распределение Вейбулла. [ нужна ссылка ]

Подобрано кумулятивное распределение Вейбулла для максимального количества осадков за один день с помощью CumFreq , см. также подбор распределения. [20]
Подогнанные кривые для данных временных рядов добычи нефти [21]
[ редактировать ]
  • Если , то переменная распределен ли Гамбель (минимум) с параметром местоположения и параметр масштабирования . То есть, .
  • Распределение Вейбулла — это обобщенное гамма-распределение, оба параметра формы которого равны k .
  • Транслированное распределение Вейбулла (или трехпараметрическое распределение Вейбулла) содержит дополнительный параметр. [11] Он имеет функцию плотности вероятности

    для и для , где параметр формы , является параметром масштаба и параметр местоположения распределения. Значение устанавливает начальное время безотказной работы перед началом обычного процесса Вейбулла. Когда , это сводится к двухпараметрическому распределению.
  • Распределение Вейбулла можно охарактеризовать как распределение случайной величины такая, что случайная величина

    — стандартное экспоненциальное распределение с интенсивностью 1. [11]
  • Это означает, что распределение Вейбулла также можно охарактеризовать как равномерное распределение : если равномерно распределен по , то случайная величина распределяется ли Вейбулл с параметрами и . Обратите внимание, что здесь эквивалентно чуть выше. Это приводит к легко реализуемой численной схеме моделирования распределения Вейбулла.
  • Распределение Вейбулла интерполирует экспоненциальное распределение с интенсивностью когда и распределение Рэлея моды когда .
  • Распределение Вейбулла (обычно достаточное в технике надежности ) представляет собой частный случай трехпараметрического возведенного в степень распределения Вейбулла , где дополнительный показатель степени равен 1. Возведенное в степень распределение Вейбулла допускает унимодальное . распределение Вейбулла в форме ванны [30] и монотонная интенсивность отказов .
  • Распределение Вейбулла является частным случаем обобщенного распределения экстремальных значений . Именно в связи с этим распределение впервые было идентифицировано Морисом Фреше в 1927 году. [31] Близко родственное распределение Фреше , названное в честь этой работы, имеет функцию плотности вероятности

  • Распределение случайной величины, определяемое как минимум из нескольких случайных величин, каждая из которых имеет различное распределение Вейбулла, является поли-распределением Вейбулла .
  • Распределение Вейбулла было впервые применено Розином и Раммлером (1933) для описания распределения частиц по размерам. Он широко используется в переработке полезных ископаемых для описания распределения частиц по размерам в процессах измельчения . В этом контексте кумулятивное распределение определяется выражением

    где
    • размер частиц
    • представляет собой 80-й процентиль распределения частиц по размерам.
    • – параметр, описывающий разброс распределения
  • Из-за своей доступности в электронных таблицах он также используется там, где основное поведение на самом деле лучше моделируется с помощью распределения Эрланга . [32]
  • Если затем ( Экспоненциальное распределение )
  • При тех же значениях k гамма-распределение принимает схожие формы, но распределение Вейбулла более платикуртное .

См. также

[ редактировать ]
  1. ^ Бауэрс и др. ал. (1997) Актуарная математика, 2-е изд. Общество актуариев.
  2. ^ Папулис, Афанасиос Папулис; Пиллаи, С. Унникришна (2002). Вероятность, случайные величины и случайные процессы (4-е изд.). Бостон: МакГроу-Хилл. ISBN  0-07-366011-6 .
  3. ^ Кизилерсу, Айше; Крир, Маркус; Томас, Энтони В. (2018). «Распределение Вейбулла» . Значение . 15 (2): 10–11. дои : 10.1111/j.1740-9713.2018.01123.x .
  4. ^ «Релеевское распределение – MATLAB и Simulink – MathWorks Australia» . www.mathworks.com.au .
  5. ^ Цзян, Р.; Мурти, DNP (2011). «Исследование параметра формы Вейбулла: свойства и значение». Проектирование надежности и системная безопасность . 96 (12): 1619–26. дои : 10.1016/j.ress.2011.09.003 .
  6. ^ Элиазар, Иддо (ноябрь 2017 г.). «Закон Линди». Физика А: Статистическая механика и ее приложения . 486 : 797–805. Бибкод : 2017PhyA..486..797E . дои : 10.1016/j.physa.2017.05.077 . S2CID   125349686 .
  7. ^ Коллетт, Дэвид (2015). Моделирование данных о выживаемости в медицинских исследованиях (3-е изд.). Бока-Ратон: Чепмен и Холл / CRC. ISBN  978-1439856789 .
  8. ^ Кэмерон, AC; Триведи, ПК (2005). Микроэконометрика: методы и приложения . п. 584. ИСБН  978-0-521-84805-3 .
  9. ^ Калбфляйш, доктор медицинских наук; Прентис, Р.Л. (2002). Статистический анализ данных о времени отказа (2-е изд.). Хобокен, Нью-Джерси: Дж. Уайли. ISBN  978-0-471-36357-6 . OCLC   50124320 .
  10. ^ Терно, Т. (2020). «Пакет для анализа выживания в R». Пакет R версии 3.1.
  11. ^ Перейти обратно: а б с Джонсон, Коц и Балакришнан, 1994 г.
  12. ^ См. ( Cheng, Tellambura & Beaulieu 2004 ) для случая, когда k является целым числом, и ( Sagias & Karagiannidis 2005 ) для рационального случая.
  13. ^ Балог, Матей; Трипуранени, Нилеш; Гахрамани, Зубин; Веллер, Адриан (17 июля 2017 г.). «Потерянные родственники трюка с Гамбелем» . Международная конференция по машинному обучению . ПМЛР: 371–379.
  14. ^ Баукхаге, Кристиан (2013). «Вычисление расхождения Кульбака-Лейблера между двумя распределениями Вейбулла». arXiv : 1310.3713 [ cs.IT ].
  15. ^ «1.3.3.30. График Вейбулла» . www.itl.nist.gov .
  16. ^ Уэйн Нельсон (2004) Прикладной анализ жизненных данных . Уайли-Блэквелл ISBN   0-471-64462-5
  17. ^ Барнетт, В. (1975). «Методы вероятностного построения графиков и статистика заказов» . Журнал Королевского статистического общества. Серия C (Прикладная статистика) . 24 (1): 95–108. дои : 10.2307/2346708 . ISSN   0035-9254 .
  18. ^ Перейти обратно: а б Коэн, А. Клиффорд (ноябрь 1965 г.). «Оценка максимального правдоподобия в распределении Вейбулла на основе полных и цензурированных выборок» (PDF) . Технометрика . 7 (4): 579–588.
  19. ^ Перейти обратно: а б Сорнетт, Д. (2004). Критические явления в естествознании: хаос, фракталы, самоорганизация и беспорядок . .
  20. ^ «CumFreq, Подбор вероятностей распределения, бесплатное программное обеспечение, накопительная частота» .
  21. ^ Перейти обратно: а б Ли, Се Юн; Маллик, Бани (2021). «Байесовское иерархическое моделирование: применение к результатам добычи в сланцах Игл Форд в Южном Техасе» . Санкхья Б. 84 : 1–43. дои : 10.1007/s13571-020-00245-8 .
  22. ^ «Распределение скорости ветра, Вейбулл – REUK.co.uk» . www.reuk.co.uk.
  23. ^ Лю, Чао; Уайт, Райен В.; Дюмэ, Сьюзен (19 июля 2010 г.). Понимание поведения при просмотре веб-страниц посредством анализа времени пребывания по Вейбуллу . АКМ. стр. 379–386. дои : 10.1145/1835449.1835513 . ISBN  9781450301534 . S2CID   12186028 .
  24. ^ Шариф, М. Наваз; Ислам, М.Назрул (1980). «Распределение Вейбулла как общая модель прогнозирования технологических изменений». Технологическое прогнозирование и социальные изменения . 18 (3): 247–56. дои : 10.1016/0040-1625(80)90026-8 .
  25. ^ Вычислительная оптимизация двигателя внутреннего сгорания, стр. 49
  26. ^ Остин, LG; Климпель, Р.Р.; Лаки, ПТ (1984). Технология уменьшения габаритов . Хобокен, Нью-Джерси: ISBN Guinn Printing Inc.  0-89520-421-5 .
  27. ^ Чандрашекар, С. (1943). «Стохастические проблемы физики и астрономии». Обзоры современной физики . 15 (1): 86.
  28. ^ ECSS-E-ST-10-12C – Методы расчета полученной радиации и ее последствий, а также политика расчетных запасов (Отчет). Европейское сотрудничество по космической стандартизации. 15 ноября 2008 г.
  29. ^ Л. Д. Эдмондс; CE Барнс; Л. З. Шейк (май 2000 г.). «8.3 Подгонка кривой». Введение в воздействие космической радиации на микроэлектронику (PDF) (Отчет). Лаборатория реактивного движения НАСА, Калифорнийский технологический институт. стр. 75–76.
  30. ^ «Эволюция систем и надежность систем» . Сысев (Бельгия). 01.01.2010.
  31. ^ Монтгомери, Дуглас (19 июня 2012 г.). Введение в статистический контроль качества . [Sl]: Джон Уайли. п. 95. ИСБН  9781118146811 .
  32. ^ Чатфилд, К.; Гудхардт, Дж.Дж. (1973). «Модель потребительских закупок с Эрлангом времени между покупками». Журнал Американской статистической ассоциации . 68 (344): 828–835. дои : 10.1080/01621459.1973.10481432 .

Библиография

[ редактировать ]
[ редактировать ]
Arc.Ask3.Ru: конец переведенного документа.
Arc.Ask3.Ru
Номер скриншота №: f77e0eb2b8ecdbc2d8dbdc1b255cbc6a__1720113660
URL1:https://arc.ask3.ru/arc/aa/f7/6a/f77e0eb2b8ecdbc2d8dbdc1b255cbc6a.html
Заголовок, (Title) документа по адресу, URL1:
Weibull distribution - Wikipedia
Данный printscreen веб страницы (снимок веб страницы, скриншот веб страницы), визуально-программная копия документа расположенного по адресу URL1 и сохраненная в файл, имеет: квалифицированную, усовершенствованную (подтверждены: метки времени, валидность сертификата), открепленную ЭЦП (приложена к данному файлу), что может быть использовано для подтверждения содержания и факта существования документа в этот момент времени. Права на данный скриншот принадлежат администрации Ask3.ru, использование в качестве доказательства только с письменного разрешения правообладателя скриншота. Администрация Ask3.ru не несет ответственности за информацию размещенную на данном скриншоте. Права на прочие зарегистрированные элементы любого права, изображенные на снимках принадлежат их владельцам. Качество перевода предоставляется как есть. Любые претензии, иски не могут быть предъявлены. Если вы не согласны с любым пунктом перечисленным выше, вы не можете использовать данный сайт и информация размещенную на нем (сайте/странице), немедленно покиньте данный сайт. В случае нарушения любого пункта перечисленного выше, штраф 55! (Пятьдесят пять факториал, Денежную единицу (имеющую самостоятельную стоимость) можете выбрать самостоятельно, выплаичвается товарами в течение 7 дней с момента нарушения.)