Taking negative numbers to fractional powers gives NaNs.... that's just how it works.
Unless you want to use complex numbers (which R does not by default): as.complex(cc$x) ^ (2.5) Michael On Wed, May 9, 2012 at 7:22 PM, casperyc <caspe...@hotmail.co.uk> wrote: > Hi all, > > I am using the 'gaussHermite' function from the 'pracma' library > > ############ CODES ########### > library(pracma) > cc=gaussHermite(10) > cc$x^2 > cc$x^5 > cc$x^4 > ############ CODES ########### > > as far so good. However, it does NOT work for any NON integer values, say > > ############ CODES ########### > cc$x^(2.5) > cc$x^(-2.5) > ############ CODES ########### > > But just think about it numberically, it should work..... > > why this is the case? > > Is there a reason for getting the "NaN"s? > > Thanks! > > > > > > > > ----- > ###################### > PhD candidate in Statistics > Big R Fan > Big LEGO Fan > Big sTaTs Fan > ###################### > > -- > View this message in context: > http://r.789695.n4.nabble.com/problem-with-Gauss-Hermite-x-and-w-tp4622115.html > Sent from the R help mailing list archive at Nabble.com. > > ______________________________________________ > 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. ______________________________________________ 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.
