newton method

9 messages
Open this post in threaded view
|

newton method

 Hi R-users, Does R has a topic on newton's method? Thank you for the info. ______________________________________________ [hidden email] mailing list https://stat.ethz.ch/mailman/listinfo/r-helpPLEASE do read the posting guide http://www.R-project.org/posting-guide.htmland provide commented, minimal, self-contained, reproducible code.
Open this post in threaded view
|

Re: newton method

 Take a look at the functionsnlm(), optim() in the stats package and   maxNR() in the maxLik package. On Mar 22, 2009, at 11:15 PM, Roslina Zakaria wrote: > Does R has a topic on newton's method? _____________________________ Professor Michael Kubovy University of Virginia Department of Psychology Postal Address:         P.O.Box 400400, Charlottesville, VA 22904-4400 Express Parcels Address:         Gilmer Hall, Room 102, McCormick Road, Charlottesville, VA 22903 Office:    B011; Phone: +1-434-982-4729 Lab:        B019; Phone: +1-434-982-4751 WWW:    http://www.people.virginia.edu/~mk9y/Skype name: polyurinsane         [[alternative HTML version deleted]] ______________________________________________ [hidden email] mailing list https://stat.ethz.ch/mailman/listinfo/r-helpPLEASE do read the posting guide http://www.R-project.org/posting-guide.htmland provide commented, minimal, self-contained, reproducible code.
Open this post in threaded view
|

Re: newton method

Open this post in threaded view
|

Re: newton method

 In reply to this post by Roslina Zakaria I'm not sure what you meant by "a topic on newton's method" (algorithm? demo?), but the demonstration in the package 'animation' might help: install.packages('animation') par(pch = 20) ani.options(nmax = 50) newton.method(function(x) 5 * x^3 - 7 * x^2 - 40 *     x + 100, 7.15, c(-6.2, 7.1)) Regards, Yihui -- Yihui Xie <[hidden email]> Phone: +86-(0)10-82509086 Fax: +86-(0)10-82509086 Mobile: +86-15810805877 Homepage: http://www.yihui.nameSchool of Statistics, Room 1037, Mingde Main Building, Renmin University of China, Beijing, 100872, China On Mon, Mar 23, 2009 at 11:15 AM, Roslina Zakaria <[hidden email]> wrote: > > Hi R-users, > > Does R has a topic on newton's method? > > Thank you for the info. > ______________________________________________ [hidden email] mailing list https://stat.ethz.ch/mailman/listinfo/r-helpPLEASE do read the posting guide http://www.R-project.org/posting-guide.htmland provide commented, minimal, self-contained, reproducible code.
Open this post in threaded view
|

Newton method again

 In reply to this post by Roslina Zakaria Hi Rolf, I would like to extend the problem that I asked you before regarding the newton method using 4 functions with 4 parameters.  My functions involve the modified bessel function of the first kind which I can type them without any problem. The big problem is the Jacobian matrix. I use Maple 11.0 to get the algebraic expression for the Jacobian matrix.  The problem is that the Jacobian matrix is so complicated and it is quite impossible to type the expression as I end up of more than 50 pages for the algebraic form of the Jacobian.  Can you suggest any method to solve the problems? Thank you so much for your attention and help. Roslina. ______________________________________________ [hidden email] mailing list https://stat.ethz.ch/mailman/listinfo/r-helpPLEASE do read the posting guide http://www.R-project.org/posting-guide.htmland provide commented, minimal, self-contained, reproducible code.
Open this post in threaded view
|

Re: Newton method again

 Hi, I do not know what your problem is, but it seems like you want to solve a system of non-linear equations. Look at the function dfsane() in the package "BB".  It doesn't require you to specify jacobians. require(BB) ?dfsane Ravi. ____________________________________________________________________ Ravi Varadhan, Ph.D. Assistant Professor, Division of Geriatric Medicine and Gerontology School of Medicine Johns Hopkins University Ph. (410) 502-2619 email: [hidden email] ----- Original Message ----- From: Roslina Zakaria <[hidden email]> Date: Tuesday, April 7, 2009 10:35 pm Subject: [R] Newton method again To: Rolf Turner <[hidden email]>, R help forum <[hidden email]> >  Hi Rolf, >   >  I would like to extend the problem that I asked you before regarding > the newton method using 4 functions with 4 parameters.  My functions > involve the modified bessel function of the first kind which I can > type them without any problem. >  The big problem is the Jacobian matrix. >  I use Maple 11.0 to get the algebraic expression for the Jacobian > matrix.  The problem is that the Jacobian matrix is so complicated and > it is quite impossible to type the expression as I end up of more than > 50 pages for the algebraic form of the Jacobian.  Can you suggest any > method to solve the problems? >   >  Thank you so much for your attention and help. >   >   >   >  Roslina. >   >   >   >   >  ______________________________________________ >  [hidden email] mailing list >   >  PLEASE do read the posting guide >  and provide commented, minimal, self-contained, reproducible code. ______________________________________________ [hidden email] mailing list https://stat.ethz.ch/mailman/listinfo/r-helpPLEASE do read the posting guide http://www.R-project.org/posting-guide.htmland provide commented, minimal, self-contained, reproducible code.
Open this post in threaded view
|