Any way to efficiently run Poisson regression thousands of times?

Thanks @Dan , the initial guess also helps!