Estimation et calibration MCMC du modèle de volatilité stochastique

Présentation du Modèle

L'analyse du modèle se fait au travers d'une représentation espace-état, ayant comme processus observé (y) les rendements financiers et comme processus caché la log-volatilité.
Le but est d'établir une méthode d'estimation de la variable latente et de réussir à calibrer le modèle.
Pour se faire, nous utiliserons les techniques MCMC notamment l'échantillonneur de Gibbs et l'algorithme de Metropolis-Hasting (MH). Ceux-ci ont été présentés auparavant sur ce blog (en mars 2013).

Echantillonneur de Gibbs

On introduit le modèle suivant :

avec x_0:T  le processus de log-volatilité, y_1:T  les données observées et  θ ={α, φ, σ} les paramètres du modèles.

Rappel

L'algorithme de Gibbs permet de générer un échantillon issu d’une distribution jointe à partir de la distribution conditionnelle de chaque variable sachant les autres paramètres et les données observées. Pour plus de détails, se référer au précédent post introduisant cette technique.

Lois conditionnelles

Avant de nous attacher à l’étude des lois conditionnelles nous écrivons la vraisemblance du modèle.

L’intégration souffre d’un problème de dimension (T), nous utilisions alors le principe d’augmentation de données pour résoudre cela (c.f. Robert & Casella, 2004).
La vraisemblance complète est donnée par l’écriture suivante :
                                               

Une fois la vraisemblance définie nous spécifions les lois a priori du modèle :

α ~ N(μ_α ¸σ_α)

φ ~ N(μ_φ¸σ_φ)

                                                               σ_η ~ σ_η^-1   

Nous pouvons maintenant établir les lois conditionnelles pour tous les paramètres et le processus latent x.

Loi conditionnelle pour α

Loi conditionnelle pour φ 

Loi conditionnelle pour σ_η

Loi conditionnelle pour x_2:T-1

Pour filtrer la variable latente nous utilisons l’algorithme MH. En effet, comme le processus est markovien de premier ordre nous devons savoir si la marche proposée améliore la vraisemblance du modèle. Nous proposons la loi suivante pour la loi q.

Loi conditionnelle pour  x₁  et x_T

Application #1

Nous appliquons l'algorithme décrit ci-dessus sur un jeu de données simulées de 500 observations générées avec  α = -0.8,  φ = 0.9,  σ_η = 0.3. Nous réalisons 30000 tirages MCMC, mais 10000 tirages sont largement suffisant.

Application #2

On considère maintenant un échantillon de 2000  données simulées avec les mêmes paramètres.

Application #3

On considère maintenant un échantillon de données provenant des performances journalières du S&P500.  Période: 03/01/2000 - 14/01/2009

Les résultats ci-dessous sont établis à partir de 30000 tirages MCMC.

Le deuxième graphique laisse apparaître la variable : log(return^2)+1.27. Cela permet d'approximer le processus latent (log-vol). En effet, en reprenant l’équation du processus observé on a :

log y_t² =  x_t + log ε_t² , où  log ε_t²  ~ X₁²

La loi peut être approximée  par une gaussienne de paramètre  E(log ε_t²) = -1.27 et V(log ε_t²) = 4.9

Donc : log y_t² - E(log ε_t²)  = x_t + N(0,V(log ε_t²))