mbco_parametric <- function(h0 = NULL,
h1 = NULL,
R = 10L,
alpha = .05,
checkHess = "No",
checkSE = "No",
optim = "NPSOL",
precision = 1e-9) {
OpenMx::mxOption(NULL, "Calculate Hessian", checkHess)
OpenMx::mxOption(NULL, "Standard Errors", checkSE)
OpenMx::mxOption(NULL, "Function precision", precision)
OpenMx::mxOption(NULL, "Default optimizer", optim)
res <- OpenMx::mxCompare(h1,
h0,
boot = TRUE,
replications = R
)
mbco_chisq <- res$diffLL[2] # Chi-square
mbco_pvalue <- res$p[2] # p-value
mbco_df <- res$df[2] - res$df[1] # df
mbcoTest <-
list(
chisq = mbco_chisq,
df = mbco_df,
p = mbco_pvalue
)
return(mbcoTest)
}
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.