Disease mapping method comparing the spatial distribution of a disease with a control disease

Small-area methods are being used in spatial epidemiology to understand the effect of location on health and detect areas where the risk of a disease is significantly elevated. Disease mapping models relate the observed number of cases to an expected number of cases per area. Expected numbers are often calculated by internal standardization, which requires both accurate population numbers and disease rates per gender and/or age group. However, confidentiality issues or the absence of high-quality information about the characteristics of a population-at-risk can hamper those calculations. Based on methods in point process analysis for situations without accurate population data, we propose the use of a case-control approach in the context of lattice data, in which an unrelated, spatially unstructured disease is used as a control disease. We correct for the uncertainty in the estimation of the expected values, which arises by using the control-disease's observed number of cases as a representation of a fraction of the total population. We apply our methods to a Belgian study of mesothelioma risk, where pancreatic cancer serves as the control disease. The analysis results are in close agreement with those coming from traditional disease mapping models based on internally standardized expected counts. The simulation study results confirm our findings for different spatial structures. We show that the proposed method can adequately address the problem of inaccurate or unavailable population data in disease mapping analysis.