I have fitted a structural equation model with 28 parameters (500 samples, 12 observed variables). there is a discrepancy between the output of MPlus to that of my own code for estimating parameter CIs.
to estimate CI for a parameter, that value of the parameter that increased the model chi-square value by 3.84 units (compared to value at original parameter estimate) was taken as the 95% confidence interval. model chi-square value was estimated as (N-1)*val, where N was sample size and val was value of the fit function (ML).
would be grateful if you could throw light on possible reasons for discrepancy in output between MPlus and my code.