# Re: [R] Odp: strange issue with "which" on "seq"

From: Joris Meys <jorismeys_at_gmail.com>
Date: Wed, 09 Jun 2010 14:31:09 +0200

use ?round is the easier way
v <- round(seq(-0.5,0.5,by=0.1),1)

On Wed, Jun 9, 2010 at 2:20 PM, Nikhil Kaza <nikhil.list_at_gmail.com> wrote:
> which(abs(v - .1) <= .Machine\$double.eps)
>
> seems to me too cumbersome to write. Any other easier way? all.equal does
> not quite work
>
> Nikhil
>
> On Jun 9, 2010, at 7:54 AM, Petr PIKAL wrote:
>
>> Hi
>>
>> r-help-bounces_at_r-project.org napsal dne 09.06.2010 13:16:40:
>>
>>> Dear R community,
>>> I am puzzled by the following lines:
>>>
>>>> v <- seq(-0.5,0.5,by=0.1)
>>>> v
>>>
>>>  [1] -0.5 -0.4 -0.3 -0.2 -0.1  0.0  0.1  0.2  0.3  0.4  0.5
>>>>
>>>> which(v == -0.4)
>>>
>>>  [1] 2
>>>>
>>>> which(v == 0)
>>>
>>>  [1] 6
>>>>
>>>> which(v == 0.1)
>>>
>>> integer(0)
>>>>
>>>> which(v == 0.2)
>>>
>>> integer(0)
>>>>
>>>> which(v == 0.3)
>>>
>>> integer(0)
>>>>
>>>> which(v == 0.4)
>>>
>>>  [1] 10
>>>>
>>>> which(v == 0.5)
>>>
>>>  [1] 11
>>>
>>> Why "which" can only match some of the values in "v"? Are the numbers
>>> generated by "seq" not exact fractional numbers?
>>> Please, help me to understand this.
>>
>> Well FAQ 7.31 was not here some time. Computing in binary results in
>> finite precision of fractional numbers.
>>
>> v <- seq(-0.5,0.5,by=0.1)
>> v[7]-0.1
>> [1] 8.326673e-17
>>
>> Regards
>> Petr
>>
>>
>>
>>>
>>> J
>>>
>>> Membrane Protein Laboratory (MPL)
>>> Diamond Light Source Ltd
>>> Diamond House
>>> Harewell Science and Innovation Campus
>>> Chilton, Didcot
>>> Oxfordshire OX11 0DE
>>>
>>>
>>> --
>>> This e-mail and any attachments may contain confidential...{{dropped:8}}
>>>
>>> ______________________________________________
>>> R-help_at_r-project.org mailing list
>>> https://stat.ethz.ch/mailman/listinfo/r-help
>>
>> http://www.R-project.org/posting-guide.html
>>>
>>> and provide commented, minimal, self-contained, reproducible code.
>>
>> ______________________________________________
>> R-help_at_r-project.org mailing list
>> https://stat.ethz.ch/mailman/listinfo/r-help
>> http://www.R-project.org/posting-guide.html
>> and provide commented, minimal, self-contained, reproducible code.
>
> ______________________________________________
> R-help_at_r-project.org mailing list
> https://stat.ethz.ch/mailman/listinfo/r-help
> and provide commented, minimal, self-contained, reproducible code.
>

```--
Joris Meys
Statistical consultant

Ghent University
Faculty of Bioscience Engineering
Department of Applied mathematics, biometrics and process control

tel : +32 9 264 59 87
Joris.Meys_at_Ugent.be
-------------------------------
Disclaimer : http://helpdesk.ugent.be/e-maildisclaimer.php

______________________________________________
R-help_at_r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help