I run a glm analysis on a dataset some days ago and obtained some relatively normal ranged (i.e from 0 to 2-3) odds ratios. I wanted to re-run the code because I had forgotten to print it and now I obtain very high and widely ranged odds ratios (i.e. OR 8 with CI from 2 to 60). Exactly same code, and same dataset. Why might this happen? I was thinking that I might have changed some scale into the R software, is there any way to verify that? Here is the code of my glm model:
MODEL <- glm(a ~ b+c+d+e, data = data, family = "binomial")