# 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)
}
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.