Форумы GIS-Lab.info

Здравствуйте. Просветите, пожалуйста, какая в R базовая и/или пакетная возможность оценить степень
соответствия эмпирического распределения распределению Вейбулла-Гнеденко?

nickleb, это оно?

Критерий Колмогорова-Смирнова, после определения параметров методом макс. правдоподобия.

antonv писал(а): ↑18 окт 2018, 14:29nickleb, это оно?

antonv, "очень даже - может быть!" спасибо!

gamm писал(а): ↑18 окт 2018, 14:37Критерий Колмогорова-Смирнова, после определения параметров методом макс. правдоподобия.

gamm, спасибо! кто-то из русскоязычных в известных изданиях об использовании этого критерия (Вейбулла-Гнеденко) на практике, м.б., написал?- почитать бы!

nickleb писал(а): ↑18 окт 2018, 14:50 использовании этого критерия (Вейбулла-Гнеденко) на практике

критерий - это Колмогорова-Смирнова, он сравнивает выборку с заданным распределением (любым).

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

Если это выборка готовых времен до отказа, то параметры можно определить используя fitdistr() из пакета MASS, fitdistr(my.data, densfun="weibull", lower = 0), https://stats.stackexchange.com/questio ... speed-data - это даст нам параметры fitted.shape, fitted.scale

Чтобы проверить, насколько хорошо определились параметры, используем ks.test(x, "pweibull", fitted.shape, fitted.scale)

gamm писал(а): ↑18 окт 2018, 18:44 https://stats.stackexchange.com/questio ... speed-data - это даст нам параметры fitted.shape, fitted.scale

- вроде бы вот в таком контексте мне и надо смотреть... - gamm, спасибо за просвещение, за уделённое время!

мне результат fitdistr() не всегда нравится, я сам пишу функцию правдоподобия, и вставляю ее в optim(). Начальное значение формы можно прикинуть по картинке (сравнив с гистограммой), а масштаб определить из среднего, см. вики https://ru.wikipedia.org/wiki/%D0%A0%D0 ... 0%BB%D0%B0, там есть формулы для MLE оценки параметров.

gamm писал(а): ↑19 окт 2018, 07:04мне результат fitdistr() не всегда нравится, я сам пишу функцию правдоподобия, и вставляю ее в optim(). Начальное значение формы можно прикинуть по картинке (сравнив с гистограммой), а масштаб определить из среднего, см. вики https://ru.wikipedia.org/wiki/%D0%A0%D0 ... 0%BB%D0%B0, там есть формулы для MLE оценки параметров.

gamm, спасибо!