[R] uniq -c
William Dunlap
wdunlap at tibco.com
Tue Oct 16 22:46:54 CEST 2012
You said you wanted the equivalent of the Unix 'uniq -c' but said
that xtab's results were roughly right and the rle might be what
you want. rle() is the equivalent of 'uniq -c', they both output the
lengths of runs of identical elements. if the data is sorted they
are equivalent to using table() or xtabs().
Since you have sorted data try the following
isFirstInRun <- function(x) UseMethod("isFirstInRun")
isFirstInRun.default <- function(x) c(TRUE, x[-1] != x[-length(x)])
isFirstInRun.data.frame <- function(x) {
stopifnot(ncol(x)>0)
retval <- isFirstInRun(x[[1]])
for(column in x) {
retval <- retval | isFirstInRun(column)
}
retval
}
i <- which(isFirstInRun(yourDataFrame))
Then I think
data.frame(Count=diff(c(i, 1L+nrow(yourDataFrame))), yourDataFrame[i,])
gives you what you want. E.g.,
> yourDataFrame <- data.frame(x1=c(1,1,2,2,1), x2=c(11,11,11,12,11))
> i <- which(isFirstInRun(yourDataFrame))
> i
[1] 1 3 4 5
> data.frame(Count=diff(c(i, 1L+nrow(yourDataFrame))), yourDataFrame[i,])
Count x1 x2
1 2 1 11
3 1 2 11
4 1 2 12
5 1 1 11
It should be pretty quick. If you have missing values in your data frame,
you will have to make some decisions about whether they should be
considered equal to each other or not and modify isFirstInRun.default
accordingly.
Bill Dunlap
Spotfire, TIBCO Software
wdunlap tibco.com
> -----Original Message-----
> From: r-help-bounces at r-project.org [mailto:r-help-bounces at r-project.org] On Behalf
> Of Sam Steingold
> Sent: Tuesday, October 16, 2012 10:46 AM
> To: r-help at r-project.org; Duncan Murdoch
> Subject: Re: [R] uniq -c
>
> > * Duncan Murdoch <zheqbpu.qhapna at tznvy.pbz> [2012-10-16 12:47:36 -0400]:
> >
> > On 16/10/2012 12:29 PM, Sam Steingold wrote:
> >> x is sorted.
> > sparseby(data=x, INDICES=x, FUN=nrow)
>
> this takes forever; apparently, it does not use the fact that x is
> sorted (even then - it should not take more than a few minutes)...
>
> --
> Sam Steingold (http://sds.podval.org/) on Ubuntu 12.04 (precise) X 11.0.11103000
> http://www.childpsy.net/ http://openvotingconsortium.org http://www.memritv.org
> http://think-israel.org http://pmw.org.il http://thereligionofpeace.com
> Save the whales, feed the hungry, free the mallocs.
>
> ______________________________________________
> 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.
More information about the R-help
mailing list