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'
}