I am not familiar with your simulator engine - therefore, I cannot explicitely list the relevant steps.
But a general question: What do you expect from this analysis ?
My recommendation is instead: Take one characteristic curve ID=f(VDS) and try to expand the slope of the curve to the left (asymptote to negative VDS values) until this line crosses the horizontal axis. Then you have the Early voltage VA which with a good approximation is independent on the VGS voltage because all these lines meet in the vicinity of VA.