# Alternative test instead of logistic regression for binary dependent outcome variables?

I have a dataset of patients who underwent an operation, and I collected information on post-surgery complications such as necrosis (binary outcome variable). Now, I would like to investigate whether independent variables such as age (continuous), gender (binary), or underlying disease (categorical) influence the likelihood of experiencing complications.

My research suggests using a logistic regression model for binary outcome variables. Running it through R also gave me promising results. It just does not feel right. Is my approach the right one?

This my code in R:

glm(necrosis_y_n ~ age_op, data = df, family = binomial)


It was a statistically significant result with an odds ratio of factor $$1.053$$ per unit increase of age. It just does not feel like the right approach statistically.

• Can you tell us what you don't like about this? Logistic regression is absolutely the first thing that comes to mind. Nov 22, 2023 at 21:58
• There are other options for the link between the linear predictor and a binary outcome (besides the logit link used for logistic regression, R also allows for probit, cauchit, log, and complementary log-log), but the default logit/logistic regression typically works well. To improve the model, you might consider more flexible fitting of continuous predictors like age, for example with a regression spline.
– EdM
Nov 22, 2023 at 22:08
• A logistic regression might not be the best way to go, but it is the first idea that will come to mind for most. I second the inquiry about what feels wrong about running a logistic regression. Until we know what you don’t like about the logistic regression approach, I do not see a way to suggest alternatives.
– Dave
Nov 23, 2023 at 1:23
• Well, because patients with one complication tend to have other complications too. not because they necessarily influence each other, but because the underlying causative diagnosis is more severe so they are more likely to develop either of the complications. Furthermore, some of the outcome variables definitely are somewhat dependent. Dec 7, 2023 at 21:25
• So you want to model the presence of multiple types of complications?
– Dave
Dec 7, 2023 at 21:29

$$\text{log} \left(\frac{P}{1-P}\right) = \beta_0 + \beta_1 x_1,$$
where $$P$$ is the probability of an event, $$\beta_0$$ is the intercept (conditional average of event happening) and $$\beta_1$$ is the slope (which can be extended to additional terms like yours).