Expand Multilevel Logistic Regression Model to include Individual Covariates

Hi Jesse, thank you for your reply. It is very helpful for a PyMC beginner. Now, I have better understanding of the dimension and math code. However, I ran into a problem when I modify the matrix (X) from 2-dims to 1-dim. Although I can produce a 1-dim array of X, I found the dimension of coords is incorrect. I read through the dimensions tutorial in PyMC but I can’t figure out the answer. Do you mind shedding some light for me?