使用具有相同比例参数的Gamma变量很容易产生具有Dirichlet分布的随机变量。如果:
然后:
问题 如果比例参数不相等会怎样?
那么这个变量的分布是什么?
对我来说,知道这种分布的期望值就足够了。
我需要一个可以由计算机非常快速地求值的近似封闭代数公式。
假设精度为0.01就足够了。
您可以假设:
注意简而言之,任务是找到该积分的近似值:
使用具有相同比例参数的Gamma变量很容易产生具有Dirichlet分布的随机变量。如果:
然后:
问题 如果比例参数不相等会怎样?
那么这个变量的分布是什么?
对我来说,知道这种分布的期望值就足够了。
我需要一个可以由计算机非常快速地求值的近似封闭代数公式。
假设精度为0.01就足够了。
您可以假设:
注意简而言之,任务是找到该积分的近似值:
Answers:
只是一个初步的说法,如果要计算速度,通常必须牺牲准确性。通常,“更多准确性” =“更多时间”。无论如何,这里是一个二阶近似值,应该改进您在上面的评论中建议的“粗略”近似值:
EDIT An explanation for the above expansion was requested. The short answer is wikipedia. The long answer is given below.
write . Now we need all the "second order" derivatives of . The first order derivatives will "cancel" because they will all involve multiples and which are both zero when taking expectations.
And so the taylor series up to second order is given by:
Taking expectations yields:
Which is the answer I gave. (although I initially forgot the minus sign in the second term)