Auf dem Weg zu einer robusten Inferenzmethode für das Galaxienbispektrum: Wahrscheinlichkeitsfunktion und Modellauswahl