Try this. apseq() sorts the input and appends a
sequence number: 0, 1, ... to successive
occurrences of each value. Apply that to both
vectors transforms it into a problem that works
with ordinary match:
> lookupTable <- c("a", "a","b","c","d","e","f")
> matchSample <- c("a", "a","a","b","d")
>
> # sort and append sequence no
> apseq <- function(x) {
+ x <- sort(x)
+ s <- cumsum(!duplicated(x))
+ paste(x, seq(s) - match(s, s))
+ }
>
> match(apseq(matchSample), apseq(lookupTable))
[1] 1 2 NA 3 5
On Sun, Jun 22, 2008 at 10:57 PM, <[EMAIL PROTECTED]> wrote:
> Hi folks,
>
> Can anyone suggest an efficient way to do "matching without
> replacement", or "one-to-one matching"? pmatch() doesn't quite provide
> what I need...
>
> For example,
>
> lookupTable <- c("a","b","c","d","e","f")
> matchSample <- c("a","a","b","d")
> ##Normal match() behaviour:
> match(matchSample,lookupTable)
> [1] 1 1 2 4
>
> My problem here is that both "a"s in matchSample are matched to the same
> "a" in the lookup table. I need the elements of the lookup table to be
> excluded from the table as they are matched, so that no match can be
> found for the second "a".
>
> Function pmatch() comes close to what I need:
>
> pmatch(matchSample,lookupTable)
> [1] 1 NA 2 4
>
> Yep! However, pmatch() incorporates partial matching, which I
> definitely don't want:
>
> lookupTable <- c("a","b","c","d","e","aaaaaaaaf")
> matchSample <- c("a","a","b","d")
> pmatch(matchSample,lookupTable)
> [1] 1 6 2 4
> ## i.e. the second "a", matches "aaaaaaaaf" - I don't want this.
>
> Of course, when identical items ARE duplicated in both sample and lookup
> table, I need the matching to reflect this:
>
> lookupTable <- c("a","a","c","d","e","f")
> matchSample <- c("a","a","c","d")
> ##Normal match() behaviour
> match(matchSample,lookupTable)
> [1] 1 1 3 4
>
> No good - pmatch() is better:
>
> lookupTable <- c("a","a","c","d","e","f")
> matchSample <- c("a","a","c","d")
> pmatch(matchSample,lookupTable)
> [1] 1 2 3 4
>
> ...but we still have the partial matching issue...
>
> ##And of course, as per the usual behaviour of match(), sample elements
> missing from the lookup table should return NA:
>
> matchSample <- c("a","frog","e","d") ; print(matchSample)
> match(matchSample,lookupTable)
>
> Is there a nifty way to get what I'm after without resorting to a for
> loop? (my code's already got too blasted many of those...)
>
> Thanks,
>
> Alec Zwart
> CMIS CSIRO
> [EMAIL PROTECTED]
>
> ______________________________________________
> R-help@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.
>
______________________________________________
R-help@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.
