|Year : 2018 | Volume
| Issue : 1 | Page : 104
Mapping of stomach, colorectal, and bladder cancers in Iran, 2004–2009: Applying Bayesian polytomous logit model
Marzieh Nasrazadani1, Mohammad Reza Maracy2, Emanuela Dreassi3, Behzad Mahaki2
1 Student Research Center, School of Health, Isfahan University of Medical Sciences, Isfahan, Iran
2 Department of Biostatistics, School of Health, Kermanshah University of Medical Sciences, Kermanshah, Iran
3 Department of Statistics G. Parenti, University of Florence, Florence, Italy
|Date of Submission||16-Jan-2017|
|Date of Acceptance||09-Jan-2018|
|Date of Web Publication||05-Dec-2018|
Department of Biostatistics, School of Public Health, Isfahan University of Medical Sciences, Hezarjerib Street, Isfahan
Source of Support: None, Conflict of Interest: None
Background: According to the last report of Iran cancer registry, stomach, colorectal, and bladder cancers are the most prevalent cancers. The present study focused on separating the latent risk surface into shared and disease-specific components. Methods: In this study, data consisting of stomach, colorectal, and bladder cancers in 30 provinces of Iran during 2004–2009 are considered. These data are analyzed by polytomous logit model. The incidence of stomach cancer acts as the reference category (the surrogate for smoking). Then, the log odds are decomposed into shared and specific structured spatial and unstructured spatial components. These latent components help to detect spatial patterns of shared and disease-specific risk factors. Results: Central, Southern, Eastern, and Southwestern provinces are supposed as high-risk regions for shared risk factor for colorectal and bladder cancers. This shared risk factor is slightly associated more with bladder than with colorectal cancer. Northern, northwestern, and central regions and also three borderline provinces in southwestern are high-risk regions for colorectal cancer. Central, eastern, southern, and western strip of the country except Ilam are found as the high-risk regions of bladder cancer. Conclusions: After considering known shared risk factor of the three cancers, it turns out that colorectal and bladder cancers have unknown shared risk factor. The significant difference in their lifestyle and eating habits could be an assumption of the risk factor.
Keywords: Colorectal neoplasms, disease mapping, polytomous logit model, shared component model, stomach neoplasms, urinary bladder neoplasms
|How to cite this article:|
Nasrazadani M, Maracy MR, Dreassi E, Mahaki B. Mapping of stomach, colorectal, and bladder cancers in Iran, 2004–2009: Applying Bayesian polytomous logit model. Int J Prev Med 2018;9:104
|How to cite this URL:|
Nasrazadani M, Maracy MR, Dreassi E, Mahaki B. Mapping of stomach, colorectal, and bladder cancers in Iran, 2004–2009: Applying Bayesian polytomous logit model. Int J Prev Med [serial online] 2018 [cited 2021 Mar 3];9:104. Available from: https://www.ijpvmjournal.net/text.asp?2018/9/1/104/246868
| Introduction|| |
Cancer is the second leading cause of deaths worldwide so that there were about 15 million incident cancer cases and over 8 million cancer deaths worldwide in 2013.,
In Iran, cancer is the third cause of death. According to the last report of Iran cancer registry, stomach, colorectal, and bladder cancers are the most prevalent cancers.
Estimating geographical distribution of cancer can lead to recognize environmental risk factors, determination of high-risk regions, and also obtain more knowledge about the process of developing cancer in the world or regional levels. For achieving precise estimation of disease distribution in recent years, different mapping methods have been used by the researchers.,
Disease mapping is a series of statistical methods in which their aim is to gain precise estimation of incidence, prevalence, survival, or mortality rate of disease and regulate them in the form of geographical maps. Most of the studies in the field of disease mapping have focused on the modeling of one disease-relative risk.
However, joint modeling of disease is preferred more than other univariate models due to its simple interpretation, considerable improvement in the deviance information criteria, and increase in estimations' preciseness. By considering the present correlation between the regions and disease, more desirable patterns of disease changes in different regions and the relation of the diseases with each other are found as well. In recent years, a number of studies have been conducted about mapping cancer in Iran among which few ones have been focused on mapping with emphasis on the role of control risk factors.,,,,,,,
The main aims of the study are to model jointly the incidence rates of stomach, colorectal, and bladder cancers in provinces of Iran during 2004–2009 and separate the latent risk surface into shared and disease-specific components. By applying polytomous logit (PL) model, we separate common and uncommon spatial patterns of colorectal and bladder cancers; where stomach cancer is considered as the reference category and also the surrogate for shared risk factor. We also estimate the relative weight of the shared component for each disease.
Knowledge of the spatial variation of the shared and specific risk factors can help politicians to evaluate healthcare system performance; better allocation of resources and implementation of prevention policies in sensitive regions.
| Methods|| |
This is an ecological study, based on data including stomach (ICD10 code C16), colorectal (C18–C20, C26), and bladder (C67) cancers, in 30 provinces of Iran during 2004–2009. The data have been collected and made available by the Iranian Ministry of Health and Medical Education.,,,,,
To joint modeling of spatial variation related to stomach, colorectal, and bladder cancers, PL model is used with Bayesian approach. It is supposed that is nominal responsive variable including the observed count of stomach, colorectal, and bladder cancers cases in 30 provinces and follows a multinomial distribution with parameters and probability vector πik = (πi1,...,πiK), where and .
yik ~ Multinomial (ni, πik = [πil,...,πiK]') k = 1,...,3
i = 1,...,30
Thus, the types of colorectal and bladder cancers are compared with stomach cancer as a reference category; each log odd is modeled as follows:
log (φik) = αk + uik + vik
φik = πik/πiK
Where αk is disease-specific intercept of every disease and uik and vik are structured and unstructured spatial terms, respectively. The amount of parameters was zero for reference category (stomach cancer) (α3 = vi3 = ui3 = 0). Considering stomach cancer as the reference, the effect of the shared risk factors on spatial changes is adjusted.
Flat noninformative distribution is specified to disease-specific intercept (αk). Clustering component uik is modeled by use of Gaussian Markov random fields (GMRFs) to cover spatial structure. The heterogeneity terms vik are independent, each vik being (λvk represents the precision parameter).
By use of shared component model, spatially structured effect is analyzed into shared and specific effects so that every spatially component includes a coefficient to assess the relative contribution of this component to each related disease.
ui1 = usi · δ1 + upi1
ui1 = usi · δ2 + upi2
The usi term presents the shared clustering component and upi1 and upi2 present colorectal and bladder cancers' specific clustering component, respectively. Both usi and upik are modeled by use of GMRFs.
Terms log δ1 and log δ2, constrained to, are assumed to be multivariate normal distributed with zero mean and variance-covariance matrix , .
Model fitting to the data was done by use of Bayesian estimation in OpenBUGS version 3.2.1 (rev 781) software. The estimation of posterior distribution was made by Monte Carlo Markov Chain for the model parameters.
| Results|| |
In [Figure 1], the map of estimating exp (usi) related to the shared structured component is shown from 2004 up to 2009 for colorectal and bladder cancers. Since this surrogate of risk factor is distributed in all the provinces, the country can be divided into high-risk and low-risk regions. Northern regions including northeastern and western provinces and also those ones at the Caspian sea's margin are known as low-risk regions. Central, southern, southwestern, and eastern provinces are supposed as high risky regions. Since the shared risk factor is adjusted, it can be correlated with other kinds of risk factors which affect the prevalence of these two diseases rather than shared risk factors in the northern regions of the country [Figure 1].
|Figure 1: Map of estimating exp (USi) representing the surrogate component of shared risk factor|
Click here to view
The posterior medians for scale parameters associated with shared component are δ1= 0.954 (0.79–1.14) and δ2= 1.047 (0.87–1.25) for colorectal and bladder cancers, respectively; δ1/δ2 which indicates that shared risk factor is slightly associated more with bladder than with colorectal cancer.
In [Figure 2], the distribution of spatially specific structured component for colorectal and bladder cancers (exp [upi1] and exp [upi1]), in all provinces, are presented. Two observed patterns are totally different due to adjusting the shared risk factor effect of two cancers. The regions with the amount more than 1 are called as the regions exposed to risk of these factors.
|Figure 2: Map of estimating exp (UPi1) and exp (UPi2) representing the surrogate of (a) colorectal and (b) bladder cancers' specific risk factors|
Click here to view
North, North-Western, and central regions of the country and also three borderline provinces were found to have the highest risk for colorectal cancer in the whole period of time. These regions are so affected by hypoactivity risk factors.
Central, eastern, southern, and western strip of the country expect Ilam are found as high-risk regions of bladder cancer. The severity of this factor worsens by getting far from the northern regions of the country [Figure 2].
| Discussion|| |
PL model was used to joint modeling of stomach, colorectal, and bladder cancers in provinces of Iran during 2004–2009 and was separated the latent risk surface into a shared and a disease-specific components.
Several studies suggest that smoking is a major known risk factor for three mentioned diseases.,,,,
The findings show the possibility of presence of an unknown risk factor, mostly, in central and southern regions of the country. The significant difference in their lifestyle and eating habits could be an assumption. However, recognizing this risk factor needs to conduct epidemiological studies in this field. The relative effect of this risk factor is more for bladder cancers rather than colorectal cancer.
Based on the previous studies, two risk factors, including, hypoactivity,,, and occupational exposure,, are considered as the specific risk of colorectal and bladder cancers, respectively.
The provinces at the margin of Caspian Sea (Golestan, Gilan, and Mazandaran), Tehran, Qazvin, and Qom are called high-risk specific risk factor for colorectal cancer. These regions are so affected by hypoactivity risk factors. The same results were found in the studies conducted in Iran previously.,
Sports development and exercise interventions are one of the most useful ways to fight colorectal cancer. The studies show that the people being active physically are caught by this disease, 20%–30% less than the others., Unfortunately, 69% of Iranian mature people do not have a daily regular physical activity.
Central, eastern, southern regions, and western strip of the country except for Ilam province are known as high-risk regions for bladder cancer and are affected by occupational exposure risk factor during the study period. Boushehr, Fars, Kerman, and South Khorasan are known as the highest risk regions.
Numerous studies show that 20% of bladder cancer are caused by occupational exposure to bladder carcinogens that occur in industrial regions processing paint, dye, metal, pharmaceutical, and petroleum products., In some developing countries, by controlling occupational exposures to bladder carcinogens, cases of bladder cancer have decreased.,
The PL model provides the possibility of analyzing the rates of some diseases jointly without having information of population at risk. Since the interpretation, in this model, is done according to log odd, the possibility of comparing and determining the priorities is provided to implement the interventions and specify the resource.
By changing the baseline category, it is probable to recognize other shared and specific risk factors for these series of cancers.
It is suggested to perform this model on smaller regions such as counties to access more precise estimations of the shared and specific effects.
It is also suggested to consider the time dimension for a larger dataset in order to recognize the temporal trend, in addition to the spatial pattern of the diseases or to account the possible difference on latency time between the first exposure to risk factors and incidence of these diseases.,
By applying PL model, without having information of risk factor and population at risk, the surrogate of shared and specific risk factors are estimated and are regulated them in the form of geographical maps. Based on the previous studies, two risk factors, including, hypoactivity and occupational exposure are considered as the specific risk of colorectal and bladder cancers, respectively. The possibility of presence of an unknown shared risk factor of colorectal and bladder cancers is reported. The significant difference in their lifestyle and eating habits could be an assumption. However, recognizing this risk factor needs to conduct epidemiological studies in this field.
Knowledge of the spatial variation of the shared and specific risk factors can help politicians to evaluate healthcare system performance; better allocation of resources and implementation of prevention policies in sensitive regions have been detected.
Since this is an ecological study and accompanies with some biases, the result is not used in personal level, but they can be followed as an appropriate goal to find the causes.
| Conclusions|| |
After considering known shared risk factor of the three cancers, it turns out that colorectal and bladder cancers have unknown shared risk factor. The significant difference in their lifestyle and eating habits could be an assumption of the risk factor.
The authors are grateful to express their sincere thanks to the office of Noncommunicable disease, Ministry of Health and Medical Education, Iran. This work was partially supported by the student Research Center, school of health, Isfahan University of Medical Sciences, Isfahan, Iran (code: 3941060). This article has been extracted from MSc thesis in Biostatistics at School of health at Isfahan University of Medical Sciences.
Financial support and sponsorship
Conflicts of interest
There are no conflicts of interest.
| References|| |
Lozano R, Naghavi M, Foreman K, Lim S, Shibuya K, Aboyans V, et al.
Global and regional mortality from 235 causes of death for 20 age groups in 1990 and 2010: A systematic analysis for the global burden of disease study 2010. Lancet 2012;380:2095-128.
Naghavi M, Wang H, Lozano R, Davis A, Liang X, Zhou M, et al
. Global, regional, and national age-sex specific all-cause and cause-specific mortality for 240 causes of death, 1990-2013: A systematic analysis for the Global Burden of Disease Study 2013. Lancet 2015;385:117-71.
Rastaghi R, Jafari-Koshki T, Mahaki B, Bashiri Y, Mehrabani K, Soleimani A. Trends and risk factors of Gastric Cancer in Iran (2005-2010). Int J Prev Med 2018. In press.
Lawson AB, Browne WJ, Rodeiro CL. Disease Mapping With WinBUGS and MLwiN. Hoboken, NJ, USA: John Wiley & Sons; 2003.
Khoshkar AH, Koshki TJ, Mahaki B. Comparison of Bayesian spatial ecological regression models for investigating the incidence of breast cancer in Iran, 2005. Asian Pac J Cancer Prev 2015;16:5669-73.
Ahmadipanahmehrabadi V, Hassanzadeh A, Mahaki B. Bivariate spatiotemporal shared component modeling: Mapping of relative death risk due to colorectal and stomach cancers in Iran provinces. Int J Prev Med 2018. In press.
Asmarian NS, Ruzitalab A, Amir K, Masoud S, Mahaki B. Area-to-area poisson kriging analysis of mapping of county-level esophageal cancer incidence rates in Iran. Asian Pac J Cancer Prev 2013;14:11-3.
Haddad-Khoshkar A, Jafari-Koshki T, Mahaki B. Investigating the incidence of prostate cancer in Iran 2005–2008 using bayesian spatial ecological regression models. Asian Pac J Cancer Prev 2015;16:5917-21.
Jafari-Koshki T, Schmid VJ, Mahaki B. Trends of breast cancer incidence in Iran during 2004-2008: A Bayesian space-time model. Asian Pac J Cancer Prev 2014;15:1557-61.
Kavousi A, Meshkani MR, Mohammadzadeh M. Spatial analysis of relative risk of lip cancer in Iran: A Bayesian approach. Environmetrics 2009;20:347-59.
Mahaki B, Mehrabi Y, Kavousi A, Akbari ME, Waldhoer T, Schmid VJ, et al.
Multivariate disease mapping of seven prevalent cancers in Iran using a shared component model. Asian Pac J Cancer Prev 2011;12:2353-8.
Rastaghi S, Jafari-Koshki T, Mahaki B. Application of bayesian multilevel space-time models to study relative risk of esophageal cancer in Iran 2005-2007 at a county level. Asian Pac J Cancer Prev 2015;16:5787-92.
Zayeri F, Kavousi A, Najafimehr H. Spatial analysis of relative risks for skin cancer morbidity and mortality in Iran, 2008-2010. Asian Pac J Cancer Prev 2015;16:5225-31.
Dreassi E. Polytomous disease mapping to detect uncommon risk factors for related diseases. Biom J 2007;49:520-9.
Report ICR. Provincial Report, Islamic Republic of Iran-2004. Ministry of Health and Medical Education, Deputy for Health Center for Diseases Control & Management Tehran; 2005.
Report ICR. Provincial Report, Islamic Republic of Iran-2005. Ministry of Health and Medical Education, Deputy for Health Center for Diseases Control & Management Tehran; 2006.
Report ICR. Provincial Report, Islamic Republic of Iran-2006. Ministry of Health and Medical Education, Deputy for Health Center for Diseases Control & Management Tehran; 2007.
Report ICR. Provincial Report, Islamic Republic of Iran-2007. Ministry of Health and Medical Education, Deputy for Health Center for Diseases Control & Management Tehran; 2009.
Report ICR. Provincial Report, Islamic Republic of Iran-2008. Ministry of Health and Medical Education, Deputy for Health Center for Diseases Control & Management Tehran; 2010.
Report ICR. Provincial Report, Islamic Republic of Iran-2009. Ministry of Health and Medical Education, Deputy for Health Center for Diseases Control & Management Tehran; 2012.
Kelsall JE, Wakefield JC. Discussion of 'Bayesian models for spatially correlated disease and exposure data', by Best et al
. Bayesian statistics. 1999;6:151.
Held L, Natário I, Fenton SE, Rue H, Becker N. Towards joint disease mapping. Stat Methods Med Res 2005;14:61-82.
Knorr-Held L, Best NG. A shared component model for detecting joint and selective clustering of two diseases. J R Stat Soc Ser A 2001;164:73-85.
Bishehsari F, Mahdavinia M, Vacca M, Malekzadeh R, Mariani-Costantini R. Epidemiological transition of colorectal cancer in developing countries: Environmental factors, molecular pathways, and opportunities for prevention. World J Gastroenterol 2014;20:6055-72.
Downing A, Forman D, Gilthorpe MS, Edwards KL, Manda SO. Joint disease mapping using six cancers in the Yorkshire region of England. Int J Health Geogr 2008;7:41.
Ebadi N, Jahed M, Mivehchi M, Majidizadeh T, Asgary M, Hosseini SA, et al.
Interleukin-12 and interleukin-6 gene polymorphisms and risk of bladder cancer in the Iranian population. Asian Pac J Cancer Prev 2014;15:7869-73.
Jing C, Huang ZJ, Duan YQ, Wang PH, Zhang R, Luo KS, et al.
Glulathione-S-transferases gene polymorphism in prediction of gastric cancer risk by smoking and Helicobacter pylori
infection status. Asian Pac J Cancer Prev 2012;13:3325-8.
Torre LA, Bray F, Siegel RL, Ferlay J, Lortet-Tieulent J, Jemal A, et al.
Global cancer statistics, 2012. CA Cancer J Clin 2015;65:87-108.
Khoubi J, Pourabdian S, Mohebbi I, Tajvidi M, Zaroorian O, Giahi O, et al.
Association between the high risk occupations and bladder cancer in Iran: A case-control study. Int J Occup Med Environ Health 2013;26:205-13.
Olfert SM, Felknor SA, Delclos GL. An updated review of the literature: Risk factors for bladder cancer with focus on occupational exposures. South Med J 2006;99:1256-63.
Chamanparaa P, Moghimbeigi A, Faradmal J, Poorolajal J. Exploring the spatial patterns of three prevalent cancer latent risk factors in Iran; using a shared component model. Int J Epidemiol Res 2015;2:68-77.
Safaee A, Fatemi SR, Ashtari S, Vahedi M, Moghimi-Dehkordi B, Zali MR, et al.
Four years incidence rate of colorectal cancer in Iran: A survey of national cancer registry data-implications for screening. Asian Pac J Cancer Prev 2012;13:2695-8.
Chan AT, Giovannucci EL. Primary prevention of colorectal cancer. Gastroenterology 2010;138:2029-43.e10.
Alikhani S, Delavari A, Alaedini F, Kelishadi R, Rohbani S, Safaei A, et al.
Aprovince-based surveillance system for the risk factors of non-communicable diseases: A prototype for integration of risk factor surveillance into primary healthcare systems of developing countries. Public Health 2009;123:358-64.
Bassi P, Pagano F. Invasive Bladder Cancer. New York, United States: Springer Science & Business Media; 2007.
Jafari-Koshki T, Arsang-Jang S, Mahaki B. Bladder cancer in Iran: Geographical distribution and risk factors. Int J Cancer Manag 2017;10:e5610.
[Figure 1], [Figure 2]