gen population1=.
summarize population
local a1=r(mean)
local a2=r(se)
replace population=1 if population<`a1'-2*`a2' &  population>`a1'-*`a2'
....