R/idx.beta2b.R

Defines functions idx.beta2b

Documented in idx.beta2b

##' Convert indices from beta to b
##'
##' @export
idx.beta2b <- function(p, q, q.i)
  {
    idx4beta <- matrix(1:(p*q), q)
    cum.q.i <- c(0, cumsum(q.i))

    out4b <- matrix(0, 1:(p*q))
    for(i in 1:length(q.i))
      {
        idx4bi <- (1 + cum.q.i[i]):cum.q.i[i+1]
        out4b[(1 + cum.q.i[i]*p):(cum.q.i[i+1]*p)] <- idx4beta[idx4bi, ]
      }
    return(out4b)
  }
thiyangt/fformpp documentation built on Jan. 5, 2024, 5:44 a.m.