206 samples :: Mass matrix contains zeros on the diagonal

[AlexAndorra]

Can you try changing the init method? --> pm.sample(init="adapt_diag").
As the GLM module’s priors are very wide, that could be a cause of your error.