Sparse array glm fit in pymc3 very long with NUTS

I think it is also worth trying to split the dataset, or ADVI with minibatch - we dont have a lot of example dealing with large data set so I would love to learn about your experience and approach