Гамма-распределение Опции

[UNDEFEAT]

Добрый день всем посвящённым в статистику

Я не посвящён, потому был бы очень благодарен, если кто-то мне поможет с решением данной проблемы.
Есть случайная величина, которая описывается гамма-распределением с неизвестными параметрами. Есть выборка значений этой СВ. Задача:
1. Найти матожидание.
2. В некой информативной форме указать характеристическое отклонение от матожидания.
Я могу найти параметры k и theta распределения, откуда я автоматически получаю среднее:
EX = k * theta
И дисперсию:
DX = k * theta * theta
Также я могу получить среднеквадратичное отклонение:
sigma = sqrt(DX)
Но оно совершенно неинформативно, потому что иногда оно больше среднего значения, потому могут получаться записи вида:
10 ± 12
Но гамма-распределение оперирует случайными величинами, которые всегда положительны. Оно несимметрично, потому "среднее отклонение влево" (если так можно выразиться) будет меньше "среднего отклонения вправо". Хотелось бы найти такое р, для которого будет простая формула для двух величин dl и dr таких, что:
1) Вероятность того, что значение СВ лежит на отрезке [EX - dl; EX] равна p.
2) Вероятность того, что значение СВ лежит на отрезке [EX; EX + dr] равна p.

Может кто уже сталкивался с такой проблемой? Моя первая идея была найти p при dr = sigma, а уже отсюда найти dl, но как-то пока не получается

Сообщение было отредактировано UNDEFEAT: 8.5.2023, 21:24

[alan]

Я со статистикой имел много дела, но уже давно.
Потому мое неэкспертное мнение:

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

С дисперсией сложнее. Дисперсия это характеристика "ширины" распределения по сути, и обычно она считается для определения доверительного интервала. Но вот если величина распределена не по Гаусу то это не лучший способ определить доверительный интервал.
В данном случае для theta > 1 дисперсия обязана быть больше среднего, так что 10 ± 12 абсолютно адекватная штука. Вопрос в том как ее интерпретировать.

Тут главный вопрос - зачем тебе нужно знать матожидание и "характеристическое отклонение от матожидания"?

Хотелось бы найти такое р, для которого будет простая формула для двух величин dl и dr таких, что:
1) Вероятность того, что значение СВ лежит на отрезке [EX - dl; EX] равна p.
2) Вероятность того, что значение СВ лежит на отрезке [EX; EX + dr] равна p.

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

[UNDEFEAT]

С дисперсией сложнее. Дисперсия это характеристика "ширины" распределения по сути, и обычно она считается для определения доверительного интервала.
...
В данном случае для theta > 1 дисперсия обязана быть больше среднего, так что 10 ± 12 абсолютно адекватная штука. Вопрос в том как ее интерпретировать.

Мне кажется, ты немного путаешь дисперсию, среднеквадратическое отклонение и ширину интервала, в который значения СВ будут попадать с конкретной вероятностью.
Ну и в примере "10 ± 12" второе число - это не дисперсия, а среднеквадратичное отклонение, sigma, корень квадратный из дисперсии.

Тут главный вопрос - зачем тебе нужно знать матожидание и "характеристическое отклонение от матожидания"?

Чтобы сравнить несколько выборок для разных СВ и выбрать самую подходящую из них. Например, если две будут иметь одинаковое матожидание, то я выберу ту, у которой меньший разброс.

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

Я понимаю, что я могу выбрать любое понравившееся мне р, а потом "подобрать" соответствующий интервал, на котором значения двух описанных выше интегралов будут равны р. Собственно, программно я это всё и делаю.

Но это запасной вариант. Я надеялся на то, что кому-то будет известна простая замкнутая форма для некого удобного значения р

Сообщение было отредактировано UNDEFEAT: 8.5.2023, 21:21