Hi,

You forgot to put the parenthesis in the way Petr told you : (-6.108576e-05)^(1/3) and the result is NaN. What do you want to preserve?


Alain




carol white wrote:
but with complex, I get complex numbers for the first and last elements:

(as.complex(tmp))^(1/3)
[1] 0.01969170+0.03410703i 0.03478442+0.00000000i 0.03285672+0.00000000i
[4] 0.08950802+0.00000000i 0.05848363+0.10129661i

whereas for the first element, we get the followings.

Moreover,
-6.108576e-05^(1/3)
 [1] -0.03938341
and
-(6.108576e-05^(1/3))
[1] -0.03938341
and
 -((6.108576e-05)^(1/3))
[1] -0.03938341

give the same results.

so using () doesn't preserve any thing

--- On Mon, 11/16/09, Petr PIKAL <[email protected]> wrote:

From: Petr PIKAL <[email protected]>
Subject: Odp: [R] ^ operator
To: "carol white" <[email protected]>
Cc: [email protected]
Date: Monday, November 16, 2009, 3:40 AM
Hi

AFAIK, this is issue of the preference of operators.
[email protected]
napsal dne 16.11.2009 11:24:59:

Hi,
I want to apply ^ operator to a vector but it is
applied to some of the
elements correctly and to some others, it generates
NaN. Why is it not able to
calculate -6.108576e-05^(1/3) even though it exists?


  tmp
[1] -6.108576e-05 4.208762e-05
3.547092e-05 7.171101e-04 -1.600269e-03
tmp^(1/3)
[1]        NaN 0.03478442
0.03285672 0.08950802        NaN

This computes (-a)^(1/3) which is not possible in real
numbers. You have to use as.complex(tmp)^(1/3) to get a result.

-6.108576e-05^(1/3)
[1] -0.03938341
this is actually
-(6.108576e-05^(1/3))

Regards
Petr


______________________________________________
[email protected]
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.






______________________________________________
[email protected] 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.


--
Alain Guillet
Statistician and Computer Scientist

SMCS - Institut de statistique - Université catholique de Louvain
Bureau c.316
Voie du Roman Pays, 20
B-1348 Louvain-la-Neuve
Belgium

tel: +32 10 47 30 50

______________________________________________
[email protected] 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.

Reply via email to