Cool.
Is one of the results the definitive function to use to mimic Octave's poly
function?
There were about 10 returns, and none clearly identified themselves as the
definitive functionality representing the poly capability, so any insight and
clarification is appreciated.
----- Original Message ----
From: David Winsemius <dwinsem...@comcast.net>
To: Jason Rupert <jasonkrup...@yahoo.com>
Cc: R-help@r-project.org
Sent: Sat, November 28, 2009 9:23:23 AM
Subject: Re: [R] R function that duplicates Octave's poly function?
On Nov 28, 2009, at 9:33 AM, Jason Rupert wrote:
>
> By any chance is anyone aware of an R function that duplicates
> Octave's poly function?
>
> Here is a description of Octave's poly function:
>
> Function File: poly (A)
> If A is a square N-by-N matrix, `poly (A)' is the row vector of
> the coefficients of `det (z * eye (N) - a)', the characteristic
> polynomial of A. As an example we can use this to find the
> eigenvalues of A as the roots of `poly (A)'.
> roots(poly(eye(3)))
> => 1.00000 + 0.00000i
> => 1.00000 - 0.00000i
> => 1.00000 + 0.00000i
> In real-life examples you should, however, use the `eig'
> function for computing eigenvalues.
>
> If X is a vector, `poly (X)' is a vector of coefficients of the
> polynomial whose roots are the elements of X. That is, of C is a
> polynomial, then the elements of `D = roots (poly (C))' are
> contained in C. The vectors C andD are, however, not equal due
> to sorting and numerical errors.
>
> Thanks again for any insights and feedback.
RSiteSearch("characteristic polynomial")
>
> ______________________________________________
> R-help@r-project.org mailing list
> https://stat.ethz.ch/mailman/listinfo/r-help
> PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
> and provide commented, minimal, self-contained, reproducible code.
David Winsemius, MD
Heritage Laboratories
West Hartford, CT
______________________________________________
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.
