Tell me more ×
Cross Validated is a question and answer site for statisticians, data analysts, data miners and data visualization experts. It's 100% free, no registration required.

Structural equation models sometimes cannot be fit due to a singular sample covariance matrix. Now some authors suggest to apply a "gentle ridging" to the diagonal which helps (this is for example done in Mplus).

  1. How does ridging exactly work? and
  2. If your preferred SEM software does not support ridging, how can you do it "by hand"?
share|improve this question

2 Answers

I can kill two birds with one stone by explaining how you get an ols software to fit ridge regression for a fixed ridging parameter. You simply augment the data set you have with one observation for each beta you intend to ridge but set the covariate pattern to be all zero except for the ridging variable which gets set to the square root of the ridging parameter.

This includes the intercept. So you need to define your own intercept variable as it will have zeros in the data augmented version.

Also you set the data augmented "response" variable to whatever you want the betas to be shrunk towards. Typically these are all zero.

If you do the maths on the augmented data set, you'll get the ridged beta estimate.

You will also find that the augmented x matrix is now full column rank. I can add the maths if you want but I don't think you need it.

Note that this only helps with understanding for fixed value of the ridging parameter. Choosing a good value is another issue.

share|improve this answer
The paucity of capital letters leads me to suspect you're posting from a mobile device, again. :-) In Hastie, et al., there is a problem which leads the reader through this data-augmentation interpretation. – cardinal Jul 30 '12 at 13:51
Your suspicion is well founded. There are now capitals :-) – probabilityislogic Jul 30 '12 at 23:32

If you have a singular sample covariance matrix, it means that your variables are perfectly collinear, or one of the variables does not vary. I would say you have to reconsider your model, rather than to try to tweak something in the black box called Mplus to make it kinda work.

share|improve this answer

Your Answer

 
discard

By posting your answer, you agree to the privacy policy and terms of service.

Not the answer you're looking for? Browse other questions tagged or ask your own question.