Re: [R] selecting values that are unique, instead of selecting unique values

From: Prof Brian Ripley <ripley_at_stats.ox.ac.uk>
Date: Wed, 25 Jun 2008 17:59:35 +0100 (BST)

On Wed, 25 Jun 2008, Gabor Csardi wrote:

> Wow, that is smart, although is seems to be overkill.....
> I guess 'duplicated' is better than O(n^2), is it really?

Yes as it hashes, but the overhead on short vectors is high since it always hashes.

>
> Gabor
>
> On Wed, Jun 25, 2008 at 05:43:30PM +0100, Prof Brian Ripley wrote:
>> On Wed, 25 Jun 2008, Marc Schwartz wrote:
>>
>>> on 06/25/2008 11:19 AM Daren Tan wrote:
>>>>
>>>> unique(c(1:10,1)) gives 1:10 (i.e. unique values), is there any
>>>> method to get only 2:10 (i.e. values that are unique) ?
>>>>
>>>
>>> The easiest might be:
>>>
>>>> Vec
>>> [1] 1 2 3 4 5 6 7 8 9 10 1
>>>
>>>> Vec[table(Vec) == 1]
>>> [1] 2 3 4 5 6 7 8 9 10
>>
>> I don't think that is right: you are relying on recycling indices. Try
>>
>> Vec <- c(1,1:10)
>> Vec[table(Vec) == 1]
>>
>> which should be the same.
>>
>> I was about to write
>>
>> tab <- table(Vec)
>> names(tab)[tab==1]
>>
>> but that gives a character vector. Here's a different way:
>>
>> Vec[rowSums(outer(Vec, Vec, "=="))==1]

-- 
Brian D. Ripley,                  ripley_at_stats.ox.ac.uk
Professor of Applied Statistics,  http://www.stats.ox.ac.uk/~ripley/
University of Oxford,             Tel:  +44 1865 272861 (self)
1 South Parks Road,                     +44 1865 272866 (PA)
Oxford OX1 3TG, UK                Fax:  +44 1865 272595

______________________________________________
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 Wed 25 Jun 2008 - 17:12:54 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 Wed 25 Jun 2008 - 17:31:14 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