2
$\begingroup$

I have a dataset composed of 61 variables a qualitative one y=(0 or 1) and 60 other quantitative variables and 40000 observations. I want to do logistic regression, Lda, svm, rpart of the model y~.

When I use the vif function of package car it shows multicollinearity. The following is a part of the output:

>  vif(glm(y~., data=mydat,family=binomial))
>>   n1                 direc1            n2        direc2            n3 
>>  2.25               1.87               3.28        2.35           3.27 
>>  direc3           P1                  P2          P3                 P4 
>> 2.49              158.87             190.24        143.28           13.74 
>>    P5            P6                   P7             P8                 P9 
>>  119.93        212.23               616.43          146.59           169.48 

Should I keep only n1 ,direc1, n2 ,direc2 , n3 ,direc3 in the model? Or there is another solution to the problem ?


In fact I found this in some internet pages: "perfect separation is related to collinearity."

I want to address both problems. I read that collinearity between variables gives wrong coefficient estimates in a logistic regression model for example.

And perfect separation gives wrong coefficients estimates also. Really I am searching the best model prediction for a month but I did not find the good one till now because of these problems! ... It is ambiguous for me.

to solve multicolinearity problem i did a PCA but i did not now how to use the reults obtained in regression . In fact the image below did not show a separation between the two classes (0 and 1) . how can i interprete it ? which command in R allow me to get the new variables to use in regressions ? thanks a lot in advance for any help enter image description here

I read plenty of articles but could not find a solution!

$\endgroup$
  • 1
    $\begingroup$ Possible duplicate of How to prevent collinearity? $\endgroup$ – AdamO Feb 6 '17 at 13:21
  • 1
    $\begingroup$ I don't think it's a dupicate of that one; this is about what to do, not how to prevent. $\endgroup$ – Peter Flom Feb 6 '17 at 13:34
  • 4
    $\begingroup$ It is difficult to see how collinearity can cause perfect separation. What, then, is the problem you really want to address? Collinearity or perfect separation? $\endgroup$ – whuber Feb 6 '17 at 16:17
  • 1
    $\begingroup$ @whuber hi, in fact i found this in some internet pages "perfect separation is related to collinearity" i want to adress both of problems .Iread that collinearity between ariables gives wrong coefficients estimates in a logistic regression model for example . And perfect separation gives wrong coefficients estimates also really I am searching the best model prediction for a month but i did not found the good one till now because of these problems !... It is ambiguous for me Iread plenty of articles bu could not find a solution ! thanks a lot in advance for any help $\endgroup$ – prep Feb 6 '17 at 16:42
  • 2
    $\begingroup$ also see my question and answer here $\endgroup$ – hxd1011 Feb 6 '17 at 20:26
3
$\begingroup$

If you intend to use your model to predict normal/abnormal status in a new set of patients, you might not have to do anything about perfect separation or multicollinearity.

Say that you had one variable that perfectly predicted normal/abnormal status. A model based on that predictor would show perfect separation, but wouldn't you still want to use it?

Your perfect separation, however, might come from the large number of predictor variables, which might make perfect separation almost unavoidable. Then, even if no particular variables are perfectly related to disease state, you will have problems in numerical convergence of your model, and the particular combinations of variables that predict perfectly in this data set might not apply well to a new one. In that case, this page provides concise help on how to proceed. Also, this question and answer by @hxd1011 shows that ridge regression (under the name of "L2 regularization" on that page) can solve the problem of perfect separation.

This page is a good introduction to multicollinearity in the logistic regression context. Multicollinearity poses problems in getting precise estimates of the coefficients corresponding to particular variables. With a set of collinear variables all related to disease status, it's hard to know exactly how much credit each of them should get individually. If you don't care about how much credit to give to each, however, they can work very well together for prediction. In general, you typically lose in predictions if you throw away predictors, even predictors that don't meet individual tests of "statistical significance." Ridge regression tends to treat (and penalize) sets of correlated variables together, providing a principled approach to multicollinearity.

Ridge regression (as provided for example by the glmnet package in R) thus could solve both the perfect-separation and the multicollinearity problems, particularly if your interest is in prediction. This answer shows an example of using glmnet functions for logistic regression. Its example is for LASSO, but you simple set the parameter value alpha=0 for ridge instead. As another example, ISLR starting on page 251 has a worked through case of ridge regression for a standard linear model; specify the parameter family="binomial" instead for logistic regression.

In my experience, however, this type of model in clinical science sometimes isn't used for predicting new cases, but rather to try to argue that certain variables are the ones most closely related to disease status in general. I think that many of the comments on your question were getting at that possibility. The temptation is that the variables included in the "best model" for explaining the present data are then taken to be the most important in general. That can be a dangerous interpretation. Follow the feature-selection tag on this site for extensive discussion.


In an update to your question, you show that the first 2 principal components of your predictor matrix do not separate the 2 groups. That's not so surprising, as these are only the first 2 dimensions of a 60-dimension space, and it's hard to know exactly where among those dimensions the perfect separation arises. I don't think that your PCA helps here at all for variable selection. Ridge regression is the best way to try to proceed. Be warned, however, that if you are looking for p-values you do not get them directly from ridge regression. If p-values are important to you, repeat the process on multiple bootstrap samples of the data.

$\endgroup$
  • $\begingroup$ hi @EdM it seems that the two classes are superposed no ? or they are separated in other space but i could not found it no ? $\endgroup$ – prep Feb 7 '17 at 14:30
  • $\begingroup$ @prep Trying to do this manually with PCA is unlikely to find the space in which the perfect separation occurs. The separation has to do with the relation of the predictor variables to the disease status, while PCA only considers relations among the predictor variables. So the superposition of classes in PCA dimensions isn't surprising. Ridge regression, or maybe elastic net (uses both L1 and L2 regularization), will probably work best; elastic net does provide some variable selection if you really need to cut down on your number of predictors. $\endgroup$ – EdM Feb 7 '17 at 17:10
  • $\begingroup$ @prep Also note that ridge regression bears a close relation to principal components regression; principal components regression does an all-or-none selection of which components to include in the final model, while ridge regression down-weights the components with the least variance. See page 79 of ESLII and the figures and equations noted there. So ridge regression might automatically do what you are trying to do with your PCA plots, in all 60 dimensions at once. $\endgroup$ – EdM Feb 7 '17 at 17:16
  • $\begingroup$ thanks a lot @EdM so as I understood I should 1- apply ridge to the data to select the best variables 2- apply logistic regression, svm, nnet,to the selected variables and choose the best model resulted ? this is what should I do or not ? thanks a lot in adavnce for any help $\endgroup$ – prep Feb 7 '17 at 17:45
  • $\begingroup$ @prep ridge regression does not select variables; it uses all of them, weighted differently. A ridge logistic regression (provided by glmnet) directly provides the types of predictions that you want. An elastic net logistic regression (also available in glmnet) provides variable selection, but it might not be wise to use the variables selected that way as the ones to use in svm, nnet, etc.; I don't have much experience with svm or nnet approaches. Validating such a process would require repeating all steps including variable selection on multiple bootstrap samples. $\endgroup$ – EdM Feb 7 '17 at 18:26

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.