P-002 · 2025 · Teamprojekt

Hierarchisches bayessches Pareto/NBD-Kundenmodell

Replikation von Abe (2009): ein hierarchisch-bayessches Pareto/NBD-Modell, per MCMC auf CDNOW geschätzt, prognostiziert Kundenaktivität, Lebensdauer und Wiederkäufe.

PythonBayesian modellingMCMCCDNOW

Das Modell hinter jedem Kunden

CDNOW-Kaufhistorien: Frequenz, Recency und Kundenalter
KaufprozessPoisson-Käufe mit kundenspezifischer Rate λ
AbwanderungsprozessExponentielle Lebensdauer mit kundenspezifischer Rate μ
Hierarchischer Priorλ und μ gemeinsam gezogen, mit und ohne Umsatzkovariate

Posterior per MCMC geschätzt; Kalibrierung und Validierung verglichen

Erwartete Lebensdauer, Überlebenswahrscheinlichkeit und prognostizierte Transaktionen
Eine Replikation von Abe (2009): Die Hierarchie erlaubt Korrelation zwischen Kauf- und Abwanderungsraten über Kunden hinweg.

Die Frage

Im nicht vertraglichen Handel sieht ein Unternehmen nie, dass ein Kunde abwandert: Er hört einfach auf zu kaufen. Das Pareto/NBD-Modell bildet das mit zwei verborgenen Prozessen je Kunde ab: wie oft er kauft, solange er aktiv ist, und wann er abwandert. Abe (2009) ergänzt einen hierarchischen Prior, sodass Kauf- und Abwanderungsraten korrelieren und mit Kundenmerkmalen verknüpft werden können. Unser Team hat dieses Modell auf dem öffentlichen CDNOW-Datensatz repliziert.

Vorgehen

Der Workflow verdichtet das CDNOW-Transaktionsprotokoll zu einer Zusammenfassung je Kunde: Anzahl der Wiederholungskäufe, Zeitpunkt des letzten Kaufs und Beobachtungsdauer. Danach werden das Standard-Pareto/NBD und zwei hierarchische Varianten per MCMC geschätzt: HB M1 ohne Kovariaten und HB M2 mit dem ersten Kaufbetrag als Kovariate. Jedes Modell wird in einer Kalibrierungs- und einer späteren Validierungsperiode geprüft, je Kunde und aggregiert.

Ergebnisse

Bei den aggregierten wöchentlichen Käufen bildete das Standard-Pareto/NBD die Daten am genauesten ab, mit einem Validierungs-MAPE von 6,54 % gegenüber 11,17 % für HB M1 und 10,16 % für HB M2. Es erreichte auch die höchste Validierungskorrelation auf Kundenebene (0,63 gegenüber 0,56). Die Kovariate verbesserte HB M2 in der Validierung gegenüber HB M1, und die geschätzte Korrelation zwischen Kauf- und Abwanderungsrate blieb nahe null. Für jeden Kunden liefern die Modelle erwartete Lebensdauer, Überlebenswahrscheinlichkeit und erwartete Transaktionen.

Meine Rolle und Grenzen

Ich habe im Team an der Implementierung und Validierung des hierarchischen Modells gearbeitet. Es handelt sich um eine Replikation auf einem Benchmark-Datensatz mit 2.357 Kunden; die Ergebnisse zeigen also das Verhalten der Modelle auf CDNOW, nicht ihre Leistung für einen produktiven Kundenbestand. Mein späterer TabPFN-3.5-Benchmark vergleicht beide Modelle erneut auf größeren Daten.

Aggregierter Prognosefehler: wöchentliche Wiederholungskäufe

Pareto/NBDHB M1HB M2
036912Pareto/NBD, Kalibrierung: 5,48%5,5HB M1, Kalibrierung: 7,10%7,1HB M2, Kalibrierung: 7,17%7,2KalibrierungPareto/NBD, Validierung: 6,54%6,5HB M1, Validierung: 11,17%11,2HB M2, Validierung: 10,16%10,2ValidierungPareto/NBD, Gepoolt: 4,07%4,1HB M1, Gepoolt: 6,81%6,8HB M2, Gepoolt: 6,54%6,5GepooltMAPE in % (niedriger ist besser)
Wöchentlicher Zeitreihen-MAPE auf CDNOW aus der Modellgüte-Tabelle des Repositorys. Das Standard-Pareto/NBD bildet die aggregierten Käufe am genauesten ab, während die Umsatzkovariate (HB M2) den Validierungsfehler gegenüber HB M1 senkt.

Öffentlichen Code und Dokumentation auf GitHub ansehen ↗

Timur Salakhetdinov ·