R/03_modelformation_matrixSetup_beta.R

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

matrixsetup_beta <- function(
  beta, # beta argument
  nNode, # Number of nodes
  nGroup, # Number of groups
  labels,
  equal = FALSE,
  sampletable,
  name = "beta",
  start,
  onlyStartSign = FALSE
){
  # Fix beta:
  beta <- fixMatrix(beta,nGroup = nGroup,nrows = nNode,ncols = nNode,equal = equal,diag0=TRUE)
  
  # For each group, form starting values:
  if (missing(start)){
    betaStart <- 0.001*(beta!=0)  
  } else {
    betaStart <- beta
    for (g in seq_along(start)){
      betaStart[,,g] <- (beta[,,g]!=0) * start[[g]]
      if (onlyStartSign){
        betaStart[,,g] <- ifelse(betaStart[,,g]!=0,0.001 * sign(betaStart[,,g]), 0)
      }
    }
  }
  
  
  # Form the model matrix part:
  list(beta,
       mat =  name,
       op =  "<-",
       rownames = labels,
       colnames = labels,
       sparse = TRUE,
       start = betaStart,
       sampletable=sampletable
  )
}