How to implement large-scale Poisson Regression in Python

Question

I am trying to implement a Poisson Regression in Python to predict rates. I am dealing with a ton of data (too much to store in a DataFrame), which means that using the stan