* ssc install itsa,replace 
use "/Users/zhangjian/COVID-19-and-AQI-in-China/data_0224.dta",clear
cap drop  xtcomfirm 


tostring year month day,replace
gen date = year + "-" + month + "-" + day
gen Date = date(date, "YMD")
format %tdCCYY-NN-DD Date

duplicates drop city date,force
encode city, gen(CITY) 

keep if  Date > date("2019-12-01","YMD")

sort CITY Date
xtset CITY Date
bys CITY : gen t = _n
tsset CITY t

* 广州53 重庆28 武汉18 成都14 长沙31。合肥8    南京5。 南昌7 杭州17  西藏2 西宁24 银川29
* 西安25 郑州27
itsa pm10,single treat(18)  trperiod(53) posttrend fig replace
itsa pm10,single treat(18)  trperiod(48) posttrend fig replace

// 以维度接近的城市作为对照组 28 14 17 1  成都 重庆 杭州 上海作为对照组
itsa pm25, treat(18) contid(14 28 17 1 ) trperiod(53) lag(1) replace fig posttrend
itsa pm25, treat(18) contid(14 28 17 1 ) trperiod(48) lag(1) replace fig posttrend

// 以临近省作为对照组 重庆 长沙 南昌 合肥 郑州 陕西
itsa pm25, treat(18) contid(28 31 8 25 27 ) trperiod(53) lag(1) replace fig posttrend

itsa pm25, treat(18) contid(28 31 8 25 27 ) trperiod(48) lag(1) replace fig posttrend