Tbf python's stats implementations can be garbage; the last time I checked you can't do multiple levels for hierarchical regression.