Как я могу рассчитать параметры и для бета-распределения, если я знаю среднее значение и дисперсию, которые я хочу иметь в распределении? Примеры команды R для этого были бы наиболее полезны.β
r
distributions
estimation
beta-distribution
Дэйв Кинкейд
источник
источник
Ответы:
Я устанавливаю и и решено для и . Мои результаты показывают, что и σ2=αβ
Я написал некоторый R-код для оценки параметров распределения Beta по заданному среднему значению mu и дисперсии var:
Там было некоторое замешательство вокруг границ и для любого данного бета-дистрибутива, поэтому давайте проясним это здесь.σ 2μ σ2
источник
estBetaParams(0.06657, 0.1)
яalpha=-0.025
,beta=-0.35
. Как это возможно?Вот общий способ решения этих типов проблем, используя Maple вместо R. Это работает и для других дистрибутивов:
что приводит к решению
Это эквивалентно решению Макса.
источник
В R бета-распределение с параметрами и имеет плотностьshape1=a shape2=b
для , и .a>0 b>0 0<x<1
В R вы можете вычислить это
В этой параметризации среднее значение а дисперсия . Итак, теперь вы можете следить за ответом Ника Сабби.E(X)=aa+b V(X)=ab(a+b)2(a+b+1)
Хорошая работа!
редактировать
Я нахожу:
а также
где и .V = V ( X )μ=E(X) V=V(X)
источник
Например, в Википедии вы можете найти следующие формулы для среднего значения и дисперсии бета-распределения с учетом альфа и бета: и Обращение этих (заполните в нижнем уравнении) должно дать вам результат, который вы хотите (хотя это может занять некоторую работу).
источник
Для обобщенного бета-распределения, определенного на интервале , у вас есть отношения:[a,b]
который можно перевернуть, чтобы получить:
где
источник
Решите уравнение для или , решив для , вы получите Затем вставьте это во второе уравнение и решите для . Таким образом, вы получаете что упрощает до Тогда закончите решение для .μ α β β β=α(1−μ)μ α σ2=α2(1−μ)μ(α+α(1−μ)μ)2(α+α(1−μ)μ+1) σ2=α2(1−μ)μ(αμ)2α+μμ σ2=(1−μ)μ2α+μ α
источник
Я искал питона, но наткнулся на это. Так что это было бы полезно для таких, как я.
Вот код Python для оценки бета-параметров (согласно приведенным выше уравнениям):
Вы можете проверить параметры и , импортировав пакет.α β
scipy.stats.beta
источник