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: [email protected] [mailto:[email protected]] On
> Behalf
> Of Sam Steingold
> Sent: Tuesday, October 16, 2012 10:46 AM
> To: [email protected]; Duncan Murdoch
> Subject: Re: [R] uniq -c
>
> > * Duncan Murdoch <[email protected]> [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.
>
> ______________________________________________
> [email protected] 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.
______________________________________________
[email protected] 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.
