Dear all,

I'm running multiple regression with the same set of controls. I use "local" command so that I don't have to repeat all control variables in each regression.

My code as follows:


local controls age female educ wealth income children province_id commune_id village_id enum_id

eststo clear
probit smoking risk1 `controls' if consistent==1, r cluster(village_id)
margins, dydx(_all) post
est store r1

probit smoking risk2 `controls' if consistent==1,r cluster(village_id)
margins, dydx(_all) post
est store r2

probit smoking risk3 `controls' if consistent==1,r cluster(village_id)
margins, dydx(_all) post
est store r3


However, the results I got are the one without controls.

Can you suggest me other ways to group the same set of control variables?

Thank you very much