Logit回归后计算P值

我用logit做100组的回归，然后需要提取出系数与系数对应的p值，
我用命令是statsby
stataby  _b _se r2=e(r2_p) n=e(df_r), clear by(county): logit y  x
后面该怎么计算？

gen b=.
gen bp=.
forvalues i=1/100 {
logit y x if country==`i'
matrix a1=e(b)
matrix  a2=e(V)
local b1=a1[1,1]
local b2=a2[1,1]
local b3=`b1'/`b2'^0.5
local b4=e(rank)
local b5=(1-normal(abs(`b3')))*2
replace b=`b1' in `i'
replace bp=`b5' in `i'
}