R/isodecre.R

Defines functions isodecre

Documented in isodecre

isodecre <-
function(a, w){
   n <- length(a)
   is <- numeric(n)
   uu <- 0
   while(uu < n) {
      v <- uu + 1
      b <- cumsum(a[v:n] )
      d <- cumsum(w[v:n])
      b <- b/d
      u <- max(b)
      m <- match(u, rev(b)) 
      m <- length(b)-m+v
      is[v:m] <- u
      uu <- m
   }
   is
}

Try the ORIClust package in your browser

Any scripts or data that you put into this service are public.

ORIClust documentation built on June 23, 2022, 9:10 a.m.