query: tag:statistics

Control Variate
モンテカルロ求積法などでVariance reduction
の為に利用されるテクニックです。

統計量[math]m[/math]が与えられた時、期待値は以下のように表します。

[math]
\mu=E[m]
[/math]

この時、期待値が[math]E[t]=\tau[/math]であり、
相関係数(correlation coefficient)が
[math]\rho_{mt}=Corr[m,t][/math]であるような統計量
[math]t[/math]があるとすると、以下の式で表される

[math]
m^\star=m-C(t-\tau)
[/math]

は、[math]\mu[/math]に対して不偏(unbiased)です。Cは任意の定数です。

ここで、
[math]\sigma_m[/math], [math]\sigma_t[/math]をそれぞれ、
[math]m[/math], [math]t[/math]の標準偏差(standard deviation)
とした時、定数Cを以下

[math]
C = \frac{\sigma_m}{\sigma_t}\rho_{mt}
[/math]

のように選ぶと、
[math]m^\star[/math]の分散は最小化され、以下のようになります。

[math]
V[m^\star]=(1-\rho_{mt}^2)V[m]
[/math]

元の統計量[math]m[/math]の分散に対して、
[math]\rho_{mt}^2[/math]の分だけ相対的に小さくなります。

すなわち、[math]m[/math]と高い相関関係
(正の相関でも負の相関でも良い)にある
[math]t[/math]が既知であるときに、[math]t[/math]を使って
[math]m^\star[/math]の期待値の推定量を求めることで、
[math]m[/math]の期待値の推定量を効率よく推定することができます。

See Also

posted by genki genki on Sat 19 Jul 2008 at 02:08 with 0 comments