Logistic Macro Output Checklist
- Are all predictors that you expected in the models?
- Does each predictor have the categories that you expected? Are the reference categories correct?
- Is the odds ratio in the direction that you expected, ie less than 1 or greater than 1? If not, double-check for reversed coding (of outcome or predictor) or use of the wrong reference category.
- Are any odds ratios extremely small or extremely large or blank? Blanks indicate detected cases of degenerate estimates, where p-values and confidence limits are not accurate. More valid P-values can be obtained by likelihood ratio tests or exact methods. CI's can be obtained by exact methods or profile likelihood (not currently implemented in SAS). There is a user-written macro called FL that uses an alternative approach for degenerate estimates; this is more convenient but not commonly used in clinical research.
- Do estimates of the odds ratio for a particular variable of interest vary a lot depending on what other variables are in the model? If so, check whether this is due to deletion of observations with missing values or due to potential confounding. To assess this, run models with and without the other variable of interest, but restricted in both cases to observations where the variables is not missing.
- Are there any small p-values for linearity checking? If so, consider modeling the predictor in categories, logarithmically transforming it, or adding additional terms. If the potentially non-linear term is not of primary interest, then additional terms can be quadratic, cubic, etc. (i.e., polynomial terms). If the variable is of primary interest, then linear spline terms may be more interpretable. If polynomial terms are used, graphical summaries of the fitted effect of the variable will be needed, as the individual coefficients are not readily interpretable.
- Is the H-L p-value small? If so, this suggests non-linearity of one or more predictor or existence of one or more interactions among the predictors. Address non-linearity as for point 6 and/or add in plausible interaction terms. If the H-L p-value cannot be raised above 0.05, then consider adjusting the SE's, CI's, and p-values for overdispersion (the Genmod procedure can do this automatically).
-- PeterBacchetti - 23 Mar 2009