我想像截图那样，求广义残差和进行W检验，请问命令是否正确？

总共才几步，真心不多，请认真看！附件为截图所用论文原文、自己写的do文档、截图（和帖子内容一样，方便下载）。

1. *第一个方程的广义残差
   
2. probit insur p1 p2  
   
3. gen p1b=_b[p1]*p1
   
4. gen p2b=_b[p2]*p2
   
5. gen XB=p1b+p2b
   
6. gen r=(normalden(XB)/normal(XB))*insur - (1-insur)*(normalden(XB)/normal(-XB))
   
7. 
   
8. *同样，求第二个方程的广义残差：c
   
9. 
   
10. *计算W统计量，看显著性
    
11. gen rc=r*c
    
12. egen rc_he=total(rc)
    
13. gen rc_he2=rc_he^2
    
14. gen r2=r^2
    
15. gen c2=c^2
    
16. gen r2c2=r2*c2
    
17. egen r2c2_he=total(r2c2)
    
18. gen W= rc_he2/r2c2_he
    
19. 
    
20. *计算W对应P分位数,看看是否显著
    
21. disp chi2tail(1,W)

复制代码