Re: [Rd] Wierd nlm behaviour in 2.10.1 and 2.12.0 [Sec=Unclassified]

2010-10-13 Thread Ravi Varadhan
easy to derive. Ravi. -Original Message- From: r-devel-boun...@r-project.org [mailto:r-devel-boun...@r-project.org] On Behalf Of Troy Robertson Sent: Wednesday, October 13, 2010 1:51 AM To: 'r-devel@R-project.org' Subject: [Rd] Wierd nlm behaviour in 2.10.1 and 2.12.0 [Sec=Unclassified

Re: [Rd] Wierd nlm behaviour in 2.10.1 and 2.12.0 [Sec=Unclassified]

2010-10-13 Thread Ravi Varadhan
'; 'r-devel@R-project.org' Subject: Re: [Rd] Wierd nlm behaviour in 2.10.1 and 2.12.0 [Sec=Unclassified] Hi Troy, I think that your problem is poorly scaled. You have variables that vary over several orders of magnitude. This means that any small changes (due to subtle differences between R versions

Re: [Rd] Wierd nlm behaviour in 2.10.1 and 2.12.0 [Sec=Unclassified]

2010-10-13 Thread Troy Robertson
Varadhan'; Troy Robertson; 'r-devel@R-project.org' Subject: RE: [Rd] Wierd nlm behaviour in 2.10.1 and 2.12.0 [Sec=Unclassified] Try this: Myfunc2 - function(rootM,Abund,Loss,OtherM) { # I have scaled the function (Loss/Abund - (rootM/(rootM+OtherM)* (1-exp(-(rootM+OtherM)^2 } nlm

Re: [Rd] Wierd nlm behaviour in 2.10.1 and 2.12.0 [Sec=Unclassified]

2010-10-13 Thread Ravi Varadhan
, October 13, 2010 9:08 pm Subject: RE: [Rd] Wierd nlm behaviour in 2.10.1 and 2.12.0 [Sec=Unclassified] To: 'Ravi Varadhan' rvarad...@jhmi.edu, 'r-devel@R-project.org' r-devel@r-project.org Thanks heaps Ravi, I will have to take the modified function to my modeller/statistician and get him

[Rd] Wierd nlm behaviour in 2.10.1 and 2.12.0 [Sec=Unclassified]

2010-10-12 Thread Troy Robertson
Hi all, When upgrading to 2.11.1 recently I noticed different results being produced by my code. After much digging I have finally narrowed it to a call to nlm(). This can be replicated by: FixedRemovals-1836180125888 AbStageInitial-2223033830403 Rates- 0.3102445

[Rd] Wierd nlm behaviour in 2.10.1 and 2.12.0 [Sec=Unclassified]

2010-10-12 Thread Troy Robertson
Hi all, When upgrading to 2.11.1 recently I noticed different results being produced by my code. After MUCH digging into my code I have finally narrowed it to a call to nlm(). The problem actually occurs in 2.10.1 and 2.12.0. This can be replicated by trying the code below in some different