Dear Statalisters,


I'm running a two-level empty mixed model, with levels individuals, city and state. The dependent variable being satisfaction with city services.

Code:
mixed satisfaction   || state: || city:
After that I'm retrieving BLUPs and standard errors for the city level:


Code:
predict eb*, reffects relevel(city)  
gen beta0 = _b[_cons] + eb1
predict eb_se*, reses relevel(city)  
gen eb_se1_l = beta0 - 1.96*eb_se1  //ci
gen eb_se1_u = beta0 + 1.96*eb_se1

What's the best way to test the statistical difference between spesific citities? Is it possible to use OLS, like this:

Code:
regress beta0 i.city
margins city, pwcompare(group sort)

Looking forward to your comments.