Re: [R] problem with Gauss Hermite ( x and w )

2012-05-10 Thread R. Michael Weylandt
Are you getting caught on order of operations? Note that unary minus has lower precedence than exponentiation (as it does in math) so -2.5^(-2.4) is x <- 2.5^(-2.4) -x Otherwise, I'm not at all sure what your question is: can you give an example of what you think you should get (and how to get

Re: [R] problem with Gauss Hermite ( x and w )

2012-05-10 Thread casperyc
Rui Barradas wrote > > > real^real is not necessarily real. > > The most well known example is (-1)^0.5 = imaginary unit. > Damn, can't believe it! It's a silly mistake! Now that something wonders me is that when applying the Gaussian Hermit, sum w f(x_i) What happens when f(x_i) does not

Re: [R] problem with Gauss Hermite ( x and w )

2012-05-10 Thread Rui Barradas
Hello, > Negative powers mean they take the reciprocal and as far as I am > concerned, real^real is just a real number. > Am I mistaking something basic? Yes, you are. real^real is not necessarily real. The most well known example is (-1)^0.5 = imaginary unit. When you say that -2.5^(-2.4) is r

Re: [R] problem with Gauss Hermite ( x and w )

2012-05-10 Thread casperyc
Hi, I know what complex number are, but I am not sure what you meant by that? ##CODES### > 2.5^(-2.4) [1] 0.1109032 > -2.5^(-2.4) [1] -0.1109032 ##CODES### works fine. Negative powers mean they take the reciprocal and as far as I am concerned, real^real

Re: [R] problem with Gauss Hermite ( x and w )

2012-05-09 Thread R. Michael Weylandt
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 wrote: > Hi all, > > I am using the 'gaussHermite' function from

[R] problem with Gauss Hermite ( x and w )

2012-05-09 Thread casperyc
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 ##