*若要各估计系数与定义式各系数的符号一致,可试以下过程(各系数均显著):
clear
set ob 1000
set se 10000
g x1=rnormal()
g x2=rnormal()
g x3=rnormal()
g w=rnormal()<0
g y=(0.8+x1+2*x2+3*x3+4*w+5*c.x1#w>=-rnormal())
*或
g z=(0.8+x1+2*x2+3*x3+4*w+5*c.x1#w>=rnormal())
prob y x* w c.x1#w
*或
prob z x* w c.x1#w


雷达卡



京公网安备 11010802022788号







