R/models__linear__etc__idx.b2beta.R

In feng-li/movingknots: Efficient Bayesian multivariate surface regression

#' Make indices from b to beta.
#'
#' Details from the paper
#' @param p  NA
#' @param q NA
#' @param q.i NA
#' @return NA
#' @author Feng Li, Department of Statistics, Stockholm University, Sweden.
#' @export
idx.b2beta <- function(p, q, q.i)
  {
    idx4b <- 1:(p*q) ## The original indices for b
    cumidx <- c(0, cumsum(q.i))

    matidx4bi <- matrix(0, q, p)
    for(i in 1:length(q.i))
      {
        idx4bi <- (1+cumidx[i]*p):(cumidx[i+1]*p)

        matidx4bi[(1+cumidx[i]):(cumidx[i+1]), ] <- matrix(idx4b[idx4bi], q.i[i])
      }
    idx4beta <- as.vector(matidx4bi)
    return(idx4beta)
  }