R/decreasing.R

In ORIClust: Order-Restricted Information Criterion-Based Clustering Algorithm

decreasing <-
function(data,x,n.rep){ 
    k <- length(n.rep)
    data <- as.numeric(data)
    n <- sum(n.rep)
    s <- x
    w <- n.rep
    a <- s*w
    is <-  numeric(k)
    uu <- 0
    while(uu < k) {
        v <- uu + 1
        b <- cumsum(a[v:k] )
        d <- cumsum(w[v:k])
        b <- b/d
        u <- max(b)
        m <- match(u, rev(b)) 
        m <- length(b)-m+v
        is[v:m] <- u
        uu <- m
    }
    sig <- 0;
    len1 <- 1;
    len2 <- 0;
    for(i in 1:k){
        muu <- is[i]
        len2 <- len2+n.rep[i]
        x1 <- data[len1:len2]
        len1 <- len1+n.rep[i]
        sig <- sig+sum((x1-muu)^2)
    }
    sig <- sig/n;
    like <- -0.5*n*log(2*pi)-0.5*n*log(sig)-0.5*n
    list( logelr=like, mu=is)
}