cancel
Showing results for 
Show  only  | Search instead for 
Did you mean: 
JMP is taking Discovery online, April 16 and 18. Register today and join us for interactive sessions featuring popular presentation topics, networking, and discussions with the experts.
Choose Language Hide Translation Bar
nmg
nmg
Level I

Help with mutivariable regression

Hello,

 

Im attempting to do a multivariable logistic regression looking at the independent association of variables for an adverse event, however the result for my variable of interest is reading as unstable. I believe this is because the variable actually has a 100% sensitvity for my outcome (adverse event) on univariate analysis. If this is this case, does anyone have suggestions on what to do?

 

Thank you!

Noelle

5 REPLIES 5

Re: Help with mutivariable regression

You have a binary response (adverse event present or not). You have multiple predictors in your model. Your estimates are unstable when you use nominal logistic regression. You might try instead to use GLM with the binomial distribution and logit link function. Be sure to select the Firth adjustment option as this will help if you have a separate problem.

nmg
nmg
Level I

Re: Help with mutivariable regression

Thank you for this answer!

 

I have a few follow up questions. When using GLM with binomial, is it alright that some of the my 2x2 tables have cells with values < 5?

Also, in this model, would you just report the estimate and SE and not an OR (like you would with nominal LR)?

One more question- does jmp have a way to report a RR if there are zeros in the 2x2? I have seen that some programs will add 0.5 or 1 to each cell.

 

Thank you again for your time

 

Noelle

nmg
nmg
Level I

Re: Help with mutivariable regression

Hello,

 

I am just wondering if anyone has any insight into this question. I have used GLM with the binomial distribution and logit link function and no longer have an "unstable" model. The primary predictor I am interested in is independently associated with my reponse variable (Adverse event vs No event) however now there are no CI reported. I have a total of 8 predictor variables for only 69 patients... When I remove some of the predictor variables from the model then the confidence intervals appear. My questions are...

 

1. Are the CI initially not being reported because there are too many predictor variables in the model

 

2. Is this issue still stemming from having 100% sensitivity for the outcome (the report states "evidence of perfect fit for some data points detected, and the Hessian matrix suggests quasi-complete separation of the data. Fit and results are of questionable value- proceed with caution".

 

3. Is there another more appropriate model I should be using?

 

Thank you again for all of you help

 

 

 

Re: Help with mutivariable regression

Just checking a few points:

  • You select Fit Model and changed the fitting personality to Generalized Linear Models.
  • You selected Binomial for the distribution and the default canonical Logit link function.
  • You selected the Firth adjustment option.

Some thoughts are:

  1. I believe that the CI are missing because of the errors that were reported and the fact that some of your LRT statistics are zero. Also, in cases where the relative likelihood function is too flat, it can't be used for the CI.
  2. Yes. If you eliminate those problems (e.g., removing some effects from the model) then the missing CI might be due to the problem with the likelihood function.
  3. You could try Generalized Regression methods if you have JMP Pro.

It could also be a problem with the data or the assumptions of the model and estimation method that you are using. Can't say for sure at this point.

Re: Help with mutivariable regression

I forgot to say that you might have a problem because of rare events, that is, response levels that are not observed for some of your predictor levels. It is related to the separation problem. The Firth adjustment to the estimation helps, but it can't always eliminate all of the problems if the rare event problem is severe.