R/models__linear__etc__X.x.delta.xi.R

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

#' Preform X multiply Dev xi
#'
#' <details>
#' @param X NA
#' @param delta.knots  NA
#' @param q.i  NA
#' @return NA
#' @references NA
#' @author Feng Li, Department of Statistics, Stockholm University, Sweden.
#' @export
X.x.delta.xi <- function(X, delta.knots, q.i)
{
  knots.comp <- names(delta.knots)
  knots.comp.len <- length(knots.comp)
  q <- sum(q.i)

  out <- list()
  for(i in 1:knots.comp.len)
    {
      out.mat <- crossprod(X, delta.knots[[i]])
      out.ary <- array(out.mat, c(q, prod(dim(out.mat))/(q*q.i[i+1]) ,q.i[i+1]))
      out.lst <- array2list(out.ary, 3)
      out[[i]] <- out.lst
    }

  names(out) <- knots.comp

  return(out)
}

##----------------------------------------------------------------------------------------
## TESTS: PASSED
##----------------------------------------------------------------------------------------

## system.time(X.x.delta.xi(X.mats, delta.knots) )