我有一个大约有500,000个观测值的数据集,对于每个观测,我试图求解一个由两个非线性方程组成的系统。我这样做如下:
proc model data=myData;
eq.equationA= some non-linear equation in x and y;
eq.equationB=some different non-linear equation in x and y;
solve x y;
run;
quit;
从日志中,我得到的确切错误是: 关于要求解的变量的方程的偏导数的牛顿方法雅可比矩阵在观察1275处是单数的,用于迭代2.方程组不能被求解。
这里的问题是当这种情况发生时,SAS会立即终止proc模型。这是令人惊讶的,因为通常SAS在遇到问题时继续运行。如何强制SAS简单地跳过不良观察(1275)并继续运行模型并进行后续观察? (我不能简单地删除观察1275,因为我知道其他观察可能会导致同样的错误)