*数据处理与清洗

egen m10 = group(m10)
egen m9 = group(m9)
xtset id year

负二项回归   跑不出，用tobit
xtnbreg y x i.m9 i.m10
est store model1

*随机效应的面板Tobit回归

gen x11=x1*x1

*y1与x1正相关；

xttobit y xa i.m9 i.m10,tobit
est store model1
esttab  model1  using y1d.rtf,stats(N r2_a) star(* 0.1 ** 0.05 *** 0.01) se ar2