Da wir nun unterschiedliche Modelle kennen gelernt haben, beschäftigen wir uns erstmal mit der Frage, wie wir Modelle Bayesianische 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 Normalverteilungssannahme im Regressionsmodell weggehen, z.B. zur Poisson-Regression. Wir lernen dabei den Metropolis-Sampler kennen:
Gibbs- und Metropolis-Samler 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. Einfach ist es jedoch, ein darauf spezialisiertes Softwarepaket zu benutzen wie
Zum Abschluss ein
Zwei Beispiele sollen die Funktionsweise der Algorithmen nochmal klar machen:
Der folgenden Abschnitte geht noch etwas genauer auf die Details der Monte Carlo-Verfahren und des Ziehen von Zufallszahlen ein: