Variable antithétique

Les variables antithétiques sont une des techniques de réduction de la variance employées dans la méthode de Monte-Carlo. Il s'agit de tirer parti de certaines symétries d'une distribution et de la corrélation négative entre deux variables aléatoires.

Sommaire

1 Principe
2 Exemple : estimation d'une intégrale
3 Notes et références
- 3.1 Bibliographie
- 3.2 Liens externes

Principe

On souhaite estimer $θ = E(h (X))$ . La méthode de Monte-Carlo de base consiste à simuler n variables iid selon la loi de X, disons $X_1, X_2, \cdots, X_n$ , puis à estimer $θ$ par

$\hat \theta = \frac{1}{n} \sum_{i=1}^n h(X_i)$ .

On peut avoir une idée de l'erreur commise en construisant un intervalle de confiance ; ce dernier nécessite un estimateur de la variance de l'estimateur $\sigma^2_{\hat \theta}$ .

Supposons que l'on dispose de deux échantillons de taille n ; le premier est noté $X_1, X_2, \cdots, X_n$ et le second $X'_1, X'_2, \cdots, X'_n$ . Pour simplifier les notations, on pose $m 1, m 2$ les estimateurs empiriques de l'espérance de $h (X)$ sur respectivement l'échantillon 1 et 2. Autrement dit, on aura

$m_1 = \frac{h(X_1) + \cdots + h(X_n)}{n}$

$m_2 = \frac{h(X'_1) + \cdots + h(X'_n)}{n}$ .

L'estimateur Monte-Carlo sur l'échantillon complet est simplement

$\hat \theta = \frac{m_1 + m_2}{2}$

et, du point de vue de la variance :

$\sigma^2_{\hat \theta} = \frac{\sigma^2_{m_1} + \sigma^2_{m_2} + 2 \mathrm{Cov}(m_1, m_2)}{4}$ .

Dans le cas iid, la covariance s'annule et $\sigma^2_{m_1}=\sigma^2_{m_2}$ , si bien que $\sigma^2_{\hat \theta} = \sigma^2_{m_1}/2$ : le facteur 2 s'explique car on a doublé la taille de l'échantillon.

La technique de la variable antithétique consiste à choisir l'échantillon 2 identiquement distribué selon la loi de X mais en renonçant à l'indépendance, plus précisément en s'arrangeant pour que $Cov(m 1, m 2) < 0$ . Il faut donc exploiter les éléments de symétrie de la loi de X afin de construire le second échantillon à partir du premier, en s'assurant de la négativité de la covariance. Se faisant, la variance sera inférieure à la variance "normale" $\sigma^2_{m_1}/2$ .

Par exemple, si la loi de X est la loi uniforme sur $[0;1]$ , le premier échantillon sera simplement $u_1, \cdots, u_n$ , où pour tout i, $u i$ est tirée selon $U (0;1)$ . On construit le second échantillon $u'_1, \cdots, u'_n$ , en posant pour tout i: $u' i = 1 - u i$ . Si les $u 1$ sont uniformes sur $[0;1]$ , alors il en va de même pour les $u' i$ . De plus, la covariance est négative, ce qui permet de réduire la variance initiale.

Un autre exemple concerne la loi normale $N (μ, s)$ . En appliquant la transformation $x' i = 2μ - x i$ , où $x i \sim N (μ, s)$ , on obtient un tirage dans $N (μ, s)$ , qui est négativement corrélé avec le premier tirage $x i$

Exemple : estimation d'une intégrale

On souhaite estimer

$I = \int_0^1 \frac{1}{1+x} \, \mathrm{d}x$ .

La valeur exacte est $I=\ln 2 \approx 0,69314718055995$ . Cette intégrale peut se voir comme l'espérance de $f (U)$ , où

$f(x) = \frac{1}{1+x}$

et U distribuée selon une loi uniforme sur $[0;1]$ .

On compare l'estimateur Monte-Carlo classique (échantillon de taille $2 n$ , avec $n = 1500$ , tiré selon la loi uniforme standard) à l'estimateur avec variable antithétique (échantillon de taille n, complété par l'échantillon transformé $1 - u i$ ). La variance se réduit comme suit

	Estimation	Variance
Méthode classique	0,69365	0,02005
Variable antithétique	0,69399	0,00063

On constate une très nette réduction de la variance dans le cas de l'utilisation d'une variable antithétique.

Notes et références

Bibliographie

M. Hammersley et K. W. Morton, A new Monte Carlo technique antithetic variates, Proc. Camb. Phil. Soc., 52, 449, 1956.

Liens externes

[PDF] un cours sur la réduction de la variance

Portail des probabilités et des statistiques
Portail de la physique
Portail de la finance

Catégories :

Wikimedia Foundation. 2010.

Contenu soumis à la licence CC-BY-SA. Source : Article Variable antithétique de Wikipédia en français (auteurs)

Игры ⚽ Нужно сделать НИР?

Regardez d'autres dictionnaires:

Variable de contrôle — Pour la méthode de Monte Carlo, une variable de contrôle peut être utilisée afin d obtenir une réduction de la variance, en exploitant la corrélation entre plusieurs statistiques. Sommaire 1 Exposé du principe 2 Exemple 3 Notes et références … Wikipédia en Français
Projet:Mathématiques/Liste des articles de mathématiques — Cette page n est plus mise à jour depuis l arrêt de DumZiBoT. Pour demander sa remise en service, faire une requête sur WP:RBOT Cette page recense les articles relatifs aux mathématiques, qui sont liés aux portails de mathématiques, géométrie ou… … Wikipédia en Français
Liste des articles de mathematiques — Projet:Mathématiques/Liste des articles de mathématiques Cette page recense les articles relatifs aux mathématiques, qui sont liés aux portails de mathématiques, géométrie ou probabilités et statistiques via l un des trois bandeaux suivants … Wikipédia en Français
Methode de Monte-Carlo — Méthode de Monte Carlo Pour les articles homonymes, voir Monte Carlo (homonymie). On appelle méthode de Monte Carlo toute méthode visant à calculer une valeur numérique, et utilisant des procédés aléatoires, c est à dire des techniques… … Wikipédia en Français
Méthode De Monte-Carlo — Pour les articles homonymes, voir Monte Carlo (homonymie). On appelle méthode de Monte Carlo toute méthode visant à calculer une valeur numérique, et utilisant des procédés aléatoires, c est à dire des techniques probabilistes. Le nom de ces… … Wikipédia en Français
Méthode de Monte Carlo — Pour les articles homonymes, voir Monte Carlo (homonymie). On appelle méthode de Monte Carlo toute méthode visant à calculer une valeur numérique, et utilisant des procédés aléatoires, c est à dire des techniques probabilistes. Le nom de ces… … Wikipédia en Français
Méthode de monte-carlo — Pour les articles homonymes, voir Monte Carlo (homonymie). On appelle méthode de Monte Carlo toute méthode visant à calculer une valeur numérique, et utilisant des procédés aléatoires, c est à dire des techniques probabilistes. Le nom de ces… … Wikipédia en Français
Réduction de la variance — La réduction de la variance regroupe l ensemble des techniques, plus ou moins simples, qui permettent de réduire la variance des estimateurs de Monte Carlo. En voici une courte liste : Variable antithétique : on introduit une seconde… … Wikipédia en Français
Simulation de Monte Carlo — Méthode de Monte Carlo Pour les articles homonymes, voir Monte Carlo (homonymie). On appelle méthode de Monte Carlo toute méthode visant à calculer une valeur numérique, et utilisant des procédés aléatoires, c est à dire des techniques… … Wikipédia en Français
Méthode de Monte-Carlo — Pour les articles homonymes, voir Monte Carlo (homonymie). Le terme méthode de Monte Carlo, ou méthode Monte Carlo, désigne toute méthode visant à calculer une valeur numérique en utilisant des procédés aléatoires, c est à dire des techniques… … Wikipédia en Français

Dictionnaires et Encyclopédies sur 'Academic'

Variable antithétique

Sommaire

Principe

Exemple : estimation d'une intégrale