 Research
 Open access
 Published:
Mapping the prevalence of cancer risk factors at the small area level in Australia
International Journal of Health Geographics volume 22, Article number: 37 (2023)
Abstract
Background
Cancer is a significant health issue globally and it is well known that cancer risk varies geographically. However in many countries there are no small arealevel data on cancer risk factors with high resolution and complete reach, which hinders the development of targeted prevention strategies.
Methods
Using Australia as a case study, the 2017–2018 National Health Survey was used to generate prevalence estimates for 2221 small areas across Australia for eight cancer risk factor measures covering smoking, alcohol, physical activity, diet and weight. Utilising a recently developed Bayesian twostage small area estimation methodology, the model incorporated surveyonly covariates, spatial smoothing and hierarchical modelling techniques, along with a vast array of small arealevel auxiliary data, including census, remoteness, and socioeconomic data. The models borrowed strength from previously published cancer risk estimates provided by the Social Health Atlases of Australia. Estimates were internally and externally validated.
Results
We illustrated that in 2017–2018 health behaviours across Australia exhibited more spatial disparities than previously realised by improving the reach and resolution of formerly published cancer risk factors. The derived estimates revealed higher prevalence of unhealthy behaviours in more remote areas, and areas of lower socioeconomic status; a trend that aligned well with previous work.
Conclusions
Our study addresses the gaps in small area level cancer risk factor estimates in Australia. The new estimates provide improved spatial resolution and reach and will enable more targeted cancer prevention strategies at the small area level. Furthermore, by including the results in the next release of the Australian Cancer Atlas, which currently provides small area level estimates of cancer incidence and relative survival, this work will help to provide a more comprehensive picture of cancer in Australia by supporting policy makers, researchers, and the general public in understanding the spatial distribution of cancer risk factors. The methodology applied in this work is generalisable to other small area estimation applications and has been shown to perform well when the survey data are sparse.
Background
In 2020, an estimated 19.3 million people were diagnosed with cancer worldwide [1], causing a huge health burden. Moreover, incidence of cancer has been shown to exhibit strong spatial disparities, which due to improved models and better data accessibility are now communicated to the public via interactive Atlas platforms. In Australia, a notable Atlas is the Australian Cancer Atlas (ACA) [2], which provides interactive maps of small area level estimates of incidence and relative survival rates for a wide range of cancer types.
Whiteman et al. [3] suggest that at least one in every three cancers in Australia can be attributed to modifiable risk factors such as tobacco smoking, obesity, poor diet, insufficient physical activity, excessive sun exposure and alcohol consumption. Understanding the prevalence of cancer risk factors is pivotal to cancer prevention.
To better assess how cancer risk factors vary by location and target interventions, many countries have generated small area estimates for their prevalence including Australia [4], the US [5], Canada [6], Iran [7], and Luxembourg [8]. When generating small area estimates, practitioners must consider the reach and resolution of their results. Reach refers to the proportion of the small areas for which estimates are available, while resolution pertains to the small area population and geographical sizes. While the need for high resolution relates to minimizing outcome heterogeneity in larger areas and populations, the need for complete coverage (or high reach) ensures policy makers have complete spatial information. If small area estimates suffer from low reach or resolution the effectiveness of targeted interventions could be affected.
In Australia, the Social Health Atlases of Australia (SHAA) [4] is the major platform providing nationwide estimates for cancer risk factors at a small area level. The estimates were derived from the 2017–2018 National Health Survey (NHS). However the reach and resolution of the SHAA estimates could be improved. The larger areal units used in the SHAA combine heterogeneous subpopulations, resulting in estimates that are averages over different populations. The limitation regarding reach meant that no estimates are provided for very remote areas. Given that health disparities tend to widen with increasing remoteness [9,10,11], generating estimates for these areas is important for targeted public health initiatives in Australia. The modelled estimates provided by the SHAA use the best data source available, so the problem cannot be solved by using a different dataset or collecting better data; the solution is to use new methods of small area estimation (SAE) [12].
SAE is a wellestablished survey method that leverages auxiliary data, such as census data, to estimate parameters of interest for small geographic areas with limited or no survey data. Modelbased SAE methods, which borrow strength across areas [13], can be applied at either the area [14] or individual level [15], with the latter requiring access to survey and census microdata.
Proportion arealevel models are commonly used [16,17,18,19]; however, they become unsuitable when some of the input data (arealevel proportion estimates) are unstable, i.e. exactly zero or one [20]. Sparse survey data and modelling rare or common population characteristics exacerbate this instability [21]. Solutions to instability include perturbing direct estimates prior to modelling [22] or excluding unstable areas [17]. Alternatively, modelling at the individual level, such as through multilevel regression and poststratification (MrP) [23], can be pursued. However, the use of individual level SAE models to derive proportion estimates are limited by the need for census microdata [24], which restricts the choice of covariates. Note that the modelling for the SHAA was conducted by the Australian Bureau of Statistics (ABS). Unfortunately, given that the published details of the ABS approach are modest [25], we can only infer the use of a individual level model.
While individual and area level models have limitations, recent work supports the utility of twostage SAE approaches, which involve separate modelling at both levels [21, 26,27,28]. Twostage approaches have many benefits that are particularly relevant for this application as they can alleviate unstable direct estimates by smoothing individual level outcomes, accommodate even severely sparse survey data thanks to multistage smoothing, and utilize more auxiliary data (e.g. surveyonly covariates), permitting more flexible models and better predictions.
In this work, we generate small area level prevalence estimates for eight cancer risk factor measures using the Bayesian twostage small area estimation methodology we developed for sparse survey data [21]. Our method considers a variety of data sources, including individual level survey data and area level auxiliary data such as census, remoteness and socioeconomic data. To assess the quality of our estimates, we used a dual validation approach whereby most SA2s are benchmarked to the substate level using fully Bayesian benchmarking [29], with the remaining SA2s (predominantly very remote areas) undergoing external validation. The results of this work will complement the current small area level estimates of cancer incidence and relative survival already available in the ACA [2].
Data
Geographical areas
Geographical location was defined according to the 2016 Australia Statistical Geography Standard (ASGS) [30]. We generated prevalence estimates at the statistical area level 2 (SA2) level, which is the lowest level of the ASGS hierarchy for which detailed census population characteristics are publicly available. SA2s are recognized as achieving the optimal balance between privacy and resolution [31]. Note that the SHAA provides estimates at a lower resolution, using population health areas (PHAs) which are constructed from single or multiple SA2s (40% and 39% of PHAs are constructed from one and two SA2s, respectively). In 2016 Australia had 1165 PHAs and 2310 SA2s [32] with median population sizes of 7500 and 15000 for SA2s and PHAs, respectively.
Throughout this analysis, we also used statistical area level 3 (SA3) and statistical area level 4 (SA4). By virtue of the hierarchical nature of the ASGS, SA2s are nested within SA3s, and SA3s are nested within SA4s. There is a median of 6 and 22 SA2s nested within each SA3 (n = 333) and SA4 (n = 88), respectively.
Of the 2310 SA2s covering Australia, SA2s with no physical location (comprising “MigratoryOffshoreShipping” and “No usual address” codes for each State and Territory) (n = 18), very remote island SA2s (Christmas Island, Cocos Island, Norfolk Island and Lord Howe Island) (n = 4), and SA2s with annual average population \(\le\) 10 (n = 67) were excluded. This left 2221 SA2s to use in the modelling. The remaining SA2s had a median (interquartile range (IQR)) population of 7859 (4483, 12753). Note that although Jervis Bay is classified as an “Other Territory” by the ABS, we included it as part of the state New South Wales.
Data sources
Survey data
The individual level survey data and sampling weights were obtained from the 2017–18 National Health Survey (NHS), which is an Australiawide populationlevel health survey conducted every 3–4 years by the ABS [33, 34]. This survey excluded very remote areas of Australia (\(\approx 0.8\)% of 2016 population), discrete Aboriginal and Torres Strait Islander communities (\(\approx 0.5\)% of 2006 population as per the ABS Community Housing and Infrastructure Needs Survey conducted only in 1999, 2001, and 2006 [35]), and nonprivate dwellings (\(\approx 2\)% of 2016 population [36]). Nonprivate dwellings include hotels and motels, hostels, boarding schools and boarding houses, hospitals, nursing and convalescent homes, prisons, reformatories and single quarters of military establishments and shortstay caravan parks. The ABS highlights that these exclusions should only have a minor effect on aggregate estimates for the states and territories of Australia.
The 2017–18 NHS data consist of 17248 sampled persons 15 years and older, with 878 persons under the age of 18. The data cover 1694 (76%) of the 2221 SA2s across Australia (see Fig. 1) and provide a median (IQR) SA2 level sample size of 8 (5, 13). The median SA3 and SA4 level sample sizes were 42 (25, 65) and 154 (101, 226), respectively. The NHS was also used to obtain daily smoking rates at the SA4 level [37]. Other sources of Australian health data are described in Section A of the Additional File 1.
Population data
Estimated Resident Population data stratified by 5year age groups (15 years and above), sex and SA2, were obtained from the ABS for both 2017 and 2018 [38]. In this study the SA2 level population counts were derived by averaging across the two years. One of the risk factors (risky waist circumference) is only appropriate for ages 18+ and so modelling excluded persons under 18. Assuming that the singleyear age distribution in this age group was uniform, we estimated that the population of 18–19year olds was 40% of the 15–19year old population.
For the SA2 level auxiliary data, we used data from the 2016 Australian census, represented as proportions (for categorical data) or averages (for continuous data) of individuals in each SA2. Census data for age, sex, nonschool level education (higher education), highschool completion status, occupation, labour force status, personal weekly income, religious affiliation, registered marital status, First Nations Australian status, and household composition were obtained from the ABS [39]. These census factors made up 84 separate variables. Like Chidumwa et al. [40], to reduce the dimension of these socioeconomic and demographic data we used Principal Components (PC) Analysis, where we retained the first six principal components as they accounted for approximately 62% of the variation (see Section C.2 of the Additional File 1 for more details).
Other data
Australian research suggests that cancer burden [31, 41], and the prevalence of cancer risk factors varies strongly by remoteness and socioeconomic status (SES) [9, 42]. Data on SA1 level remoteness were provided by ABS and based on the Accessibility and Remoteness Index of Australia (ARIA+) [43], and converted to SA2 using population proportions. Remoteness is divided into five groups  major cities, inner regional, outer regional, remote, and very remote  based on a measure of relative geographic access to services. Given that very remote areas of Australia were intentionally excluded during data collection for the 2017–18 NHS, we followed the approach of Das et al. [37], and collapsed the outer regional, remote and very remote categories to a single remoteness group. Of the SA2s with sample data, 69% were major cities. The SA2 sample sizes tended to be larger for outer regional to very remote areas (median of 11 and IQR of 6 to 21) than for major city areas (median of 8 and IQR of 4 to 12).
SA2 level SES was sourced from the ABS SocioEconomic Indexes for Areas product [44]. Like other Australian health studies [37, 42, 45], we used the Index of Relative SocioEconomic Disadvantage (IRSD). The IRSD is a general SES index constructed using principal components analysis that summarises the economic and social conditions of individuals and households within a given area in order to determine the area’s overall relative disadvantage. A low IRSD score indicates a large proportion of relatively disadvantaged individuals in a given SA2 [45].
In this work, we used IRSD national deciles as a categorical variable with 10 groups, where 1 represents the most disadvantaged or lowest SES group and was used as the reference group. Although the IRSD can be used as a continuous variable, it is recommended to use deciles [44], and this also gave superior model performance. There were 44 of the 2221 SA2s without an IRSD value provided, so these had the closest IRSD decile assigned according to their corresponding PC1 (principal component 1) values.
We also obtained prevalence estimates and measures of uncertainty for risky alcohol consumption (more than 2 standard drinks a day on average), adequate fruit intake, obesity, overweight, current smokers and inadequate physical activity from the SHAA [4] at the Primary Health Network (PHN) and PHA level for adults. These data were downloaded as agestandardised rates per 100 people with 95% confidence intervals. Definitions and details are available in Section C and D of the Additional File 1 and the online SHAA platform [4].
Risk factors
Broad risk factor groups were selected by consulting three sources: a wide range of experts in the fields of public health, epidemiology and oncology; literature, specifically evidence for casual associations [46] with, and population attributable fractions [3, 41, 47] for, cancer incidence; and the availability of data in the 2017–18 NHS. In this work we selected the following five broad risk factor groups: tobacco smoking, alcohol, diet, weight and physical activity. According to the 2015 Australian Burden of Disease study [41] these were attributable to 22.1%, 4.5%, 4.2%, 7.8% and 2.9% of the total cancer burden, respectively.
We explored a variety of possible measures and corresponding definitions for each of the five broad risk factor groups, placing priority on the definitions and recommendations used in the SHAA [4], the work by Whiteman et al. [3], Cancer Council Australia [48] and those provided by Australian government agencies such as Cancer Australia [45], the National Health and Medical Research Council (NHMRC), the Australian Institute of Health and Welfare (AIHW) [49] and the Australian Department of Health and Aged Care (DOH).
Table 1 summarizes the five broad risk factor groups and the eight corresponding measures and definitions. Table 2 gives direct estimates for these measures stratified by the eight states and territories of Australia. The risk factor measures proposed are designed to be crosssectional and strike a natural balance between being specific to cancer while maintaining applicability to a variety of other health conditions [41]. Note that some risk factor groups, for example weight, required several differing measures.
We defined the risk factor measures as binary where a survey individual received a value of one if they did not meet guidelines, or were in the unhealthy category. Unlike the SHAA which provides agestandardised rates by PHAs [4], we used proportions (prevalence) due to their common use in both the literature [12, 16, 53] and other digital Atlases [5]. Furthermore, deriving agestandardised rates requires prevalence estimates by area and age. This level of disaggregation is possible at the PHA level, but not feasible at the SA2 level.
We provide further details and the motivation for the selected risk factor measure definitions in Section B of the Additional File 1.
Statistical models
Bayesian model
Given the sparse nature of the available data for this SAE analysis, we used the Bayesian twostage logistic normal (TSLN) approach we proposed recently [21]. Our previous study showed that the TSLN approach could outperform commonly used area [17, 18, 54] and individual level [55] models both in a simulation study focusing on sparse survey data and an application using the 2017–18 NHS data. The twostage structure of the TSLN approach includes an individual level stage 1 model, followed by an area level stage 2 model.
The same TSLN approach, with very similar components, was chosen to be applied to all eight risk factor measures. The selection of fixed and random effect structures for the two models was guided by the goal of achieving a balance between parsimony across risk factor measures and predictive performance. We followed the advice by Goldstein [56] and initially used frequentist algorithms to select fixed and random effects, with fully Bayesian inference via Markov chain Monte Carlo (MCMC) for final model checking. Further details regarding model selection are given in Section E of the Additional File 1.
Let \(y_{ij} \in \{0,1\}\) be the binary value from the NHS for sampled individual \(j = 1, \dots , n_i\) in SA2 \(i = 1, \dots , m\), where \(n_i\) is the sample size in SA2 i. Further, let \(m = 1694\) and \(M = 2221\) be the number of sampled and total number of SA2s, respectively. The goal of this analysis is to generate estimates of the true proportions of each risk factor measure, \(\varvec{\mu } = \left( \mu _1, \dots , \mu _M \right)\).
In this analysis, we used two versions of the survey weights, \(w^{\text {raw}}_{ij}\), provided by the ABS [55, 57] to correct for sampling bias and promote designconsistency. The first, \(w_{ij}\), was used for direct estimation and the second, \(\tilde{w}_{ij}\), was used in the stage 1 model (see Section C.1 in the Additional File 1). Using the survey weights, small area proportion estimates can be computed using the Hajek [58] direct estimator,
with an approximate sampling variance of [54, 59],
Direct estimators, such as Eqs. (1) and (2), have low variance and are designunbiased for \(\mu _i\) when \(n_i\) is large, but have high variance when \(n_i\) is small [13].
Stage 1: Individual level model
The stage 1 model is a Bayesian pseudolikelihood logistic mixed model. Let \(\pi _{ij}\) be the probability of \(y_{ij} = 1\) for sampled individual j in SA2 i. Following the notation of Parker et al. [55], we represent the pseudolikelihood for a probability density, p(.), as \(p\left( y_{ij} \right) ^{\tilde{w}_{ij}}\). Pseudolikelihood is used to ensure the predictions from the logistic model are approximately unbiased under the sample design [60, 61]. Thus, the stage 1 model likelihood is given by,
where \(\text {logit}\left( \pi _{ij} \right)\) is modelled using a generic linear predictor that is applicationspecific. In this work, we used several unique components summarised in Fig. 2. The linear predictor included eight individual level categorical covariates and seven area level covariates as fixed effects. Unstructured individual and SA2 level random effects were also applied. In addition to these, borrowing ideas from MrP [62], we included two hierarchical random effects based on categorical covariates that were themselves derived from the interaction of numerous individual level demographic and health covariates. A discussion of the priors used is given on the subsequent page. More details can be found in Section C of the Additional File 1.
Stage 2: Area level model
After fitting the stage 1 model, the individual level predictions are aggregated to the area level, producing stage 1 (S1) proportion estimates \(\hat{\mu }^{\text {S1}, (t)}_i\) using Eq. (1), and sampling variances, \(\psi ^{\text {S1}, (t)}_i = \widehat{\text {v}} \left( \hat{\mu }_i^D \right) + \widehat{\text {v}} \left( \hat{B}^{(t)}_i \right)\), for all posterior MCMC draws, \(t=1, \dots , T\) [21], where the function to compute the sampling variance, \(\widehat{\text {v}}(.)\), is given in Eq. (2) and \(\hat{B}^{(t)}_i = n_i^{1} \left( \sum _{j=1}^{n_i} w_{ij} \left( \pi ^{(t)}_{ij}  y_{ij} \right) \right)\) quantifies the level of smoothing achieved by using \(\pi _{ij}\) instead of \(y_{ij}\).
Using the common logistic transformation [18, 54], let
thereby permitting the use of a Gaussian likelihood in the second stage model. Let \(\bar{\tau }_i^{\text {S1}}\) be the empirical posterior mean of \(\tau _i^{\text {S1}}\) and \(\widehat{\text {v}} \left( \hat{\theta }_i^{\text {S1}} \right)\) be the empirical posterior variance of \(\hat{\theta }_i^{\text {S1}}\). Finally, by selecting a random subset of the posterior draws, say \(\widetilde{T}\), let \(\hat{\varvec{\theta }}^{\text {S1}}_i = \left( \hat{\theta }^{\text {S1}, (1)}_i, \dots , \hat{\theta }^{\text {S1}, (\widetilde{T})}_i \right)\).
The stage 2 model is a Bayesian spatial FayHerriot [14] model. Unlike previous twostage approaches [26, 27], we accommodate some of the uncertainty inherent in fitting the stage 1 model by using the vector \(\hat{\varvec{\theta }}^{\text {S1}}_i\) as input to the stage 2 model. The stage 2 model likelihood for the posterior draws from the stage 1 model is,
where \(\theta _i\) is modelled using a generic linear predictor that is problem specific. The final proportion/prevalence estimate for the ith SA2, denoted \(\mu _i\), is given by the posterior distribution of \(\text {logit}^{1} \left( \theta _i \right)\). To ensure that posterior uncertainty remains unaffected by the choice of \(\widetilde{T}\), we downscale the likelihood contribution by \(1/\widetilde{T}\).
In this work, we used several unique components for the linear predictor of \(\theta _i\) which are summarised in Fig. 3. The linear predictor included the SES index deciles and remoteness as standard fixed effects. In addition, PC1 to PC6 were used as fixed effects with coefficients varying according to remoteness. The linear predictor also included an external latent field constructed from the SHAA’s estimates and a BYM2 spatial random effect [63] at the SA2 level. Given we did not include SA3 level census covariates, an unstructured random effect at the SA3 level was employed. To smooth unstable variances we used the generalized variance function [12, 64, 65] described in Section C.4.6 of the Additional File 1. More details can be found in Section C of the Additional File 1.
Priors
The Bayesian models described above are completed by the specification of priors. Given the complexity of the two models, in this work generic weakly informative priors were adopted based on preliminary analysis of the data [66]. In both models, all fixed effect coefficients were given \(\text {N}\left( 0, 2^2 \right)\) priors with intercepts given a student\(t\left( 0, 2^2, \text {df} = 3 \right)\). We used \(\text {N}^{+}\left( 0, 1^2 \right)\) and \(\text {N}^{+}\left( 0, 2^2 \right)\) priors for all standard deviation terms in the stage 1 and stage 2 models, respectively. The mixing parameter in the BYM2 [63] random effect was given a \(\text {Uniform}\left( 0,1 \right)\) prior (see Section C of the Additional File 1).
We conducted sensitivity analysis by using more, \(\text {N}\left( 0, 1^2 \right)\), and less, \(\text {N}\left( 0, 100^2 \right)\), informative priors for the fixed effects in both models. We also experimented with using exponential priors with rates of 0.5 and 1 for standard deviation terms. Finally, we examined model fit when using an informative Beta prior for the mixing parameter. We found that the model fit and prevalence estimates were unaffected by these prior changes. The chosen priors gave superior sampling efficiency and convergence.
Validation
For validation of the small area estimates, we adopted a dual approach, using both internal and external methods. See Section C.5 in the Additional File 1 for details.
Internal benchmarking
Internal validation involved a fully Bayesian benchmarking procedure [29] that adjusts the results obtained in the stage 2 model by penalizing discrepancies between modelled and direct estimates. Unlike previous benchmarking approaches that adjust the point estimates only [13, 67], Bayesian benchmarking adjusts the entire posterior — automatically accounting for benchmarkinginduced uncertainty.
In this work we simultaneously enforced two benchmarks referred to as “state” and “majorbystate”. The state benchmark had seven groups which were composed of the states and territories of Australia (except the Northern Territory, which was not benchmarked due to ABS instruction [57]).
The majorbystate benchmark had 12 groups, composed of the interaction of the states and territories of Australia (except the Northern Territory) and dichotomous remoteness (major city vs nonmajor city). Thus, for each state, apart from Tasmania (where all areas were nonmajor city), and the Australian Capital Territory (where 96% of areas were major city), each SA2 was benchmarked differently depending on whether the area was in a major city or not.
External validation
External validation was performed by comparing the estimates to those from the SHAA at the PHA level and the overall trends observed in the modelled results with the general findings from other Australian health surveys conducted on specific subpopulations, such as states [68] or First Nations Australians [69]. Although this validation affirmed the validity and reliability of our estimates in general, it was particularly helpful in assessing the credibility of estimates for areas that could not be benchmarked.
Computation
We used fully Bayesian inference using MCMC via the R package rstan Version 2.26.11 [70]. Where possible we used the nonmean centered parameterization for random effects and the QR decomposition for fixed effects [71]. The stan code for the stage 1 and stage 2 models can be found on GitHub [72].
For the stage 1 model we used 1000 warmup and 1000 postwarmup draws for each of the four chains, feeding a random subset of 500 posterior draws from the stage 1 to the stage 2 model. For the stage 2 model we used 3000 warmup and 3000 postwarmup draws for each of four chains. For storage reasons we thinned the final posterior draws from the stage 2 model by 2, resulting in 6000 useable posterior draws.
Convergence of the models was assessed using trace and autocorrelation plots, effective sample size and \(\hat{R}\) [73]. While convergence ranged slightly between risk factors, all the proportion parameters, \(\varvec{\mu } = \left( \mu _1, \dots , \mu _M \right)\), had \(\hat{R} < 1.03\), with 96% having effective sample sizes \(>1000\) and 99% having \(\hat{R} < 1.01\).
Summaries and visualisation
Estimates from the benchmarked stage 2 model were reported in a variety of forms, including absolute, relative and classification measures. For point estimates we used posterior medians and for uncertainty intervals we used 95% highest posterior density intervals (HPDIs). We used the modelled proportions as the absolute indicator and odds ratios (ORs) as the relative indicator. The ORs for the tth posterior draw were derived as,
with \(\hat{\mu }^D\) being the national prevalence estimate for the risk factor measure. An OR above one indicates that the SA2 has a prevalence higher than the national average.
In addition to using point estimates and credible intervals to summarize the ORs, we also used the exceedence probability (EP) [31, 53, 74].
Generally an EP above 0.8 (or below 0.2) is considered to provide evidence that the proportion in the corresponding SA2 was substantially higher (or lower) than the national average, respectively [75]. Note that the exceedance probabilities calculated using either ORs or prevalence are identical.
To facilitate decisionmaking, we classified SA2s by assessing whether their individual and neighbor values (i.e. clusters [76, 77]) were significantly different to the national average. In this work, these classifications were called evidence classifications. Any area classified as HC, H, L, or LC has an exceedance probability suggesting that the modelled prevalence is significantly different to the national average; HC or H denotes higher, while L or LC denotes lower. The difference between HC and H (or LC and L) is that the former provides an indication of clustering of areas, while the later only indicates significance of the area itself. If an area is not classified according to the criteria above (defined as None (“N”)) the modelled estimate is not sufficiently different to the national average. See details in Section D.3 of the Additional File 1.
Code to produce subsequent plots is available on GitHub [72].
Results
Prevalence
Large spatial variation in the proportion of cancer risk factors across Australia can be clearly observed in Figs. 4,5,6 and Section H of the Additional File 1. Slightly more heterogeneity of the point estimates was observed within major cities as a result of the much greater socioeconomic variation within these areas. For example, the range of principal component 1 (a proxy for SES that is unique to the SES index) was largest in major cities and inner regional areas, but 50% the size in remote and very remote areas.
Stratifying by risk factor, the results highlight interesting patterns and trends. A more thorough discussion of the result is given in Section F of the Additional File 1.

Current smoking (Section H.2 in the Additional File 1): Spatial patterns show lower prevalence in major cities and less disadvantaged areas. Although very high prevalence estimates are observed in the very remote regions in the middle of the country, these estimates come with substantial uncertainty.

Risky alcohol consumption (Section H.3 in the Additional File 1): The spatial patterns were inconsistent with the other factors, particularly in terms of the relationship between (higher) socioeconomic status and healthy behaviours. The results suggest that less disadvantaged areas have higher prevalence, which generally manifests in higher prevalence in major cities. Unlike other risk factors where prevalence estimates exhibit relative homogeneity within the SES index deciles and remoteness groups (see Section G of the Additional File 1), for risky alcohol consumption the estimates exhibit far greater heterogeneity for more disadvantaged areas in major cities.

Inadequate diet (Section H.4 in the Additional File 1): The spatial patterns suggest less dependence on the SES index and remoteness than the other risk factors. Inadequate diet exhibits the lowest heterogeneity of the risk factors considered in this work.

Body weight (Sections H.5 to H.7 in the Additional File 1): Similar spatial patterns are observed for the three measures. The prevalence was very strongly tied to remoteness with substantially lower prevalence almost exclusively occurring in major cities. Furthermore, the most notable differences in patterns between the estimates for obese and overweight/obese are found in major cities.

Physical activity (Sections H.8 to H.9 in the Additional File 1): Similar spatial patterns are observed for the two measures. Lower prevalence of inadequate activity is observed in major cities and less disadvantaged areas.
The estimates demonstrate reliability, as around 97% of them possess coefficients of variation (CV) below 25% — a widely accepted threshold for reliability [25]. Furthermore, the modelled estimates show considerable stability improvements over the SA2 direct estimates with a reduction in variability (measured by standard deviation) across Australia by an average factor of 3.3. The estimate uncertainty varied by risk factor, with current smoking having the highest median CV (17.4) and inadequate activity (leisure) having the lowest (2.6). The distribution of CV also varied by remoteness; the median CV for major cities (62% of the survey data) was, on average, 1.8 to 3.1 times smaller than that for very remote areas.
To investigate the impact of the finer resolution, we derived PHA level CVs by taking the population weighted mean of the SA2 level estimates. The CVs of point estimates at the SA2 level range from 5% to 34% larger than point estimates at the PHA level across the risk factors. Similarly, by calculating and summarising the heterogeneity of SA2s within each PHA, we find that across the risk factors, the median PHA CV is between 1.5% to 9.1%. Of the PHAs composed of multiple SA2s, 10% have CVs \(>15\)%. The large CVs indicate that the corresponding PHAs were highly heterogeneous, highlighting the benefits of using higher resolution estimates. Given the similar definitions for the obese risk factor measure, Fig. 7 compares the estimates used in this work and that of the SHAA, indicating strong agreement.
Section G of the Additional File 1 provides more plots describing the modelled results, including how they vary by the SES index and remoteness. An interactive exploration of the modelled results will be made available in the Australian Cancer Atlas 2.0 [2], planned for release in early 2024.
Evidence classifications
Table 3 summarises the number of evidence classifications for each risk factor measure. Figure 8 stratifies these by remoteness. A similar stratified plot for the SES index is found in Section G of the Additional File 1.
Across most risk factors many more HC areas are identified than LC areas. For example, for risky waist circumference, around 783 SA2s are classed as HC, while only around 397 are classed as LC. We observed that HC or H evidence classifications are generally found in the most disadvantaged areas, while L or LC areas are more likely in the least disadvantaged areas in major cities.
The evidence classifications revealed several interesting trends. For the physical activity risk factor measures, a larger proportion of the areas in major cities were classified as HC or H as opposed to LC or L. For inadequate physical activity, the HC classifications favour the most disadvantaged areas. The weight risk factor measures exhibited different trends with a relatively even distribution of evidence classifications in major cities. Furthermore, almost all areas classified as LC or L occurred in less disadvantaged areas. As mirrored in the maps, the evidence classifications for smoking suggest a very strong correlation with remoteness and SES; almost all the LC or L classifications occur in major cities and less disadvantaged areas. Inadequate diet has the smallest number of evidence classifications (1155 out of 2221), with the largest proportion of them being LC areas in major cities and less disadvantaged areas. The results for risky alcohol consumption suggest that less disadvantaged areas have higher proportions of risky alcohol consumption; a trend unique to this risk factor measure.
Discussion
This work improves the spatial resolution and reach of previously published cancer risk factor estimates in Australia. While the estimates highlight broadly similar findings as those from the SHAA, they provide greater resolution and reach allowing for more granular exploration of the spatial disparities (see Fig. 7). This is particularly pertinent due to the heterogeneity of the component SA2s within each PHA in terms of population size, socioeconomic status and remoteness.
By improving the reach of the previously published cancer risk factor estimates, the estimates in this work uniquely enable the exploration of spatial disparities in very remote areas of Australia. As expected, the very remote areas have far greater uncertainty than those in major cities (CVs greater than 3 times higher). Nevertheless, by utilising the estimates and their uncertainty measures policy makers will have the capability to more effectively allocate health interventions and resources to these disadvantaged areas and triage areas where more data should be collected in the future to improve the quality of small area estimates.
The cancer risk factor estimates generated in this work reveal substantial spatial disparities in cancer risk behaviours across Australia, with higher prevalence of high risk behaviours generally occurring in more remote areas. While the prevalence of most cancer risk factors is higher in areas of lower SES, the spatial patterns for risky alcohol consumption demonstrated the opposite effect. Point estimates for risky alcohol consumption and current smoking exhibited the most heterogeneity across Australia, while those from the physical activity measures exhibit the least. The distribution of the point estimates are mostly consistent across states and territories of Australia.
Although generating prevalence estimates and their uncertainty intervals are useful in a variety of applications, using them to visualize which areas are substantially different to the national average can be difficult as the two components must be considered jointly. By further classifying the estimates according to their posterior probabilities, we were able to streamline this process. Classifications, such as those used in this work, are pivotal in developing targeted interventions as they enable policymakers to quickly identify areas, or groups of areas, with substantially higher (or lower) prevalence.
Our Bayesian methodology, along with its associated exceedance probabilities and evidence classifications, provides insights that cannot easily be attained via the estimates from the SHAA. Although the spatial patterns of the evidence classifications vary by risk factor, a consistent pattern was that areas with lower than average prevalence of risk factors (classified as LC or L) were almost exclusively located in major cities. Although there were areas with higher than average prevalence (HC or H) in major cities these were often less common, except for the physical activity risk factors where about half were higher and lower than the national prevalence.
Although this applied work represents a significant step in the ongoing improvements in cancer prevention in Australia, it has some limitations. Firstly and most critically, like previous research [47, 78], most of the risk factor measures used were based on data derived from selfreports which are highly susceptible to various biases [79]. Furthermore, some 2017–18 NHS questions focused on behaviour from the previous week (e.g. alcohol, physical activity), while others on a usual week (e.g. fruit and vegetables consumption, smoking).
Given the nature of the survey questions, caution must be exercised in using the risk factor measures presented herein. While the estimates provide insights into the spatial variation, due to the ecological fallacy [80] and the often varying lag time between exposure (to a risk factor) and a cancer diagnosis [3], the estimates here cannot be used to establish individuallevel associations between risk factors and cancer incidence. Moreover, as these estimates are derived from crosssectional data, they do not enable inference into lifetime risky behaviour or causal relationships with cancer.
The second limitation, relevant to any spatial analysis of lattice data, is the modifiable areal unit problem (MAUP) [81]. The MAUP refers to the sensitivity of estimates to the specified definition of a small area (e.g. choice of partitioning and resolution). While we have presented our estimates for SA2s, which offers wide applicability, we acknowledge that this particular partitioning of Australia represents just one of countless possible configurations, each yielding unique results. Thus, the conclusions drawn from our estimates are inherently entwined with the choice of partitioning and resolution of the small areas we employed [82].
Thirdly, the accuracy of our estimates are conditional on the 2017–18 NHS exclusions (very remote areas, discrete Aboriginal and Torres Strait Islander communities and nonprivate dwellings [57]). Without data for these subpopulations, there is currently no way to assess the impact of these exclusions on modelled estimates from this survey.
Next, while the SHAA provides estimates by sex [4], our study, constrained by the sparsity of the survey data at the SA2 level, did not allow for a similar disaggregation. Given the evidence that health behaviours can depend on sex, the non sexspecific estimates generated in this work may suffer from inadvertent smoothing toward the mean.
The final limitation is that the quality, in terms of both bias and variance, of small area estimates can always be improved by using larger surveys. Although we used the best survey data available, in the future, linkage of multiple surveys could provide much larger sample sizes across Australia, enabling the production of higher resolution estimates.
In terms of future research directions, one approach could involve developing distinct models for each of the eight risk factor measures. That is the linear predictor for each risk factor measure could have different sets of covariates, random effect structures or even include nonlinear relationships via splines. Alternatively, future work could model the numerous risk factors jointly by leveraging univariate stage 1 models, followed by a multivariate spatial stage 2 model [83].
Conclusions
Using a Bayesian twostage small area estimation model we have, for the first time, generated and validated point estimates of the prevalence of eight cancer risk factors, and measures of their uncertainty, at the SA2 level across Australia. By aggregating the estimates, we have shown that they are very similar to those given by the SHAA [4], external surveys [84,85,86,87,88,89] and previous research on how area level socioeconomic status and remoteness relate to healthy behaviours [42]. The new estimates provide improved spatial resolution and reach and will enable more targeted cancer prevention strategies at the small area level. Furthermore, by including the results in the next release of the Australian Cancer Atlas [2], this work promises to provide a more comprehensive picture of cancer in Australia. Since the health factors used in this study are also common risk factors for other diseases, the prevalence estimates generated here may be useful in other disease modelling applications both in Australia and internationally.
Data availability
The 2017–18 National Health Survey microdata cannot be shared publicly due to the ABS privacy policy. An application can be made to the ABS directly to gain access to the data. The statistical analysis has been conducted in the secure ABS DataLab online computing environment. The modelled estimates are available on GitHub [72].
Abbreviations
 ABS:

Australian Bureau of Statistics
 ACA:

Australian Cancer Atlas
 ACT:

Australian Capital Territory
 AIHW:

Australian Institute of Health and Welfare
 ASGS:

Australian Statistical Geography Standard
 BMI:

Body Mass Index
 DOH:

Department of Health
 CCQ:

Cancer Council Queensland
 EP:

Exceedance probability
 HPDI:

Highest posterior density interval
 IQR:

Interquartile range
 IRSD:

Index of Relative SocioEconomic Disadvantage
 MCMC:

Markov Chain monte carlo
 MrP:

Multilevel regression and poststratification
 NHMRC:

National Health and Medical Research Council
 NHS:

National Health Survey
 NSW:

New South Wales
 NT:

Northern Territory
 PC:

Principal component
 PHN:

Primary health network
 QLD:

Queensland
 SA:

South Australia
 SA2:

Statistical area level 2
 SA3:

Statistical area level 3
 SA4:

Statistical area level 4
 SEIFA:

SocioEconomic Indexes for Areas
 SES:

Socioeconomic status
 SHAA:

Social Health Atlases of Australia
 TAS:

Tasmania
 TSLN:

Twostage logisticnormal
 VIC:

Victoria
 WA:

Western Australia
References
Sung H, Ferlay J, Siegel RL, Laversanne M, Soerjomataram I, Jemal A, Bray F. Global cancer statistics 2020: Globocan estimates of incidence and mortality worldwide for 36 cancers in 185 countries. CA Cancer J Clin. 2021;71(3):209–49. https://doi.org/10.3322/caac.21660.
Cancer Council Queensland, Queensland University of Technology, FrontierSI: Australian Cancer Atlas (2023). https://gis.cancer.gov/canceratlas/app/
Whiteman DC, Webb PM, Green AC, Neale RE, Fritschi L, Bain CJ, Parkin DM, Wilson LF, Olsen CM, Nagle CM, Pandeya N, Jordan SJ, Antonsson A, Kendall BJ, Hughes MCB, Ibiebele TI, Miura K, Peters S, Carey RN. Cancers in australia in 2010 attributable to modifiable factors: introduction and overview. Australian New Zealand J Public Health. 2015;39(5):403–7. https://doi.org/10.1111/17536405.12471.
Public Health Information Development Unit: Social Health Atlases of Australia (2018). https://phidu.torrens.edu.au/socialhealthatlases
Centre for Disease Control: PLACES: Local Data for Better Health (2021). https://www.cdc.gov/places
AbdelRahman O. Disparities in modifiable cancer risk factors among canadian provinces, territories, and health regions. Curr Med Res Opin. 2021. https://doi.org/10.1080/03007995.2021.1971184.
Mansori K, SolaymaniDodaran M, MosaviJarrahi A, Motlagh AG, Salehi M, Delavari A, AsadiLari M. Spatial inequalities in the incidence of colorectal cancer and associated factors in the neighborhoods of tehran, iran: Bayesian spatial models. J Prevent Med Public Health. 2018;51:33–40. https://doi.org/10.3961/jpmph.17.167.
Samouda H, RuizCastell M, Bocquet V, Kuemmerle A, Chioti A, Dadoun F, Kandala NB, Stranges S. Geographical variation of overweight, obesity and related risk factors: findings from the european health examination survey in luxembourg, 2013–2015. PLoS ONE. 2018. https://doi.org/10.1371/journal.pone.0197021.
Australian Institute of Health and Welfare: The relationship between health risk factors and the neighbourhood environment. Report, AIHW (2022). https://doi.org/10.25816/jpx51f37
Australian Institute of Health and Welfare: Australia’s health 2018. Report, AIHW (2018). https://doi.org/10.25816/5ec1e56f25480 . https://www.aihw.gov.au/reports/australiashealth/australiashealth2018
Australian Institute of Health and Welfare: Australia’s health 2020: in brief. Report, AIHW (2020). https://www.aihw.gov.au/getmedia/2aa9f51bdbd64d568dd406a10ba7cae8/aihwaus232.pdf.aspx?inline=true
Rao JNK, Molina I. Small Area Estimation. 2nd ed. Hoboken, New Jersey: Wiley series in survey methodology; 2015.
Pfeffermann D. New important developments in small area estimation. Stat Sci. 2013;28(1):40–68.
Fay RE, Herriot RA. Estimates of income for small places: an application of jamesstein procedures to census data. J Am Stat Assoc. 1979;74(366):269–77. https://doi.org/10.2307/2286322.
Battese GE, Harter RM, Fuller WA. An errorcomponents model for prediction of county crop areas using survey and satellite data. J Am Stat Assoc. 1988;83(401):28–36. https://doi.org/10.1080/01621459.1988.10478561.
Fuglstad GA, Li ZR, Wakefield J. The two cultures for prevalence mapping: Small area estimation and spatial statistics (2021) arXiv:2110.09576
Janicki R. Properties of the beta regression model for small area estimation of proportions and application to estimation of poverty rates. Communications in Statistics  Theory and Methods 2020;49:(9), 2264–2284 https://doi.org/10.1080/03610926.2019.1570266
Mercer L, Wakefield J, Chen C, Lumley T. A comparison of spatial smoothing methods for small area estimation with sampling weights. Spatial Stat. 2014;8(1):69–85. https://doi.org/10.1016/j.spasta.2013.12.001.
Moura FAS, Migon HS. Bayesian spatial models for small area estimation of proportions. Stat Model. 2002;2(3):183–201. https://doi.org/10.1191/1471082x02st032oa.
Paige J, Fuglstad GA, Riebler A, Wakefield J. Designand modelbased approaches to smallarea estimation in a lowand middleincome country context: comparisons and recommendations. J Surv Stat Methodol. 2022;10(1):50–80.
Hogg, J., Cameron, J., Cramb, S., Baade, P., Mengersen, K.: A twostage bayesian small area estimation method for proportions. arXiv preprint arXiv:2306.11302
Liu B, Lahiri P, Kalton G. Hierarchical bayes modeling of surveyweighted small area proportions. Surv Methodol. 2014;40:1–13.
Gelman A, Little TC. Poststratification into many categories using hierarchical logistic regression. Surv Methodol. 1997;23:2713.
Leemann L, Wasserfallen F. Extending the use and prediction precision of subnational public opinion estimation. Am J Polit Sci. 2017;61(4):1003–22. https://doi.org/10.1111/ajps.12319.
Australian Bureau of Statistics. Modelled estimates for small areas based on the 2017–18 National Health Survey. Australian Bureau of Statistics: Report; 2019.
Gao PA, Wakefield J. Smoothed modelassisted small area estimation of proportions. Can J Stat. 2023. https://doi.org/10.1002/cjs.11787.
Das S, Brakel J, Boonstra HJ. Haslett S multilevel time series modelling of antenatal care coverage in Bangladesh at disaggregated administrative levels. Surv Methodol. 2022;48(2):1.
Honaker J, Plutzer E. Small area estimation with multiple overimputation. Chicago: Midwest political science association; 2011.
Zhang JL, Bryant J. Fully bayesian benchmarking of small area estimation models. J Official Stat. 2020;36(1):197–223. https://doi.org/10.2478/jos20200010.
Australian Bureau of Statistics: Australian Statistical Geography Standard (ASGS) (2011). https://www.abs.gov.au/websitedbs/d3310114.nsf/home/australian+statistical+geography+standard+(asgs)
Duncan EW, Cramb SM, Aitken JF, Mengersen KL, Baade PD. Development of the Australian cancer Atlas: spatial modelling, visualisation, and reporting of estimates. Int J Health Geogr. 2019;18(1):1–12. https://doi.org/10.1186/s1294201901859.
Public Health Information Development Unit: Population health areas: Overview (2021). https://phidu.torrens.edu.au/helpandinformation/aboutourdata/geographicalstructures/phaoverview
Australian Bureau of Statistics: National Health Survey: First Results methodology (2018). https://www.abs.gov.au/methodologies/nationalhealthsurveyfirstresultsmethodology/201718
Australian Bureau of Statistics: Microdata: National Health Survey 201718 [DataLab]. Australian Bureau of Statistics (2017)
Australian Bureau of Statistics: 4710.0  Housing and Infrastructure in Aboriginal and Torres Strait Islander Communities, Australia, 2006 (2007). https://www.abs.gov.au/ausstats/abs@.nsf/mf/4710.0
Australian Bureau of Statistics: Household and Family Projections, Australia (2023). https://www.abs.gov.au/statistics/people/population/householdandfamilyprojectionsaustralia/latestrelease#citewindow1
Das S, Baffour B, Richardson A, Cramb S, Haslett S. Daily smoking prevalence for small domains in Australia. Research Square preprint (2023) https://doi.org/10.21203/rs.3.rs2829471/v1
Australian Bureau of Statistics: ERP by SA2 (ASGS 2016), Age and Sex, 2001 Onwards (2023). https://www.abs.gov.au/about/dataservices/dataexplorer
Australian Bureau of Statistics: 2016 Census of Population and Housing, Canberra (2016). https://www.abs.gov.au/websitedbs/censushome.nsf/home/2016
Chidumwa G, Maposa I, Kowal P, Micklesfield LK, Ware LJ. Bivariate joint spatial modeling to identify shared risk patterns of hypertension and diabetes in south africa: evidence from who sage South Africa wave 2. Int J Environ Res Public Health. 2021;18(1):359.
Australian Institute of Health Welfare: Australian burden of disease study: impact and causes of illness and death in australia 2015. Report, Australian Institute of Health Welfare (2019). https://doi.org/10.25816/5ebca2a4fa7dc . https://www.aihw.gov.au/reports/burdenofdisease/burdendiseasestudyillnessdeath2015
Patterson KAE, Cleland V, Venn A, Blizzard L. Gall S a crosssectional study of geographic differences in health risk factors among young australian adults: the role of socioeconomic position. BMC Public Health. 2014. https://doi.org/10.1186/14712458141278.
Australian Bureau of Statistics: 1270.0.55.005  Australian Statistical Geography Standard (ASGS): Volume 5  Remoteness Structure, July 2016 (2016). https://www.abs.gov.au/ausstats/abs@.nsf/Latestproducts/1270.0.55.005Main%20Features15July%202016?opendocument &tabname=Summary &prodno=1270.0.55.005 &issue=July%202016 &num= &view=
Australian Bureau of Statistics: Technical Paper: SocioEconomic Indexes for Areas (SEIFA) (2016)
Cancer Australia: Lifestyle risk factors and the primary prevention of cancer (2015)
World Cancer Research Fund and American Institute for Cancer Research: Exposures, risk factors and cancer (2018). https://www.wcrf.org/dietandcancer/exposures/
Rezende LFM, Murata E, Giannichi B, Tomita LY, Wagner GA, Sanchez ZM, CelisMorales C, Ferrari G. Cancer cases and deaths attributable to lifestyle risk factors in chile. BMC Cancer. 2020;20(1):693. https://doi.org/10.1186/s12885020071874.
Cancer Council Australia: Maintain a Healthy Weight (2023). https://www.cancer.org.au/cancerinformation/causesandprevention/dietandexercise/maintainahealthyweight
Australian Institute of Health and Welfare: Risk factors to health (2017). https://www.aihw.gov.au/reports/riskfactors/riskfactorstohealth
National Health and Medical Research Council: Australian guidelines to reduce health risks from drinking alcohol. Report, National Health and Medical Research Council 2020. 9781864960716
Health N, Council MR. Australian dietary guidelines. National Health and Medical Research Council: Report; 2013.
Department of Health: Physical activity and exercise guidelines for all Australians (2014). https://www.health.gov.au/resources/collections/collectionofphysicalactivityandsedentarybehaviourguidelinesforallages
Quick H, Terloyeva D, Wu Y, Moore K, Diez oux AV. Trends in tractlevel prevalence of obesity in philadelphia by raceethnicity, space, and time. Epidemiology 2020;1:1
Cassy SR, Manda S, Marques F, Martins M Accounting for sampling weights in the analysis of spatial distributions of disease using health survey data, with an application to mapping child health in malawi and mozambique. International Journal of Environmental Research and Public Health 19(10) (2022) https://doi.org/10.3390/ijerph19106319
Parker PA, Janicki R, Holan SH Unit level modeling of survey data for small area estimation under informative sampling: A comprehensive overview with extensions. arXiv preprint arXiv:1908.10488 (2019)
Goldstein H. Multilevel statistical models. United Kingdom: John Wiley and Sons; 2011.
Australian Bureau of Statistics: 4363.0  National Health Survey: Users’ Guide, 201718 (2017)
Hajek, J.Comment on “an essay on the logical foundations of survey sampling, part one”. The Foundations of Survey Sampling (1971)
Vandendijck Y, Faes C, Kirby RS, Lawson A, Hens N. Modelbased inference for small area estimation with sampling weights. Spatial Stat. 2016;18(1):455–73. https://doi.org/10.1016/j.spasta.2016.09.004.
Binder A. On the variances of asymptotically normal estimators from complex surveys. International Statistical Review, 1983;79–292
Savitsky TD, Toth D. Bayesian estimation under informative sampling. Electron J Stat. 2016;10(1):1677–708.
Ghitza Y, Gelman A. Deep interactions with mrp: election turnout and voting patterns among small electoral subgroups. Am J Polit Sci. 2013;57(3):762–76. https://doi.org/10.1111/ajps.12004.
Riebler A, Sørbye SH, Simpson D, Rue H. An intuitive bayesian spatial model for disease mapping that accounts for scaling (2016) arXiv:1601.01180
Wolter KM. Introduction to Variance Estimation. Springer, New York, NY (2007). https://doi.org/10.1007/9780387350998
Hidiroglou M, You Y. Comparison of unit level and area level small area estimators. Surv Methodol. 2016;42(1):41–61.
Stan Development Team: Prior Choice Recommendations. GitHub repository (2023). https://github.com/standev/stan/wiki/PriorChoiceRecommendations
Bell WR, Datta GS, Ghosh M. Benchmarking small area estimators. Biometrika. 2013;100(1):189–202. https://doi.org/10.1093/biomet/ass063.
Ministry of Health: New South Wales population health surveys (2023). https://www.health.nsw.gov.au/surveys/Pages/default.aspx
Australian Bureau of Statistics: National Aboriginal and Torres Strait Islander Health Survey (2019). https://www.abs.gov.au/statistics/people/aboriginalandtorresstraitislanderpeoples/nationalaboriginalandtorresstraitislanderhealthsurvey/201819#keystatistics
Stan Development Team: Stan. https://mcstan.org (2023)
Stan Development Team: The qr reparameterization. In: Stan User’s Guide, (2022). https://mcstan.org/docs/stanusersguide/QRreparameterization.html
Hogg, J.: ACAriskfactors (2023). https://github.com/JamieHoggdepo/ACAriskfactors
Aki V, Andrew G, Daniel S, Bob C, PaulChristian B. Ranknormalization, folding, and localization: an improved \({\widehat{R}}\) for assessing convergence of mcmc (with discussion). Bayesian Anal. 2021;16(2):667–718. https://doi.org/10.1214/20BA1221.
Dong TQ, Wakefield J. Modeling and presentation of vaccination coverage estimates using data from household surveys. Vaccine. 2021;39(18):2584–94.
Richardson S, Thomson A, Best N, Elliott P. Interpreting posterior relative risk estimates in diseasemapping studies. Environ Health Perspect. 2004;112(9):1016–25. https://doi.org/10.1289/ehp.6740.
Gramatica M, Congdon P, Liverani S. Bayesian modelling for spatially misaligned health areal data: a multiple membership approach. J Royal Stat Soc Series C. 2021;70(3):645–66. https://doi.org/10.1111/rssc.12480.
Congdon, P.: Assessing persistence in spatial clustering of disease, with an application to drug related deaths in scottish neighbourhoods. Epidemiology Biostatistics and Public Health (2020)
Reijneveld SA, Verheij RA, De Bakker DH. The impact of area deprivation on differences in health: Does the choice of the geographical classification matter? J Epidemiol Commun Health. 2000;54(4):306–13. https://doi.org/10.1136/jech.54.4.306.
Zhang X, Holt JB, Yun S, Lu H, Greenlund KJ, Croft JB. Validation of multilevel regression and poststratification methodology for small area estimation of health indicators from the behavioral risk factor surveillance system. Am J Epidemiol. 2015;182(2):127–37. https://doi.org/10.1093/aje/kwv002.
Wakefield, J., Lyons, H.: Spatial aggregation and the ecological fallacy. Chapman and Hall/CRC handbooks of modern statistical methods 2010, 541–558 (2010) https://doi.org/10.1201/9781420072884c30
Openshow, S.: A million or so correlation coefficients, three experiments on the modifiable areal unit problem. Statistical applications in the spatial science, 127–144 (1979)
Roquette R, Painho M, Nunes B. Spatial epidemiology of cancer: a review of data sources, methods and risk factors. Geospatial Health. 2017;12(1):23–35. https://doi.org/10.4081/gh.2017.504.
Gelfand AE, Vounatsou P. Proper multivariate conditional autoregressive models for spatial data analysis. Biostatistics. 2003;4(1):11–5.
Ministry of Health: HealthStats NSW (2021). https://www.healthstats.nsw.gov.au/#/home
Queensland Health: About the preventive health survey and Queensland survey analytic system (2021). https://www.health.qld.gov.au/researchreports/populationhealth/preventivehealthsurveys/about
South Australia Health: South Australian Population Health Survey (2023). https://www.sahealth.sa.gov.au/wps/wcm/connect/public+content/sa+health+internet/about+us/health+statistics/sa+population+health+survey
Department of Health Tasmania. Report on the Tasmanian Population Health Survey 2019. Department of Health Tasmania: Report; 2020.
Australian Bureau of Statistics: Results from the 201819 National Aboriginal and Torres Strait Islander Health Survey (NATSIHS). Report, Australian Bureau of Statistics, (2019). https://apo.org.au/sites/default/files/resourcefiles/202006/aponid306475.pdf
Australian Institute of Health and Welfare: National drug strategy household survey 2019. Report, AIHW, Australian Government, Canberra (2020). https://doi.org/10.25816/e42pa447
Acknowledgements
We thank the Australian Bureau of Statistics (ABS) for designing and collecting the National Health Survey data and making it available for analysis in the DataLab. The views expressed in this paper are those of the authors and do not necessarily reflect the policy of QUT, CCQ or the ABS.
Funding
JH was supported by the Queensland University of Technology (QUT) Centre for Data Science and Cancer Council QLD (CCQ) Scholarship. SC receives salary and research support from a National Health and Medical Research Council Investigator Grant (#2008313).
Author information
Authors and Affiliations
Contributions
JH led the project conception, modelling, analysis, computation and writing. All other authors were equally involved in discussion, interpretation and review.
Corresponding author
Ethics declarations
Ethics approval and consent to participate
This study has received ethical approval from the Queensland University of Technology Human Research Ethics Committee (Project ID: 4609) for the project entitled “Statistical methods for small area estimation of cancer risk factors and their associations with cancer incidence”. Ethics approval was received for the inclusion of the modelled estimates from the 2017–18 NHS into the Australian Cancer Atlas (Griffith University Human Research Ethics Committee (EC00162) Ref:2018/052). Approval was received to access the secure ABS Datalab for analyses of the NHS data (Project ID: 2021033 QUT).
Competing interests
The authors declare that they have no competing interests.
Additional information
Publisher's Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Supplementary Information
Additional file 1.
Additional material containing further details of the data and model, and more plots, maps, and results.
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.
About this article
Cite this article
Hogg, J., Cameron, J., Cramb, S. et al. Mapping the prevalence of cancer risk factors at the small area level in Australia. Int J Health Geogr 22, 37 (2023). https://doi.org/10.1186/s12942023003525
Received:
Accepted:
Published:
DOI: https://doi.org/10.1186/s12942023003525