The model including one binary outcome (0/1; incident rate ~1.2%), one main exposure, and 13 covariates. The whole model is significant and the goodness-of-fit is OK. However, model diagnostic is completely questionable (See Fig). Almost all the observations fall into the category of y=1 had residuals larger than 3. I wonder whether it is due to the low event rate of the outcome, and can anyone give me some advice as to how to deal with this problem? If the purpose of the model is inference rather than prediction, can I just ignore the model diagnostic results?