2380 shaares
marrant
L'histoire est géniale.
Une approche pour bien sélectionner la proposal en MCMC
Dans SAS, la phase de tuning de la proposal dans l'algorithme de Metropolis consiste dans le cas multivarié, à estimer la matrice de covariance d'après les réalisations déjà obtenues. Le principe est le suivant: on a des tuning loops de 500 itérations par défaut chacunes. Lors de chaque loop, on lance une chaîne avec une proposal multinormale caractérisée par une matrice de covariance donnée. À la loop suivante, même principe, mais la matrice de covariance est mise à jour comme une moyenne pondérée de l'ancienne matrice de covariance et de la nouvelle matrice obtenue d'après les réalisations. En plus, un side product, c'est qu'on trouve le MAP de la posterior, dont on peut se servir comme valeurs initiales.
Très intéressant: en MCMC lorsque la proposal (e.g. une gaussienne) n'a pas le même support que la cible (e.g. une variable strictement positive), on est tenté de répéter l'échantillonnage de la proposal jusqu'à ce que la valeur proposée soit dans le support de la cible. Il est bien connu que cette approche n'est pas valide quand elle est implémentée de façon bourrine. Mais l'auteur montre ici qu'elle peut être valide, à condition d'utiliser la bonne distribution dans le calcul du rapport pour déterminer la proba d'acceptation.
M'a l'air super intéressant ce blog...
M'a l'air super intéressant ce blog...
C'est marrant ça. Breiman était profondément anti-bayésien, parce qu'à l'époque (milieu années 1990), les bayésiens étaient surtout des théoriciens qui ne mettaient jamais les mains dans les données. Les pbs réels trop compliqués étaient inaccessibles aux bayésiens... Les choses ont bien changé!
N'empêche, j'attends quand même le jour ou qqn inventera une méthode qui nous permet de nous passer du MCMC et de tous ses pbs de tuning...
N'empêche, j'attends quand même le jour ou qqn inventera une méthode qui nous permet de nous passer du MCMC et de tous ses pbs de tuning...
Très utile. Ya même une formule de politesse pour écrire au pape.
Sinon, yen a qui servent.
Sinon, yen a qui servent.
Ah ben ya une illustration de l'analogie de Gelman sur les modélisateurs tentant d'estimer des effets fins avec des données pourries. Me la stocke celle-là.
Marrant. Enfin, pas pour elle.
Un article pas inintéressant sur l'intérêt d'utiliser dplyr pour la manip de données.
Yen a plusieurs que je ne connaissais pas...
À regarder donc
À regarder donc
Pour presque tout savoir sur la situation politique en turquie.
Marrant
Une explication plus détaillée
Putain ils veulent virer Houiller. Dommage que la pétition ne doive être signée que par le personnel INRA, j'aurais signé direct...
Edit: j'ai signé, apparemment, les signataires non-INRA sont nombreux.
Edit: j'ai signé, apparemment, les signataires non-INRA sont nombreux.
Intéressant.
Oui, j'avais entendu dire ça aussi...
via Mathieu.
via Mathieu.
Bon ben maintenant, de John Kerry aux experts de sciences po, tout le monde est d'accord... Et la temporisation de Cameron, finalement acceptée de tous depuis hier, tend à confirmer le truc.
hmm
Ben décidément, j'aime bien ce gars-là, il est très rare que je sois en désaccord.
Et ça fait du bien d'avoir une voix un peu dissonnante du discours convenu des médias
Et ça fait du bien d'avoir une voix un peu dissonnante du discours convenu des médias