Sparse array glm fit in pymc3 very long with NUTS

I would give this a go on AWS. It seems to me like you’ve hit a limitation on your computer. Let me know how it goes.