chyihung.hsu@pharma.Novartis.com
2001-Mar-19 21:52 UTC
[Rd] A limitation for polyroot ? (PR#880)
Dear R Development Team, I have encountered the following difficulty in using the function polyroot under either NT4.0 (R version 1.2.1) or linux (R version 0.90.1). In the provided example, the non-zero root of c(0,0,0,1) depends on the results of the previous call of polyroot. R : Copyright 2001, The R Development Core Team Version 1.2.1 (2001-01-15) R is free software and comes with ABSOLUTELY NO WARRANTY. You are welcome to redistribute it under certain conditions. Type `license()' or `licence()' for distribution details. R is a collaborative project with many contributors. Type `contributors()' for more information. Type `demo()' for some demos, `help()' for on-line help, or `help.start()' for a HTML browser interface to help. Type `q()' to quit R. [Previously saved workspace restored]> polyroot(c(6,-5,1))[1] 2+1.615587e-27i 3-1.615587e-27i> polyroot(c(0,0,0,1))[1] 0+0.000000e+00i 0+0.000000e+00i 3-1.615587e-27i> polyroot(c(4,-4,1))[1] 2+1.110223e-16i 2-1.110223e-16i> polyroot(c(0,0,0,1))[1] 0+0.000000e+00i 0+0.000000e+00i 2-1.110223e-16i>Thank you for your attention. Best regards, Chyi-Hung Hsu -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-devel mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html Send "info", "help", or "[un]subscribe" (in the "body", not the subject !) To: r-devel-request@stat.math.ethz.ch _._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._
>Dear R Development Team,I think most of the core team is at a conference, but I seem to recall that this problem was fixed in 1.2.2. You should probably upgrade. Below is what I get. Paul Gilbert _____ R : Copyright 2001, The R Development Core Team Version 1.2.2 (2001-02-26) R is free software and comes with ABSOLUTELY NO WARRANTY. You are welcome to redistribute it under certain conditions. Type `license()' or `licence()' for distribution details. R is a collaborative project with many contributors. Type `contributors()' for more information. Type `demo()' for some demos, `help()' for on-line help, or `help.start()' for a HTML browser interface to help. Type `q()' to quit R.> polyroot(c(6,-5,1))[1] 2+1.615587e-27i 3-1.615587e-27i> polyroot(c(0,0,0,1))[1] 0+0i 0+0i 0+0i> polyroot(c(4,-4,1))[1] 2+1.110223e-16i 2-1.110223e-16i> polyroot(c(0,0,0,1))[1] 0+0i 0+0i 0+0i>-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-devel mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html Send "info", "help", or "[un]subscribe" (in the "body", not the subject !) To: r-devel-request@stat.math.ethz.ch _._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._
On Mon, 19 Mar 2001 chyihung.hsu@pharma.Novartis.com wrote:> Dear R Development Team, > > I have encountered the following difficulty in using the function polyroot > under either NT4.0 (R version 1.2.1) or linux (R version 0.90.1). In the > provided example, the non-zero root of c(0,0,0,1) depends on the results of > the previous call of polyroot.This is AFAIK already solved in 1.2.2. Please upgrade, at least that Linux version ....> > R : Copyright 2001, The R Development Core Team > Version 1.2.1 (2001-01-15) > > R is free software and comes with ABSOLUTELY NO WARRANTY. > You are welcome to redistribute it under certain conditions. > Type `license()' or `licence()' for distribution details. > > R is a collaborative project with many contributors. > Type `contributors()' for more information. > > Type `demo()' for some demos, `help()' for on-line help, or > `help.start()' for a HTML browser interface to help. > Type `q()' to quit R. > > [Previously saved workspace restored] > > > polyroot(c(6,-5,1)) > [1] 2+1.615587e-27i 3-1.615587e-27i > > polyroot(c(0,0,0,1)) > [1] 0+0.000000e+00i 0+0.000000e+00i 3-1.615587e-27i > > polyroot(c(4,-4,1)) > [1] 2+1.110223e-16i 2-1.110223e-16i > > polyroot(c(0,0,0,1)) > [1] 0+0.000000e+00i 0+0.000000e+00i 2-1.110223e-16i > > > > Thank you for your attention. > > Best regards, > > Chyi-Hung Hsu > > > -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- > r-devel mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html > Send "info", "help", or "[un]subscribe" > (in the "body", not the subject !) To: r-devel-request@stat.math.ethz.ch > _._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._ >-- Brian D. Ripley, ripley@stats.ox.ac.uk Professor of Applied Statistics, http://www.stats.ox.ac.uk/~ripley/ University of Oxford, Tel: +44 1865 272861 (self) 1 South Parks Road, +44 1865 272860 (secr) Oxford OX1 3TG, UK Fax: +44 1865 272595 -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-devel mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html Send "info", "help", or "[un]subscribe" (in the "body", not the subject !) To: r-devel-request@stat.math.ethz.ch _._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._