I'm using the robustbase package to run a glm estimation. However when I do it, I get the following error:

Error in solve.default(crossprod(X, DiagB * X)/nobs, EEq) : 
  system is computationally singular: reciprocal condition number = 1.66807e-16

What does this mean/indicate? And how can I debug it?

PS. If you need anything (the formula/specification or data) to answer, I'll gladly provide it.


It means your design matrix is not invertible and therefore can't be used to develop a regression model. This results from linearly dependent columns, i.e. strongly correlated variables. Examine the pairwise covariance (or correlation) of your variables to investigate if there are any variables that can potentially be removed. You're looking for covariances (or correlations) >> 0. Alternatively, you can probably automate this variable selection by using a forward stepwise regression.

This can also result from having more variables than observations, in which case your design matrix is probably not full rank. This is a bit trickier to fix, but there are ways. I believe lasso regression is supposed to work well when the data is "wider" than it is "long."

Keep in mind: if you decide to try lasso or stepwise selection, your doing much more (in terms of variable selection) than just handling multicolinearity.

  • $\begingroup$ Thanks, two questions: 1. Do you know of any procedure in r, that might work well for me in debugging this? 2. I can see, that when I run the exact same model on the exact same data with glm() instead of glmrob() I don't get any errors - does that make sense, or does it imply an error in the glmrob() function? $\endgroup$ – NK1 Nov 14 '13 at 9:11

Spent a long time on this error but checkout this page which solved it for me


Basically you can just change the method mice uses. I was using ppn and changed it to cart

  • 1
    $\begingroup$ Try to make the answer self-contained within the text of your response, please. That way, people are more likely to read it and less likely to delete it as too short and not on-topic. $\endgroup$ – Carl Sep 24 '17 at 3:02

Your Answer

By clicking “Post Your Answer”, you agree to our terms of service, privacy policy and cookie policy

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