Without getting into the gory details, I was having problems getting the results from the scipy minimize routine to match output from Nonlinear and Life Distribution. I update both to JMP 19 (running jpip to installed needed packages) and my OS to Tahoe. Now when I run the same Python script in JMP 18 or JMP 19 everything matches. Any ideas what might be going on? I'm attaching the data set and code if anyone who hasn't updated to Tahoe is willing to give it a try. The values you should get are approx. 2.2325 and 13716.73. I was getting values of about 3.01 and 4700 before the updates.