残差分析代码