What happens if you run modelE with mxRun() instead? Also, what are the free parameter values of modelE before it is run (you can just apply summary() or coef() to modelE)?
You could also try doing incomplete "test runs" of modelE, for example:
Thanks for the answer. I applied coef () and summary () to model E. The findings are as follows and I tried using Mxrun as you mentioned, but I got an error like this. What do I have to do in this situation?
Thank you so much for your time and help!
> fitE <- mxRun( modelE, onlyFrontEnd=T)
Error: mxRun does not accept ... arguments. The first parameter in ... was named 'onlyFrontEnd'with value 'TRUE'> fitE <- mxRun( modelE, useOptimizer=F)
Running twoEvj with6 parameters
Error: The job for model 'twoEvj' exited abnormally with the error message: fit is not finite (Ordinal covariance is not positive definite indata'DZ.data'row64(loc1))>coef(modelE)
meanno_2_FA t1thno_2_FA t1thcannabis_dummy VA11 VA21 VA22
0.510300807890.980878911800.959140079260.000610607070.007859280831.76357137423
VE11 VE21 VE22
0.00245181741-0.00282031911-0.76357137423Summary of twoEvj
The model does not satisfy the first-order optimality conditions to the required accuracy, and no improved point for the merit function could be found during the final linesearch (Mx status RED)
Your ordinal model may converge if you reduce mvnRelEps
try: model <- mxOption(model, 'mvnRelEps', mxOption(model, 'mvnRelEps')/5)
free parameters:
name matrixrowcol Estimate lbound ubound
1 meanno_2_FA meanG 110.510300807892 t1thno_2_FA thinG 110.98087891180-33 t1thcannabis_dummy thinG 120.95914007926-34 VA11 VA 110.000610607075 VA21 VA 120.007859280836 VA22 VA 221.763571374237 VE11 VE 110.002451817418 VE21 VE 12-0.002820319119 VE22 VE 22-0.76357137423
confidence intervals:
lbound estimate ubound note
twoAEvj.US[1,7]-0.0796416110.199386820.44356784
twoAEvj.US[1,9]0.0000000000.000000000.00000000!!!
twoAEvj.US[1,11]0.5564322140.800613181.07964766
twoAEvj.US[2,7] NA 1.55970243 NA !!!
twoAEvj.US[2,9]0.0000000000.000000000.00000000!!!
twoAEvj.US[2,11] NA -0.55970243 NA !!!
twoAEvj.US[2,8]1.3783592391.763571551.93985753
twoAEvj.US[2,10]0.0000000000.000000000.00000000!!!
twoAEvj.US[2,12]-0.939845346-0.76357155-0.37836482
To investigate missing CIs, run summary() again, with verbose=T, to see CI details.
Model Statistics:| Parameters | Degrees of Freedom | Fit (-2lnL units)
Model:9407-404.64357
Saturated: NA NA NA
Independence: NA NA NA
Number of observations/statistics:105/416
Constraint 'Var1' contributes 1 observed statistic.
Information Criteria:|df Penalty | Parameters Penalty | Sample-Size Adjusted
AIC:-1218.6436-386.64357-384.74884
BIC:-2298.8054-362.75793-391.19060
CFI: NA
TLI:1(also known as NNFI)
RMSEA:0[95% CI (NA, NA)]
Prob(RMSEA <=0.05): NA
> fitE <- mxRun( modelE, onlyFrontEnd=T)
Error: mxRun does not accept ... arguments. The first parameter in ... was named 'onlyFrontEnd'with value 'TRUE'
My mistake. The argument is onlyFrontend, not onlyFrontEnd.
Based on what I see in your post, try setting 'VE21' to zero before running the E-only model.
Thanks, it worked. Another question I have is this: I am trying to add another continuous variable to the same model. In this case, I will have two continuous and an ordinal variable. I could not find a script with multiple variables containing both ordinal and continuous variables.
Can I modify it to this model? If possible, should the oc<- (0,0,1) and SvMe <- (0, 10,15) variables be 3 digits like this?
I'm just trying to learn, sorry for so many questions.
Printer-friendly version
What happens if you run
modelEwithmxRun()instead? Also, what are the free parameter values ofmodelEbefore it is run (you can just applysummary()orcoef()tomodelE)?You could also try doing incomplete "test runs" of
modelE, for example:Thanks for the answer. I applied coef () and summary () to model E. The findings are as follows and I tried using Mxrun as you mentioned, but I got an error like this. What do I have to do in this situation?
Thank you so much for your time and help!
My mistake. The argument is
onlyFrontend, notonlyFrontEnd.Based on what I see in your post, try setting 'VE21' to zero before running the E-only model.
Thanks, it worked. Another question I have is this: I am trying to add another continuous variable to the same model. In this case, I will have two continuous and an ordinal variable. I could not find a script with multiple variables containing both ordinal and continuous variables.
Can I modify it to this model? If possible, should the oc<- (0,0,1) and SvMe <- (0, 10,15) variables be 3 digits like this?
I'm just trying to learn, sorry for so many questions.
Thank you in advance!