In diesem Abschnitt vergleichen wir Modelle mit dem Bayes-Faktor und lernen an der Poisson-Regression den Metropolis- und den Metropolis-Hastings-Algorithmus kennen. Danach schätzen wir dieselben Modelle mit der Software JAGS.
Da wir nun unterschiedliche Modelle kennen gelernt haben, beschäftigen wir uns erstmal mit der Frage, wie wir Modelle Bayesianisch vergleichen können. Das Grundprinzip ist dabei einfach: Das Modell ist ein weiterer uns unbekannter Parameter. Wir müssen uns also die Posterioriwahrscheinlichkeit des Modells ansehen. Das führt uns zum
Mit dem Bayes-Ansatz kann man relativ leicht von der Normalverteilungsannahme im Regressionsmodell weggehen, z.B. zur Poisson-Regression. Wir lernen dabei den Metropolis-Sampler kennen:
Gibbs- und Metropolis-Sampler lassen sich auch kombinieren. Wir wenden die Methode auf ein weiteres Poisson-Modell an:
In den einzelnen Abschnitten sehen wir uns immer die Implementierung “von Hand” an. Einfacher ist es jedoch, ein darauf spezialisiertes Softwarepaket zu benutzen wie
Zum Abschluss ein
Zwei Beispiele sollen die Funktionsweise der Algorithmen nochmal klar machen:
Der folgende Abschnitt geht noch etwas genauer auf die Details der Monte Carlo-Verfahren und des Ziehens von Zufallszahlen ein: