R/02_algebrahelpers_vectorize.R

In SachaEpskamp/psychonetrics: Structural Equation Modeling and Confirmatory Network Analysis

# some helper functions:
Vec <- function(x){
  if (is.null(dim(x))){
    return(matrix(x,1,1))
  }
  
  if (is(x,"matrix")){
    return(matrix(c(x)))
  } else {
    # Assume Matrix package
    if (is(x,"sparseMatrix")){
      return(c(as.matrix(x))) # FIXME: This should be better...
      # if (length(x@x) == 0){
      #   return(Matrix(0, nrow=nrow(x)*ncol(x),ncol=1))
      # } else {
      #   as(c(as.matrix(x)), "Matrix") # FIXME: This should be better...
      # }
    } else {
      return(matrix(x@x))  
    }
  }
}

Vech <- function(x,diag=TRUE){
  if (!is(x,"matrix")){
    x <- as.matrix(x)
  }
  return(x[lower.tri(x,diag=diag)])
}