1. regress logg1018 a2000 labor2 children retire hh_income house houseloan e2002 c7001 savings a2003 age a2012 a2015 a2024 f2021 work a3003 workyear worktime f1001 f2001 f3001 cinsurance region a2022rural[pweight=swgt]
2. xi: regress logg1018 a2000 labor2 children retire hh_income house houseloan e2002 c7001 savings a2003 age a2012 a2015 a2024 f2021 work a3003 workyear worktime f1001 f2001 f3001 cinsurance region i.rural*a2022[pweight=swgt]
3. regress logg1018 a2000 labor2 children retire hh_income house houseloan e2002 c7001 savings a2003 age a2012 a2015 a2024 f2021 work a3003 workyear worktime f1001 f2001 f3001 cinsurance region a2022#rural[pweight=swgt]