Re: [R] why doesn't ifelse work ?

From: Joshua Wiley <jwiley.psych_at_gmail.com>
Date: Thu, 28 Apr 2011 22:56:33 -0700

Hi Eric,

On Thu, Apr 28, 2011 at 6:46 PM, eric <ericstrom_at_aol.com> wrote:
> I have the following lines of code:
>
> ind <- rollapply(GSPC, 200, mean)
> signal <- ifelse(diff(ind, 5) > 0 , 1 , -1)

If you had looked at signal here, you would see that it is logical.

> signal[is.na(signal)] <- 0

but this is not doing what you almost certainly expect it is. Consider:

ind[is.na(ind)]
ind[!is.na(ind)]

nothing is selected. Gabor or Achim would know for sure, but I believe this is because zoos have methods for extraction (i.e., `[.zoo`), which does not work with matrices, and is.na(ind) returns a matrix. So all you are doing with:

signal[is.na(signal)] <- 0

is converting the class from logical to numeric

which is why you get 0s and 1s instead of TRUEs and FALSEs.

>
> I never get a value of -1 for signal even though I know diff(ind , 5) is
> less than zero frequently. It looks like when diff(ind , 5) is less than
> zero, signal gets set to 0 instead of - 1. Any ideas why ?  Here's some
> information on ind and diff(ind, 5) :
>
>> mode(diff(ind, 5) >0)
> [1] "logical"
>> class(diff(ind, 5) >0 )
> [1] "zoo"

This should have set off red flags. You are using a core, nongeneric function, ifelse, with a very special class of object. In fact, even your test is not a regular old logical object, it is a zoo. You can get lucky and everything still work as expected, and often times authors have written methods that accompany their classes so everything is handled automagically in the background without you ever needing to know about it. But in this case, ifelse is not generic, so there are not any special methods written for it (side note, there is an ifelse.zoo, but it also does not quite do what I think you want).

If you look at the code for ifelse, you will see that it returns the results of the test, AFTER updating any values that meet (or do not meet) the test AND are not NA. But wait, we already saw earlier what the results of: ind[!is.na(ind)] were, and you can tell from the class of your test (class(diff(ind, 5) > 0) that it is a zoo itself. This all seems to end up with ifelse() not seeing any values that match the criteria, so it simple returns the results of the the test itself (i.e., diff(ind, 5) > 0). Then you try to replace any NA values with 0, which converts the logical test results to numeric and because of how logical data is stored (FALSE = 0, TRUE = 1), you get a bunch of 0s and 1s.

David already alluded to this solution, but you can use coredata (see ?zoo for details) to get what I suspect are the results you are after:

signal <- ifelse(diff(coredata(ind), 5) > 0, 1, -1)

mind you "signal" will not be a zoo anymore.

Cheers,

Josh

>> str(diff(ind, 5) > 0 )
> ‘zoo’ series from 1990-05-31 to 2010-12-02
>  Data: logi [1:5171, 1] FALSE FALSE FALSE FALSE FALSE FALSE ...
>  - attr(*, "dimnames")=List of 2
>  ..$ : NULL
>  ..$ : chr "GSPC.Adjusted"
>  Index:  Date[1:5171], format: "1990-05-31" "1990-06-01" "1990-06-04"
> "1990-06-05" "1990-06-06" "1990-06-07" "1990-06-08" "1990-06-11" ...
>> class(ind)
> [1] "zoo"
>> mode(ind)
> [1] "numeric"
>> str(ind)
> ‘zoo’ series from 1990-05-23 to 2010-12-02
>  Data: num [1:5176, 1] 339 339 338 338 338 ...
>  - attr(*, "dimnames")=List of 2
>  ..$ : NULL
>  ..$ : chr "GSPC.Adjusted"
>  Index:  Date[1:5176], format: "1990-05-23" "1990-05-24" "1990-05-25"
> "1990-05-29" "1990-05-30" "1990-05-31" "1990-06-01" "1990-06-04"
>
> --
> View this message in context: http://r.789695.n4.nabble.com/why-doesn-t-ifelse-work-tp3482680p3482680.html
> Sent from the R help mailing list archive at Nabble.com.
>
> ______________________________________________
> R-help_at_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.
>

-- 
Joshua Wiley
Ph.D. Student, Health Psychology
University of California, Los Angeles
http://www.joshuawiley.com/

______________________________________________
R-help_at_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.
Received on Fri 29 Apr 2011 - 06:10:11 GMT

Archive maintained by Robert King, hosted by the discipline of statistics at the University of Newcastle, Australia.
Archive generated by hypermail 2.2.0, at Fri 29 Apr 2011 - 06:20:33 GMT.

Mailing list information is available at https://stat.ethz.ch/mailman/listinfo/r-help. Please read the posting guide before posting to the list.

list of date sections of archive