High-accuracy spatial distribution estimation is crucial for cancer prevention and control. Due to their complicated pathogenic factors, the distributions of many cancers' mortalities appear blocky, and spatial heterogeneity is common. However, most of the commonly used cancer mapping methods are based on spatial autocorrelation theory. Sandwich estimation is a new method based on spatial heterogeneity theory. A modified sandwich estimation method suitable for the estimation of cancer mortality distribution is proposed in this study. The variances of cancer mortality data are used to fuse sandwich estimation results from various auxiliary variables, the feasibility of which in estimating cancer mortality distributions is explained theoretically. The breast cancer (BC) mortality of the Chinese mainland in 2005 was taken as a case, and the accuracy of the modified sandwich estimation method was compared with that of the Hierarchical Bayesian (HB), the Co-Kriging (CK) and the Ordinary Kriging (OK) methods. The accuracy of the modified sandwich estimation method was better than the HB, the CK and the OK methods, and the estimation result from the modified sandwich estimation method was more likely to be acceptable. Therefore, this study represents an attempt to apply the sandwich estimation method to the estimation of cancer mortality distributions with strong spatial heterogeneity, which holds great potential for further application.
