Open Access Research article

Geographic remoteness, area-level socioeconomic disadvantage and inequalities in colorectal cancer survival in Queensland: a multilevel analysis

Peter D Baade1235*, Paramita Dasgupta1, Joanne F Aitken134 and Gavin Turrell2

Author Affiliations

1 Cancer Council Queensland, Brisbane, Australia

2 School of Public Health and Social Work, Queensland University of Technology, Brisbane, Australia

3 Griffith Health Institute, Griffith University, Gold Coast, Australia

4 School of Population Health, University of Queensland, Brisbane, Australia

5 Senior Research Fellow, Cancer Council Queensland, PO Box 201, Spring Hill, QLD 4001, Australia

For all author emails, please log on.

BMC Cancer 2013, 13:493  doi:10.1186/1471-2407-13-493

The electronic version of this article is the complete one and can be found online at: http://www.biomedcentral.com/1471-2407/13/493


Received:24 June 2013
Accepted:16 October 2013
Published:24 October 2013

© 2013 Baade et al.; licensee BioMed Central Ltd.

This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract

Background

To explore the impact of geographical remoteness and area-level socioeconomic disadvantage on colorectal cancer (CRC) survival.

Methods

Multilevel logistic regression and Markov chain Monte Carlo simulations were used to analyze geographical variations in five-year all-cause and CRC-specific survival across 478 regions in Queensland Australia for 22,727 CRC cases aged 20–84 years diagnosed from 1997–2007.

Results

Area-level disadvantage and geographic remoteness were independently associated with CRC survival. After full multivariate adjustment (both levels), patients from remote (odds Ratio [OR]: 1.24, 95%CrI: 1.07-1.42) and more disadvantaged quintiles (OR = 1.12, 1.15, 1.20, 1.23 for Quintiles 4, 3, 2 and 1 respectively) had lower CRC-specific survival than major cities and least disadvantaged areas. Similar associations were found for all-cause survival. Area disadvantage accounted for a substantial amount of the all-cause variation between areas.

Conclusions

We have demonstrated that the area-level inequalities in survival of colorectal cancer patients cannot be explained by the measured individual-level characteristics of the patients or their cancer and remain after adjusting for cancer stage. Further research is urgently needed to clarify the factors that underlie the survival differences, including the importance of geographical differences in clinical management of CRC.

Keywords:
Colorectal cancer; Epidemiology; Survival; Inequalities; Multilevel

Background

Worldwide, colorectal cancer (CRC) was the second most common invasive cancer in 2008 and the fourth most deadly form of cancer [1]. Advances in cancer prevention, screening, and management over recent decades [2] have contributed to the ongoing improvements in CRC survival in developed countries [1] with Australia having one of the highest survival rates globally [3]. However not all patients have benefited equally from these advances, with international studies consistently reporting survival inequalities by area disadvantage and heath care access, [4-6] with evidence that these inequalities may be widening [7]. Australians living outside major cities, in socioeconomically disadvantaged regions or further away from radiation facilities also have poorer survival after a diagnosis of colon or rectal cancer [8-11]. Inequities in oncology services and general health care provision with increasing geographical isolation in Australia have been well documented [12] and acknowledged to be contributing factors to the greater burden for remote cancer patients.

Nonetheless, relatively few studies have quantified the impact that area-level factors have on geographical inequalities in CRC survival while specifically considering the effect of the underlying nested geographical structure [5,13]. Multilevel models enable us to simultaneously estimate the impact of both individual- and area-level explanatory variables on the total variation in individual outcomes while accounting for the clustering of observations within the same geographical location. Improved computing capacity has led to the increasing adoption of sophisticated multilevel techniques for large-scale population-based studies to quantify geographical inequalities in cancer survival and explore underlying causes [5,13-15].

A recent study examined the extent of spatial variation in CRC relative survival across small areas in Queensland [16]. However, that study was designed primarily to measure the impact that spatial variations had on premature mortality and so utilized data aggregated over each region and combinations of covariates. This removed any opportunity to simultaneously examine the impact that area- and individual-level factors had on differences in survival between individual patients.

In this study we explore whether geographical remoteness and socioeconomic characteristics of the area where a CRC patient resides at diagnosis are associated with their survival, independently of the characteristics of the individual patients themselves. To the best of our knowledge it is the first Australian study to quantitatively assess the independent associations between the characteristics of geographical areas and the characteristics of individuals in those areas with survival.

Specifically we aimed to:

i. assess whether all-cause and CRC-specific survival varied with a patient’s area of residence while controlling for within-area variation in individual effects and between-group variation in area-level factors;

ii. explore the independent impact of remoteness and area disadvantage on survival after adjusting for individual characteristics;

iii. identify individual-level factors influencing CRC survival; and

iv. explore the effect of interactions between area-level factors on survival.

Being able to quantify geographical variations in survival and identify associations with these disparities may assist advocates and health planners to develop strategic plans and public health interventions to reduce these inequalities.

Methods

Ethical approval to conduct this study was obtained from the University of Queensland Social and Behavioral Sciences Ethical Review Committee. Queensland Health provided legislative approval to access routinely collected population-based cancer data in Queensland.

Study cohort

All incident cases of invasive CRC (ICD-O3 codes C18 to C20, C218) diagnosed between January 1, 1997 and December 31, 2007 (inclusive) were extracted from the state-wide population-based Queensland Cancer Registry to which all confirmed invasive cancers diagnosed among Queensland residents must be legally reported. Data quality is high as evidenced by the low percentage (1.4%) of death certificate notifications only and high percentage (92.1%) of histologically verified cases in 2007 [17]. We restricted our cohort to those aged between 20 and 84 years at diagnosis since CRC is relatively rare among younger age groups, while among older patients death certificates are less accurate [18] and their clinical management is different [19,20]. Cases were excluded if they were notified by death certificate only, were first identified at autopsy or could not be geocoded. For patients with multiple primary colon or rectal cancers, only the tumor with the most advanced stage was considered. Variables extracted (categories in Table 1) included year and age at diagnosis, gender, occupation, marital status, country of birth, CRC site (colon C18; rectum C19-C20,C218), differentiation and Indigenous status, with the latter being considered sufficiently complete for analysis [21].

Table 1. Cohort description and unadjusted five year estimates of all-cause and colorectal cancer-specific outcomes for colorectal cancer patients aged 20–84 in Queensland, 1997–2007

Geocoding and travel distance calculations

Residential addresses were geocoded using full street address (98.0% of cases), a street at the center of the suburb (1.8%) or the post code (0.2%) at diagnosis. Radiotherapy facilities in Queensland are concentrated in larger cities and typically affiliated to major cancer care centers; hence these distances are a proxy measure of access to optimum cancer treatment. Geographical Information System software and a street network database were used to calculate road travel distances from each patient’s geocoded location to the closest radiotherapy facility as described previously [8]. These road travel distances represent the minimum distance, since it is possible that some patients may not have chosen the closest facility for treatment.

Survival data

The study cohort was followed up to 31st December 2010. The Queensland Cancer Registry routinely all incident cases to the Registrar of Births, Deaths and Marriages and the National Death Index to ascertain mortality status for all cancer patients diagnosed in Queensland [17]. Additional data from hospitals and pathology records are used to finalize the cause of death thereby providing a high degree of accuracy; although as with all population-based registries cause of death misclassification remains a possibility. Survival was measured in years from date of diagnosis to death or the study end point. Deaths from other causes were censored when estimating CRC-specific survival. The follow up time for patients who survived more than five years after diagnosis was censored at five years.

Extraction of clinical data from pathology forms

Stage at diagnosis

As has been previously described, [8] information extracted from pathology forms [22] was used to categorize stage at diagnosis into four groups ranging from Stage I (least advanced) to IV (metastatic) based on the TNM system [23].

Surgical margins

The recorded information on the distance between the tumor and outer edge of tissue sample removed during biopsy or CRC resection was used to categorize patients as having clear (no cancer cells at outer edge of sample), positive (cancer cells present at or close to the edge) or unknown surgical margins. Cancers that are recorded with a clear margin are deemed to have been completely excised which has been shown to be associated with lower recurrence and better survival [24].

Geographical area

Statistical Local Areas (SLA, n = 478) were used as the geographical unit for this study as they are deemed to be relatively homogenous with respect to population characteristics and socioeconomic status. In 2006 there were 478 SLAs in Queensland with a median population of 5,810. Cancer incidence data across all years were mapped to the 2006 SLA boundaries based on geocoded location at diagnosis thereby removing any impact of temporal changes in geographic boundaries.

Area-level socioeconomic disadvantage

Each patient was assigned to a quintile of area disadvantage (of increasing advantage from Quintile 1) based on the Australian Bureau of Statistics-derived Index of Relative Socioeconomic Advantage and Disadvantage (IRSAD). The IRSAD is determined from census measures related to both advantage and disadvantage, such as the proportion of tertiary educated residents and the proportion of low income households in a SLA [25]. This index was chosen as it does not include Indigenous status in its derivation.

Geographic remoteness

The address of usual residence at CRC diagnosis was grouped according to their level of geographic remoteness using the Australian Standard Geographical Classification Remoteness Index [26], which is a purely geographic measure of remoteness based on road distances from population centers to various levels of service provision (see Table 1 for categories).

Statistical analysis

The five-year all-cause and CRC-specific survival rates were assessed using Kaplan-Meir analysis and estimates compared across patient-sub groups with the log rank test.

Discrete-time multilevel logistic survival models

We carried out a full discrete-time multilevel logistic survival analysis that retained the underlying nested structure. This approach differs markedly from frailty models [27] where random effects are used to model clustering effects by adjusting the standard errors to account for non-independence of data within clusters. Frailty models offer a more restricted approach than multilevel methods [28,29]. Multilevel models use measures of clustering and variance in informative ways and can simultaneously model and partition the observed variation across individual and area-levels [28,30,31]. They are increasingly the method of choice when analysing data with a clear hierarchical structure [32].

Discrete-time multilevel logistic survival models are fitted to an expanded person-period dataset, containing a sequence of binary responses for each individual from each event time (in years). This variable is coded as 1 if an individual dies during a time interval t (measured in years in current analysis) and zero otherwise. Therefore a patient who is censored is indicated by a sequence of zeros for each t while one who dies will be denoted by value 1 for year in which death occurred and zero for each previous year. Once an individual dies data collection stops for that person. Hence the discrete response for a person who died during the third year following diagnosis would be (0, 0, 1), whereas it would be (0, 0, 0) if an individual was censored that year. Discrete-time multilevel survival models are thus equivalent to fitting a logistic regression model to the expanded dataset [14,33]. The restriction of follow-up to 5 years enabled the efficient computation of these complex models, which can be problematic when using this approach for analyzing large population-based datasets with long follow up intervals [33].

Although continuous time models remain the most popular for survival analysis, discrete-time methods have several advantages, especially in the multilevel framework when using large public health data sets [33,34]. Generating an expanded person-time dataset using months or days, instead of years, would increase the size of the dataset by more than ten or 300-fold [34]. Given the size of our initial data set, this additional expansion was not feasible. It is for this reason that discrete-time methods, using years as the time variable, are preferred in the multilevel framework.

The hazard function for a discrete-time multilevel logistic survival model is the conditional probability of death in interval t given that no death has occurred in the previous intervals [33,34]. When the hazard is modelled using the logit link, the exponentiated regression parameters are interpreted as the odds ratios (OR) rather than hazard ratios. Although the baseline hazard can be modelled using dummy variables for each time interval, in practice efficient estimates of model parameters can be obtained using low order polynomials for the time [35].

When the hazard is small, which is often the case if time intervals are narrow or probability of death occurring in time interval t is low (i.e. death is a rare event), the parameter estimates from the logit and Cox models are likely to comparable. Hence discrete-time logistic regression may be considered as an approximation to the Cox model [33,34]. An additional table demonstrates this. (see Additional file 1, Additional file 2).

Additional file 1. Parameter estimates from Cox and logistic regression models for colorectal cancer patients aged 20–84 in Queensland, 1997–2007.

Format: XLSX Size: 14KB Download fileOpen Data

Additional file 2. Comparison of Parameter estimates from Cox model and a logistic regression model.

Format: PDF Size: 5KB Download file

This file can be viewed with: Adobe Acrobat ReaderOpen Data

Discrete-time multilevel logistic survival analysis was used to quantify the effect of area disadvantage and geographic remoteness on all-cause and CRC-specific survival after adjusting for individual-level factors. Models were fitted using Markov chain Monte Carlo (MCMC) [36] simulations in MLwiN version 2.26 [37] (University of Bristol, United Kingdom) interfaced with Stata (StataCorp, Texas) [38]. Convergence was assessed by trace and density plots, the autocorrelation of model parameters from posterior distributions and diagnostic tests [36] with none indicating non-convergence. After a burn in period of 40,000 iterations, parameter estimates were obtained from a further 80,000 iterations (with every 10th iteration kept). The underlying hazard was described with a second-order polynomial (i.e. time (years) and time-squared) [14].

A systematic three-step approach was used for each outcome. First we estimated null models that comprised individuals nested in SLAs without covariates. A significant area-level random term (based on the Wald χ2) [37] suggested that the modeled survival rates vary across areas in Queensland. We then added individual covariates before including area-level remoteness and neighborhood disadvantage (separately or simultaneously) with the full model being simultaneously adjusted for all explanatory variables on both levels. Interactions were tested (Wald χ2) by including both second-order terms and main-effects of scrutinized variables in the models.

Model comparison

Models were compared using the Bayesian deviance information criterion (DIC) [39] with smaller values (with a difference of at least 7 units) indicating an improvement in model fit [39].

All available covariates (Table 1) were initially used for the multivariate discrete-time multilevel logistic survival analysis. We initially ran a series of MLwiN models based on the likelihood method [37] to determine those variables that were not significant (p > 0.20) and so excluded from the final models. To explore the impact of unknown stage at diagnosis on model fit and summary measures, sensitivity analyses were carried out by repeating the all-cause and CRC-specific survival analyses under three different assumptions; (a) all unstaged cases being reclassified as Stage I, b) reclassified as Stage IV or c) equally distributed over all four stage categories.

Fixed parameter estimates are presented as odds ratios (OR) with their 95% credible intervals (CrI). Joint chi-square tests were used to assess the contribution of each variable to model fit.

The median odds ratio

The median odds ratio (MOR) [40,41] is a measure of the variation between the mortality rates of different SLAs that is not explained by the modeled risk factors. It is expressed in terms of the odds ratio scale. If the MOR is equal to 1 there is no difference between areas. Larger values indicate greater geographical variation in survival. The MOR was calculated for the discrete-time multilevel logistic survival models as:

<a onClick="popup('http://www.biomedcentral.com/1471-2407/13/493/mathml/M1','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2407/13/493/mathml/M1">View MathML</a>

where Ζ0.75 is the 75th percentile of the normal distribution and σ2 is the estimated area-level variance from the MCMC simulations. A 95% CrI for the MOR was generated from the posterior distribution of the variance [30].

The interval odds ratio

In multilevel modeling, the interpretation of an area-level risk factor such as remoteness or area disadvantage should be interpreted as the effect of the risk factor given a comparison between two SLAs of identical values of the random effect whose mortality probabilities differ only in terms of the risk factor under consideration [41]. Therefore, to interpret the area-level risk factors more generally, the unexplained between-area variability also needs to be taken into account. This is achieved using the 80% interval Odds Ratio (IOR) [30], which shows the impact of area-level risk factors on mortality when comparing SLAs with different area-level characteristics. The IOR is calculated as:

<a onClick="popup('http://www.biomedcentral.com/1471-2407/13/493/mathml/M2','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2407/13/493/mathml/M2">View MathML</a>

where β is the regression coefficient for the area-level variable, σ2 is the area-level variance and Z0.10 and Z0.90 are the 10th and 90th centiles respectively of the standard normal distribution. If the IOR does not include 1.0 it indicates that the effect of the area-level variable is large relative to the clustering effect of the SLAs.

Results

Study population

Between 1997 and 2007 there were 25,788 invasive CRC cases in Queensland. Of these 23,634 were aged 20–84 years at diagnosis who initially comprised the study cohort. The exclusion of cases that had incomplete address at diagnosis information (n = 723), were identified at autopsy (n = 33), had death certificate notification only (n = 126) or who survived for less than one day (n = 25) gave the final cohort of 22,727 cases.

Among the final cohort (Table 1), approximately 37% of cancers were diagnosed at advanced stage of which one third (31%) had metastatic (Stage IV) disease. There were 9,337 (41.1%) deaths during the first five years after diagnosis of which 7,221 were attributed to CRC.

Bivariate Kaplan-Meier survival analysis

The unadjusted 5-year all-cause and CRC-survival rates were 58.1% (95% CI: 57-58%) and 66.3% (95% CI: 66-67%) respectively (Figure 1). For both survival measures there was a difference of about 6–8 percentage points between people living in the most remote areas and those from major cities, and also between residents of the most and least disadvantaged areas. (Table 1; Figure 2) All-cause and CRC-specific survival decreased with increasing age, longer travel distances, poorer tumor differentiation or higher stage at diagnosis with poorer survival also seen for patients who were Indigenous, blue collar workers, unmarried, males or born in non-English-speaking countries.

thumbnailFigure 1. Kaplan-Meir survival curves for the cumulative probability of all-cause and colorectal cancer-specific survival five years from diagnosis for colorectal cancer patients aged 20–84 in Queensland, 1997–2007.

thumbnailFigure 2. Kaplan-Meir five-year survival curves (from diagnosis) for colorectal cancer patients aged 20–84 in Queensland, 1997–2007 by geographic remoteness (early: n = 13,155; inner regional: n = 5,139; outer regional: n = 3,485; remote: n = 948) and area socio-economic disadvantage which was categorized into 5 quintiles of increasing advantage from Quintile 1 (Quintile 1: n = 3,023; 2: n = 5,335; 3: n = 6,075; 4: 5,101; 5: 3,193). a) all-cause survival by remoteness b) colorectal cancer-specific survival by remoteness c) all-cause survival by area disadvantage d) colorectal cancer-specific survival by area disadvantage.

Discrete-time multilevel logistic survival analysis

Development of final all-cause survival model

Based on the DIC measure, model fit was markedly improved by adding the individual effects to the null model for all-cause survival (Model 2). Adding in remoteness (Model 3) or area disadvantage (Model 4) further reduced the DIC by at least 7 units. Comparing the DIC statistic of these models with the fully adjusted main-effects model (Model 5) suggested that Model 5 provided an improved fit (Table 2). The additional introduction of the area-level interaction term (Model 6) did not reduce the DIC statistic; hence we retained Model 5 as the final model for all-cause survival. Parameter estimates presented here refer to this model.

Table 2. Measures of model fit and estimates of geographical variations in all-cause and colorectal cancer-specific survival in Queensland, 1997–2007

Development of final CRC-specific survival model

The DIC statistic indicated that adjusting for individual effects (Model 8) significantly improved fit over the null model (Model 7). The DIC was further reduced by at least 7 units on introduction of remoteness (Model 9) or area disadvantage (Model 10). Based on DIC criteria model fit was further improved for the fully adjusted main-effects Model 11 (Table 2) while overall fit of the interaction model (Model 12) was poorer than its main-effects counterpart. Therefore we considered model 11 to be best fitting model for these CRC-survival data and used it for the remainder of this analysis.

Area-level interactions

Interactions between geographic remoteness and area disadvantage were also not statistically significant for all-cause (Wald χ2 = 12.22, df = 11, p = 0.347) and CRC-specific (Wald χ2 = 8.83, df = 11, p = 0.638) survival, implying that the impact of socioeconomic disadvantage on both all-cause and CRC-survival were similar for urban and rural CRC patients.

Area-level variance

The null models indicated significant evidence of geographical variation in both all-cause (Model 1; p < 0.001) and CRC-specific (Model 7; p = 0.001) survival across 478 SLAs in Queensland (Table 2). However, when successively adding the individual-level and area-level variables to the models, the amount of unexplained geographical variation decreased, to which point it became non-significant for the final model for both all-cause (Model 5, p = 0.118) and CRC-specific survival (Model 11; p = 0.231). This lack of statistical significance was reflected in the relatively low values (i.e. close to one) of the MORs in Table 2.

Impact of area-level covariates on area-level variation

All of the IOR-80 intervals (Table 3) by area-disadvantage were relatively ‘narrow’ and did not contain 1, suggesting that the impact that area-level disadvantage quintiles had on survival was large relative to the clustering effect of the SLAs. The impact of remoteness was less clear; while there was no evidence (IOR interval contained 1) that the difference between major cities and inner regions had an impact on the area-level variation; there was some difference between major city and remote areas.

Table 3. Interval odds ratios (80%) for the influence of area disadvantage or remoteness on geographical variations in all-cause and colorectal cancer survival in Queensland, 1997–2007

Fixed parameter estimate

Independent of individual-level factors, both area disadvantage (p = 0.004) and geographic remoteness (p < 0.001) were significantly associated with CRC-cancer specific survival (Table 4). Statistically significant associations were also evident between area disadvantage (p < 0.001) and remoteness (p =0.003) with all-cause survival.

Table 4. Geographic remoteness, area-disadvantage and the adjusted odds of all-cause and colorectal cancer mortality in Queensland, 1997–2007

Compared to CRC patients from the least disadvantaged quintile (Quintile 5), residents of the remaining four quintiles had worse CRC-specific (OR 1.12, 1.15, 1.20, 1.23 for Quintiles 4, 3, 2 and 1 respectively) and all-cause survival (OR ranging from 1.14 to 1.25 by quintiles of increasing disadvantage). Compared to those from major cities, living in outer regional and remote areas was also associated with significantly poorer all-cause (outer regional: OR 1.09, remote: OR 1.15) and CRC-specific survival (outer regional: OR 1.15, remote: OR 1.24).

In addition to increasing age and stage, all other individual-level clinical and socio-demographic factors (except site were independent predictors of both all-cause and CRC-specific survival in multivariate analysis. (Table 4) Finally CRC survival decreased with time (with a quadratic relationship between survival and years of follow-up). There was little difference in the parameter estimates for all the fixed effects across each set of models (full results not shown).

The sensitivity analyses for stage (full results not shown) suggested that the independent association of area-level remoteness and disadvantage with both all-cause and CRC-specific survival remained regardless of the proposed assumptions for the true distribution of cancer stage at diagnosis.

Discussion

In this large population-based study of CRC patients in Queensland we found that survival outcomes depended on where patients lived at diagnosis, and that this disparity remains after adjustment for important individual-level socio-demographic and clinical factors. Specifically our results demonstrated that residents of more disadvantaged and remote areas had significantly lower all-cause and CRC-specific survival five years after diagnosis of CRC, irrespective of their individual characteristics and irrespective of the clinical characteristics of their cancers, including cancer stage at diagnosis.

There are a number of potential explanations for the observed survival disparity including possible differences in management patterns, although without more information these remain speculative. Geographical barriers and poorer health infrastructure [12] have previously been associated with lower receipt of multimodal therapies and lower survival of CRC patients in rural or (outer) regional Australia compared to major cities, [8,42,43] and, internationally, the impact on patient outcomes of variations in hospital volume and clinical experience are well-documented [44-46] in that higher caseloads and increased specialization generally improve CRC-related outcomes. All major centers of oncological care in Queensland are located in metropolitan areas and this is likely to be part of the explanation for the relatively better outcomes in major city and inner regional areas.

However we also found strong evidence of lower survival outcomes for people living in socioeconomically disadvantaged areas, irrespective of their remoteness and other individual and clinical characteristics, meaning that distance is not the only explanation. While Australia does have universal free hospital cover, previous research has shown that colorectal cancer patients who seek medical care in private hospitals have experienced better outcomes [10]. It is also possible that these area-level effects may at least partially reflect geographical differences in the distribution of other important patient characteristics that are known to influence prognosis, such as overweight, physical inactivity, smoking, dietary patterns, comorbidities and general health status as well as treatment [4-6,9,42]. For example people living in socioeconomically disadvantaged or rural areas are more likely to engage in high risk behaviors such as smoking and decreased physical activity [47]. In addition, people living in rural and remote areas of Australia have a higher prevalence of comorbidities such as diabetes and cardiovascular conditions that can significantly impact the clinical management and prognosis for CRC patients [48]. The impact of selected individual risk factors on geographical differences in all-cause mortality for Australia has been recently quantified and suggests that interventions targeted at modifiable health factors could translate to a substantial reduction (around 36-45%) in the regional mortality differentials [49].

Since the vast majority of people diagnosed in Queensland over the current study period (1996–2007) were symptomatic [50] and the gradual implementation of the National Bowel Cancer Screening Program [51] only began in late 2006; any influence that screening had on current results would be limited. However geographical differences in participation rates may impact CRC outcomes in the future.

Important strengths of this study include the population-based coverage and high quality of incidence data from the Queensland cancer registry [17] and inclusion of all routinely available covariates in the analysis. Disease stage was clinically coded from pathological forms. All-cause survival was estimated, along with CRC-specific survival, thereby avoiding the intrinsic dependence of cause-specific survival on cause-of death coding [52]. The multilevel design takes into account the hierarchical data structure and allows the simultaneous estimation of both individual- and area-level effects on survival, something that is not possible in ecological studies. The estimated random effects from MCMC simulations were quantified using MOR [30,53] to assess the magnitude of geographic variation in a meaningful way. In addition the IOR-80 interval which integrates area-level fixed and random terms was used to quantify area-level covariate effects in comparison to the unexplained variation [30].

However the Queensland cancer registry does not collect information on potential confounders including but not limited to treatment, life style, comorbidities, family history, ongoing surveillance, stress, inflammation and other measures of individual socioeconomic status (income and education) [4,5,54-56]. Different measures of socioeconomic status are not interchangeable and can have a diverse impact on health outcomes [57]. The occupation measure used for the current study was limited in its sensitivity and precision since it was not possible to disaggregate the ‘Not stated/Inadequately described category into more homogenous groupings such as ‘home duties’, ‘retired’ or unemployed’ based on available information. Around 16% of cases could not be staged and are likely to be fairly advanced at diagnosis (based on 5-year survival estimates by stage) however sensitivity analyses confirmed that the results were consistent under the various assumptions of the missing stage information.

Conclusions

This study has demonstrated that people diagnosed with CRC in rural and disadvantaged areas have significantly poorer survival than those living in urban and affluent areas, independently of cancer stage and other individual-level characteristics. Addressing this survival disadvantage will require as a first step a commitment of resources to clarify and quantify the main causes for this disparity, and it is hoped that these results provide the necessary motivation and impetus for this to happen. The causes of these inequalities are likely to be complex and difficult to unravel, however, a better understanding is essential to inform the development of interventions to improve survival in rural and disadvantaged areas to the level of the rest of the population.

Abbreviations

CRC: Colorectal cancer; IRSAD: Index of relative socioeconomic advantage and disadvantage; SLA: Statistical local area; MCMC: Markov chain Monte Carlo; MOR: Median odds ratio; CrI: Credible interval; DIC: Deviance information criterion; OR: Odds ratio; IOR: Interval odds ratio.

Competing interests

The authors declare they have no competing interests.

Authors’ contributions

PB, GT and JA conceived the study, PD conducted the data analysis, PB and PD drafted the manuscript and PD, PD, JA and GT refined and approved the final version of the paper. All authors read and approved the final manuscript.

Acknowledgments

We thank staff in the Queensland Cancer registry for providing us with the data extracts used in this analysis.

The NHMRC is an external funding agency that provided funds to conduct this research. They had no input into the content or conclusions of this paper.

Funding

This work was supported by a research grant from the (Australian) National Health and Medical Research Council (NHMRC) (ID561700). Associate Professor Peter Baade is supported by an NHMRC Career Development Fellowship (ID1005334); Professor Gavin Turrell is supported by a NHMRC Senior Research Fellowship (ID1003710).

Financial disclosures

The authors have no financial disclosures to report.

References

  1. Ferlay J, Shin HR, Bray F, Forman D, Mathers C, Parkin DM:

    GLOBOCAN 2008: Cancer Incidence and Mortality Worldwide. 2011.

    http://globocan.iarc.fr/ webcite

    OpenURL

  2. Cunningham D, Atkin W, Lenz H-J, Lynch HT, Minsky B, Nordlinger B, Starling N: Colorectal cancer.

    Lancet 2010, 375(9719):1030-1047. PubMed Abstract | Publisher Full Text OpenURL

  3. Coleman MP, Forman D, Bryant H, Butler J, Rachet B, Maringe C, Nur U, Tracey E, Coory M, Hatcher J, McGahan CE, Turner D, Marrett L, Gjerstorff ML, Johannesen TB, Adolfsson J, Lambe M, Lawrence G, Meechan D, Morris EJ, Middleton R, Steward J, Richards MA: Cancer survival in Australia, Canada, Denmark, Norway, Sweden, and the UK, 1995–2007 (the international cancer benchmarking partnership): an analysis of population-based cancer registry data.

    Lancet 2011, 377(9760):127-138. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  4. Aarts MJ, Lemmens VE, Louwman MW, Kunst AE, Coebergh JW: Socioeconomic status and changing inequalities in colorectal cancer? A review of the associations with risk, treatment and outcome.

    Eur J Cancer 2010, 46(15):2681-2695. PubMed Abstract | Publisher Full Text OpenURL

  5. Haas JS, Brawarsky P, Iyer A, Fitzmaurice GM, Neville BA, Earle C: Association of area sociodemographic characteristics and capacity for treatment with disparities in colorectal cancer care and mortality.

    Cancer 2011, 117(18):4267-4276. PubMed Abstract | Publisher Full Text OpenURL

  6. Lejeune C, Sassi F, Ellis L, Godward S, Mak V, Day M, Rachet B: Socio-economic disparities in access to treatment and their impact on colorectal cancer survival.

    Int J Epidemiol 2010, 39(3):710-717. PubMed Abstract | Publisher Full Text OpenURL

  7. Lyratzopoulos G, Barbiere JM, Rachet B, Baum M, Thompson MR, Coleman MP: Changes over time in socioeconomic inequalities in breast and rectal cancer survival in England and Wales during a 32-year period (1973–2004): the potential role of health care.

    Ann Oncol 2011, 22(7):1661-1666. PubMed Abstract | Publisher Full Text OpenURL

  8. Baade PD, Dasgupta P, Aitken JF, Turrell G: Distance to the closest radiotherapy facility and survival after a diagnosis of rectal cancer in Queensland.

    Med J Aust 2011, 195(6):350-354. PubMed Abstract | Publisher Full Text OpenURL

  9. Kelsall HL, Baglietto L, Muller D, Haydon AM, English DR, Giles GG: The effect of socioeconomic status on survival from colorectal cancer in the Melbourne collaborative cohort study.

    Soc Sci Med 2009, 68(2):290-297. PubMed Abstract | Publisher Full Text OpenURL

  10. Morris M, Iacopetta B, Platell C: Comparing survival outcomes for patients with colorectal cancer treated in public and private hospitals.

    Med J Aust 2007, 186(6):296-300. PubMed Abstract | Publisher Full Text OpenURL

  11. Yu XQ, O’Connell DL, Gibberd RW, Armstrong BK: A population-based study from New South Wales, Australia 1996–2001: area variation in survival from colorectal cancer.

    Eur J Cancer 2005, 41(17):2715-2721. PubMed Abstract | Publisher Full Text OpenURL

  12. Underhill C, Bartel R, Goldstein D, Snodgrass H, Begbie S, Yates P, White K, Jong K, Grogan P: Mapping oncology services in regional and rural Australia.

    Aust J Rural Health 2009, 17(6):321-329. PubMed Abstract | Publisher Full Text OpenURL

  13. Lian M, Schootman M, Doubeni CA, Park Y, Major JM, Stone RAT, Laiyemo AO, Hollenbeck AR, Graubard BI, Schatzkin A: Geographic variation in colorectal cancer survival and the role of small-area socioeconomic deprivation: a multilevel survival analysis of the NIH-AARP diet and health study cohort.

    Am J Epidemiol 2011, 174(7):828-838. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  14. Dasgupta P, Baade P, Aitken J, Turrell G: Multilevel determinants of breast cancer survival: association with geographic remoteness and area-level socioeconomic disadvantage.

    Breast Cancer Res Treat 2012, 132(2):701-710. PubMed Abstract | Publisher Full Text OpenURL

  15. Dejardin O, Bouvier AM, Faivre J, Boutreux S, De Pouvourville G, Launoy G: Access to care, socioeconomic deprivation and colon cancer survival.

    Aliment Pharmacol Ther 2008, 27(10):940-949. PubMed Abstract | Publisher Full Text OpenURL

  16. Cramb SM, Mengersen KL, Turrell G, Baade PD: Spatial inequalities in colorectal and breast cancer survival: premature deaths and associated factors.

    Health Place 2012, 18(6):1412-1421. PubMed Abstract | Publisher Full Text OpenURL

  17. Queensland Cancer Registry: Cancer in Queensland: Incidence, Mortality, Survival and Prevalence,1982 to 2007. Brisbane: QCR, Cancer Council Queensland and Queensland Health; 2010. OpenURL

  18. Grulich AE, Swerdlow AJ, Silva IDS, Beral V: Is the apparent rise in cancer mortality in the elderly real? Analysis of changes in certification and coding of cause of death in England and Wales, 1970–1990.

    Int J Cancer 1995, 63(2):164-168. PubMed Abstract | Publisher Full Text OpenURL

  19. Audisio RA, Papamichael D: Treatment of colorectal cancer in older patients.

    Nat Rev Gastroenterol Hepatol 2012, 9(12):716-725. PubMed Abstract | Publisher Full Text OpenURL

  20. Patel SS, Nelson R, Sanchez J, Lee W, Uyeno L, Garcia-Aguilar J, Hurria A, Kim J: Elderly patients with colon cancer have unique tumor characteristics and poor survival.

    Cancer 2013, 119(4):739-747. PubMed Abstract | Publisher Full Text OpenURL

  21. Zhang X, Condon JR, Rumbold AR, Cunningham J, Roder DM: Estimating cancer incidence in Indigenous Australians.

    Aust N Z J Public Health 2011, 35(5):477-485. PubMed Abstract | Publisher Full Text OpenURL

  22. Krnjacki LJ, Baade PD, Lynch BM, Aitken JF: Reliability of collecting colorectal cancer stage information from pathology reports and general practitioners in Queensland.

    Aust N Z J Public Health 2008, 32(4):378-382. PubMed Abstract | Publisher Full Text OpenURL

  23. Sobin LH, Wittlekind C: TNM Classification of Malignant Tumours (6th Edition). 6th edition. New York: John Wiley and Sons; 2002. OpenURL

  24. National Institute for Health and Clinical Excellence:

    Colorectal Cancer:The Diagnosis and Management of Colorectal Cancer. NICE Clinical Guideline 131. 2013.

    http://www.nice.org.uk/nicemedia/live/13597/56957/56957.pdf webcite

    OpenURL

  25. Australian Bureau of Statistics: Information Paper: An Introduction to Socio-Economic Indexes for Areas (SEIFA), 2006. Canberra: ABS; 2008. [ABS Cat. no. 2039.0] OpenURL

  26. AIHW: Rural, Regional and Remote Health: A Guide to Remoteness Classifications. Canberra: Australian Institute of Health and Welfare AIHW Cat. No. PHE 53; 2004. OpenURL

  27. Banerjee S, Wall MM, Carlin BP: Frailty modeling for spatially correlated survival data, with application to infant mortality in Minnesota.

    Biostatistics 2003, 4(1):123-142. PubMed Abstract | Publisher Full Text OpenURL

  28. Diez Roux AV: A glossary for multilevel analysis.

    J Epidemiol Community Health 2002, 56(8):588-594. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  29. Merlo J: Multilevel analytical approaches in social epidemiology: measures of health variation compared with traditional measures of association.

    J Epidemiol Community Health 2003, 57(8):550-552. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  30. Merlo J, Chaix B, Ohlsson H, Beckman A, Johnell K, Hjerpe P, Rastam L, Larsen K: A brief conceptual tutorial of multilevel analysis in social epidemiology: using measures of clustering in multilevel logistic regression to investigate contextual phenomena.

    J Epidemiol Community Health 2006, 60(4):290-297. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  31. Merlo J, Viciana-Fernandez FJ, Ramiro-Farinas D, Res Grp L: Bringing the individual back to small-area variation studies: a multilevel analysis of all-cause mortality in Andalusia, Spain.

    Soc Sci Med 2012, 75(8):1477-1487. PubMed Abstract | Publisher Full Text OpenURL

  32. Goldstein H: Multilevel Statistical Models. 4th edition. Hoboken: Wiley; 2010. OpenURL

  33. Steele F:

    Event History Analysis: A National Centre for Research Methods Briefing Paper. 2013.

    http://eprints.ncrm.ac.uk/88/1/MethodsReviewPaperNCRM-004.pdf webcite

    OpenURL

  34. Stewart CH: Multilevel Modelling of Event History Data: Comparing Methods Appropriate for Large Datasets. PhD Thesis. University of Glasgow; 2010.

    http://theses.gla.ac.uk/2007/ webcite

    OpenURL

  35. Goldstein H: Multilevel Statistical Models. 3rd edition. London: Arnold; 2003. OpenURL

  36. Browne WJ: MCMC Estimation in MLwiN, v2.26. University of Bristol: Centre for Multilevel Modelling; 2012. OpenURL

  37. Rasbash J, Steele F, Browne WJ, Goldstein H: A Users Guide to MLwiN, v2.26. University of Bristol: Centre for Multilevel Modelling; 2012. OpenURL

  38. Leckie G, Charlton C: Runmlwin- a program to Run the MLwiN multilevel modelling software from within stata.

    J Stat Softw 2013, 52(11):1-40. PubMed Abstract | PubMed Central Full Text OpenURL

  39. Spiegelhalter DJ, Best NG, Carlin BP, Van Der Linde A: Bayesian measures of model complexity and fit.

    J R Stat Soc Ser B Stat Methodol 2002, 64(4):583-639. Publisher Full Text OpenURL

  40. Larsen K, Merlo J: Appropriate assessment of neighborhood effects on individual health: integrating random and fixed effects in multilevel logistic regression.

    Am J Epidemiol 2005, 161(1):81-88. PubMed Abstract | Publisher Full Text OpenURL

  41. Sanagou M, Wolfe R, Forbes A, Reid CM: Hospital-level associations with 30-day patient mortality after cardiac surgery: a tutorial on the application and interpretation of marginal and multilevel logistic regression.

    BMC Med Res Methodol 2012, 12:28. PubMed Abstract | BioMed Central Full Text | PubMed Central Full Text OpenURL

  42. Morris M, Platell C, Fritschi L, Iacopetta B: Failure to complete adjuvant chemotherapy is associated with adverse survival in stage III colon cancer patients.

    Br J Cancer 2007, 96(5):701-707. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  43. Young JM, Leong DC, Armstrong K, O’Connell D, Armstrong BK, Spigelman AD, Ackland S, Chapuis P, Kneebone AB, Solomon MJ: Concordance with national guidelines for colorectal cancer care in New South Wales: a population-based patterns of care study.

    Med J Aust 2007, 186(6):292-295. PubMed Abstract | Publisher Full Text OpenURL

  44. Borowski DW, Bradburn DM, Mills SJ, Bharathan B, Wilson RG, Ratcliffe AA, Kelly SB: Volume–outcome analysis of colorectal cancer-related outcomes.

    Br J Surg 2010, 97(9):1416-1430. PubMed Abstract | Publisher Full Text OpenURL

  45. Archampong D, Borowski D, Wille-Jorgensen P, Iversen LH: Workload and surgeon’s specialty for outcome after colorectal cancer surgery.

    Cochrane Database Syst Rev 2012., 3

    CD005391

    OpenURL

  46. Sacerdote C, Baldi I, Bertetto O, Dicuonzo D, Farina E, Pagano E, Rosato R, Senore C, Merletti F, Ciccone G: Hospital factors and patient characteristics in the treatment of colorectal cancer: a population based study.

    BMC Public Health 2012, 12:775. PubMed Abstract | BioMed Central Full Text | PubMed Central Full Text OpenURL

  47. Australian Institute of Health and Welfare: Australia’s Health 2012. Canberra: AIHW; 2012. [Australias Health series no.13] OpenURL

  48. Lee L, Cheung WY, Atkinson E, Krzyzanowska MK: Impact of comorbidity on chemotherapy use and outcomes in solid tumors: a systematic review.

    J Clin Oncol 2011, 29(1):106-117. PubMed Abstract | Publisher Full Text OpenURL

  49. Stevenson CE, Mannan H, Peeters A, Walls H, Magliano DJ, Shaw JE, McNeil JJ: The effect of modifiable risk factors on geographic mortality differentials: a modelling study.

    BMC Public Health 2012, 12:79. PubMed Abstract | BioMed Central Full Text | PubMed Central Full Text OpenURL

  50. Lynch BM, Baade P, Fritschi L, Leggett B, Owen N, Pakenham K, Newman B, Aitken JF: Modes of presentation and pathways to diagnosis of colorectal cancer in Queensland.

    Med J Aust 2007, 186(6):288-291. PubMed Abstract | Publisher Full Text OpenURL

  51. Australian Institute of Health and Welfare, Australian Government Department of Health and Ageing: National Bowel Cancer Screening Program: Annual Monitoring Report. Canberra: AIHW; 2009. [Cancer series no. 49. Cat. no. CAN 45] OpenURL

  52. German RR, Fink AK, Heron M, Stewart SL, Johnson CJ, Finch JL, Yin D: The accuracy of cancer mortality statistics based on death certificates in the United States.

    Cancer Epidemiol 2011, 35(2):126-131. PubMed Abstract | Publisher Full Text OpenURL

  53. Chaix B, Rosvall M, Lynch J, Merlo J: Disentangling contextual effects on cause-specific mortality in a longitudinal 23-year follow-up study: impact of population density or socioeconomic environment?

    Int J Epidemiol 2006, 35(3):633-643. PubMed Abstract | Publisher Full Text OpenURL

  54. Bass AJ, Meyerhardt JA, Chan JA, Giovannucci EL, Fuchs CS: Family history and survival after colorectal cancer diagnosis.

    Cancer 2008, 112(6):1222-1229. PubMed Abstract | Publisher Full Text OpenURL

  55. Frederiksen BL, Osler M, Harling H, Ladelund S, Jorgensen T: Do patient characteristics, disease, or treatment explain social inequality in survival from colorectal cancer?

    Soc Sci Med 2009, 69(7):1107-1115. PubMed Abstract | Publisher Full Text OpenURL

  56. Vrieling A, Kampman E: The role of body mass index, physical activity, and diet in colorectal cancer recurrence and survival: a review of the literature.

    Am J Clin Nutr 2010, 92(3):471-490. PubMed Abstract | Publisher Full Text OpenURL

  57. Braveman PA, Cubbin C, Egerter S, Chideya S, Marchi KS, Metzler M, Posner S: Socioeconomic status in health research: one size does not fit all.

    JAMA 2005, 294(22):2879-2888. PubMed Abstract | Publisher Full Text OpenURL

Pre-publication history

The pre-publication history for this paper can be accessed here:

http://www.biomedcentral.com/1471-2407/13/493/prepub