I am using the sem package to create a SEM model. Though the model runs, I get a warning:
Warning message:
In eval(expr, envir, enclos) :
Could not compute QR decomposition of Hessian.
Optimization probably did not converge.
When I call summary(model_test), it returns:
Error in summary.objectiveML(model_test) :
coefficient covariances cannot be computed
In addition: Warning message:
In vcov.sem(object, robust = robust, analytic = analytic.se) :
singular Hessian: model is probably underidentified.
My SEM model examines relationships between Market Orientation (MO), Government Incentives (GovInc), Entrepreneurial Orientation (EO), and Firm Performance (FirmPerf), with various constructs under each. Here’s the simplified code:
model_def <- specifyModel()
GovInfo -> GovInfo1, NA, 1
GovInfo -> GovInfo2, lamd1, NA
GovInfo -> GovInfo3, lamd2, NA
...
GovInfo -> GovInfo7, lamd6, NA
...
Perf -> Econ1, NA, 1
...
Perf -> Econ12, lamd23, NA
...
MO -> GovInfo, NA, 1
...
EO <-> FirmPerf, psi1, NA
...
cfa_test <- sem(model_def, data_matrix, N = 500)
summary(cfa_test)
I need guidance on understanding these errors and warnings, and how to address them. Your input would be very helpful.