做二值因变量回归后，有一段程序，请帮忙把各个步骤解释一下。还有最终得到的是什么？
probit y x1 x2  x3 x4(注：这里x1 x2是连续变量，  x3 x4是虚拟变量)

predict xb if e(sample),xb
predict p if e(sample),p

local cont "x1,x2"
local dummies "x3 x4"

gen t_var=normalden(xb)

foreach c of local cont{
       gen pe_'c'=t_var*_b['c']
}
foreach d of local dummies{
       gen pe_'d'=p-norm(xb-_b['d'])
       replace pe_'d'=norm(xb+_b['d'])-p if 'd'==0
}
summ ‘cont’ 'dummies' ,detail     [注：最后这一步回归出现错误提示，不知道‘cont’ 'dummies' 代表什么意思]