Next, I do not know if you understand this or not, but if you interept the answer (even though its running on cluster), you ought to be in a position to obtain the NR residuals for the final three (I usually use 3 but you can enter even more) iterations. https://a-solution.info/