R/04_generalFit_fitfunction.R

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

# General fit function!
psychonetrics_fitfunction <- function(x, model){

  # Prepare model:
  if (model@cpp){
    prep <- prepareModel_cpp(x, model) # <- upated!
  } else {
    prep <- prepareModel(x, model)  
  }
  
  # What estimator to use:
  estimator <- model@estimator
  
  # What distribution to use:
  distribution <- model@distribution

  
  # Obtain estimator:
  if (model@cpp){
    estFun <- switch(
      estimator,
      "ML" = switch(distribution,
                    "Gaussian" = maxLikEstimator_Gauss_cpp, # <- updated!
                    "Ising" = maxLikEstimator_Ising_cpp # <- updated!
                    ),
      "ULS" =  switch(distribution,"Gaussian" = ULS_Gauss_cpp), # <- updated
      "DWLS" = switch(distribution,"Gaussian" = ULS_Gauss_cpp), # <- updated
      "WLS" = switch(distribution,"Gaussian" = ULS_Gauss_cpp), # <- updated
      "FIML" = switch(distribution,"Gaussian" = fimlestimator_Gauss_cpp) # <- updated!
    )
  } else {
    estFun <- switch(
      estimator,
      "ML" = switch(distribution,"Gaussian" = maxLikEstimator_Gauss, "Ising" = maxLikEstimator_Ising),
      "ULS" =  switch(distribution,"Gaussian" = ULS_Gauss),
      "DWLS" = switch(distribution,"Gaussian" = ULS_Gauss),
      "WLS" = switch(distribution,"Gaussian" = ULS_Gauss),
      "FIML" = switch(distribution,"Gaussian" = fimlEstimator_Gauss)
    )
  }
  
  
  # Run and return:
  estFun(prep)
}