of that condition/ diagnosis was reported. A series of submodels examining the statistical significance of A, C and E were compared to the full model to derive the best-fitting trivariate model. Submodels were tested by calculating the difference between the −2 log likelihood fit of the full model and the nested submodel, which is distributed as χ2 for the given degrees of freedom.