This might work:

> tmp <- c(24,25,29,35,36,37,38,39,40,41,42,43,44,45,46,47,68,69,70,71)
> # generate breaks
> group <- c(0, cumsum(diff(tmp) != 1))
> tapply(tmp, group, summary)
$`0`
   Min. 1st Qu.  Median    Mean 3rd Qu.    Max.
  24.00   24.25   24.50   24.50   24.75   25.00

$`1`
   Min. 1st Qu.  Median    Mean 3rd Qu.    Max.
     29      29      29      29      29      29

$`2`
   Min. 1st Qu.  Median    Mean 3rd Qu.    Max.
     35      38      41      41      44      47

$`3`
   Min. 1st Qu.  Median    Mean 3rd Qu.    Max.
  68.00   68.75   69.50   69.50   70.25   71.00

>
>



On 7/24/06, Kevin J Emerson <[EMAIL PROTECTED]> wrote:
>
> Hello R-helpers!
>
> I have a question concerning extracting sequence information from a
> vector.  I have a vector (representing the bins of a time series where
> the frequency of occurrences is greater than some threshold) where I
> would like to extract the min, median and max of each group of
> consecutive numbers.
>
> For Example:
>
> tmp <- c(24,25,29,35,36,37,38,39,40,41,42,43,44,45,46,47,68,69,70,71)
>
> I would like to have the max,min,median of the following groups:
>
> 24,25
> 29
> 35,36,37,38,39,40,41,42,43,44,45,46,47,
> 68,69,70,71
>
> I would like to be able to perform this for many time series so an
> automated process would be nice.  I am hoping to use this as a peak
> detection protocol.
>
> Any advice would be greatly appreciated,
> Kevin
>
> -----
> -----
> Kevin J Emerson
> Center for Ecology and Evolutionary Biology
> 1210 University of Oregon
> Eugene, OR 97403
> USA
> [EMAIL PROTECTED]
>
> ______________________________________________
> R-help@stat.math.ethz.ch 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.
>



-- 
Jim Holtman
Cincinnati, OH
+1 513 646 9390

What is the problem you are trying to solve?

        [[alternative HTML version deleted]]

______________________________________________
R-help@stat.math.ethz.ch 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.

Reply via email to