R/02_algebrahelpers_expectedmodel.R

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

expectedmodel <- function(x){
  # FIXME: If the distribution is Ising, just return:
  if (x@distribution == "Gaussian"){
    start <- parVector(x)
    prep <- prepareModel(start, x)
    for (g in 1:nrow(x@sample@groups)){
      x@sample@means[[g]] <- prep$groupModels[[g]]$mu
      x@sample@covs[[g]] <- prep$groupModels[[g]]$sigma
      
      if (length(x@sample@fimldata) > 0){
        nPat <- length(x@sample@fimldata[[g]])
        for (i in seq_len(nPat)){
          x@sample@fimldata[[g]][[i]]$means <- as.matrix(prep$groupModels[[g]]$mu[!x@sample@fimldata[[g]][[i]]$pattern,drop=FALSE])
          if (!all(x@sample@fimldata[[g]][[i]]$S == 0)){
            x@sample@fimldata[[g]][[i]]$S <- as.matrix(prep$groupModels[[g]]$sigma[!x@sample@fimldata[[g]][[i]]$pattern,!x@sample@fimldata[[g]][[i]]$pattern, drop = FALSE])
          }
        }
      }
    }
  }
  

  return(x)
}