Temporal Trend and Spatial Distribution of Drug Poisoning in Semnan Province: A Retrospective Longitudinal Study

authors:

avatar Masoudeh Babakhanian ORCID 1 , avatar Khadijeh Mamashli 2 , avatar Faezeh Ansariniya 3 , avatar Somayeh Rezaie ORCID 4 , avatar Hamed Azadi 5 , avatar Masumeh Ghazanfarpour ORCID 6 , avatar Sajjad Rahimi Pordanjani ORCID 1 , 7 , *

Social Determinants of Health Research Center, Semnan University of Medical Sciences, Semnan, Iran
Faculty of Health, Golestan University of Medical Sciences, Gorgan, Iran
Velayat Hospital of Damghan, Semnan University of Medical Sciences, Semnan, Iran
School of Nursing and Midwifery, Shahroud University of Medical Science, Shahroud, Iran
Department of Anesthesiology, School of Allied Medical Sciences, Ilam University of Medical Sciences, Ilam, Iran
Nursing Research Center, Kerman University of Medical Sciences, Kerman, Iran
Assistant Professor of Epidemiology, School of Medicine, Department of Epidemiology and Biostatistics, Semnan University of Medical Sciences, Semnan, Iran

how to cite: Babakhanian M, Mamashli K, Ansariniya F, Rezaie S, Azadi H, et al. Temporal Trend and Spatial Distribution of Drug Poisoning in Semnan Province: A Retrospective Longitudinal Study. Middle East J Rehabil Health Stud. 2023;10(1):e129352. https://doi.org/10.5812/mejrh-129352.

Abstract

Background:

Substance abuse and the risks of its consequences are the major public health problems in the world and Iran.

Objectives:

The objective of the present study is determination the epidemiological status and spatial distribution of high-risk and low-risk clusters along with estimating the relative risk of drug poisoning in Semnan in 8 years.

Methods:

A descriptive-analytical and retrospective longitudinal study was conducted that examined all 3360 poisoned patients who referred to the hospitals of Semnan province from 2011 - 2018. The total number of visits was received and analyzed by year with the relevant diagnostic code in ICD-10 from hospitals of the province, and data were analyzed using Joinpoint regression analysis, global Moran’s index and Kuldorff's spatial statistic scan at a significant level of 0.05.

Results:

The number of clients varied from 219 in 2011 to 570 in 2015. The patients' mean age was 26.7 years. The average annual number of cases per one hundred thousand individuals (annual cases/100000) was 598.4. The AAPC of the incidence of drug poisoning was 9.8 and significant in the whole province, demonstrating that the incidence of poisoning increased by an average of 9.8% per year in the whole province from 2011 to 2018. Semnan County with LLR = 533.07 was the most likely low-risk cluster, and Shahroud County was the most likely high-risk cluster with LLR = 219.94.

Conclusions:

Drug policies need to be adapted to the geographical differences in the place death. After clustering regions in poisoning incidents, the necessary interventions can be concentrated for these specific geographical locations.

1. Background

Substance abuse and the risks of its consequences are the major health problems in the world and an important health problem in Iran (1). Intentional drug poisoning means the intentional consumption of more than the prescribed or advised amount of any medicine, recreational drugs, non-ingestible drug, or excess alcohol for self-harming purposes (2). The drug poisoning with different intentions is a common method of suicide attempt (3), and intentional drug poisoning or self-harm through poisoning is a common reason for patients' visits to emergency departments (4). The suicide attempt is a significant predictor of future suicide attempts, and if individuals receive psychiatric care, they are significantly less likely to develop subsequent similar behaviors (3). Substance abuse and the risks of its consequences are the major health problems in the world and an important health problem in Iran (1). Intentional drug poisoning means the intentional consumption of more than the prescribed or advised amount of any medicine, recreational drugs, non-ingestible drug, or excess alcohol for self-harming purposes (2). The drug poisoning with different intentions is a common method of suicide attempt (3), and intentional drug poisoning or self-harm through poisoning is a common reason for patients' visits to emergency departments (4). The suicide attempt is a significant predictor of future suicide attempts, and if individuals receive psychiatric care, they are significantly less likely to develop subsequent similar behaviors (3).

Addiction to drugs and alcohol is closely associated with defective socioeconomic symptoms. Social and economic factors shape health disparities among drug users (5). They affect on health indirectly by shaping individual drug-use behaviors; they directly affect on health by affecting on the availability of sources (6). Adverse socioeconomic status and low social support are interactively consistent with suicidal ideation (5). The interventions in Public health should to improving the health of drug users and address the social and economical factors that exacerbate the health consequences of illicit drug use (6).

Accordingly, the present study examined the hypothesis that differences in the economic and social environments of the counties of Semnan Province could lead to differences in people's health-related behaviors. In recent years, Semnan Province has had a cumulative incidence rate higher than the national average in deaths from substance abuse (7) The results of one study examining the process of referral of drug addicts to hospitals in Semnan indicated that the referral process had an upward trend in the last decade (8). For proper preventive planning, we need to know how drug poisoning is distributed in different periods in a place to assess the effectiveness of environmental factors in increasing or decreasing the incidence of some drug poisoning.

Addiction to drugs and alcohol is closely related to defective socioeconomic symptoms. Social and economic factors shape health disparities among drug users (5). They affect on health indirectly by shaping individual drug-use behaviors; they directly affect on health by affecting on the availability of sources (6). Adverse socioeconomic status and low social support are interactively consistent with suicidal ideation (5). The interventions in Public health should to improving the health of drug users and address the social and economical factors that exacerbate the health consequences of illicit drug use (6).

Accordingly, the present study examined the hypothesis that differences in the economic and social environments of counties of Semnan Province could lead to differences in people's health-related behaviors. In recent years, Semnan Province has had a cumulative incidence rate higher than the national average in deaths due to substance abuse (7). Since, in the last decade, Iran has a significant increase in substances than in opioid derivatives, leading to an increase in the rate of poisoning and mortality owing to substance abuse in addition to their physical and psychiatric complications (9). In a study, which examined the process of referral of drug addicts to hospitals in Semnan, the results indicated that the referral process had been upward in recent years (8). For proper preventive planning, we need to know how drug poisoning are distributed in different periods in a place to assess the effectiveness of environmental factors in increasing or decreasing the incidence of some drug poisoning.

1.1. Literature Review

As drug poisoning data include deaths from illicit and non-illicit sources, intentional and not, examining spatial variation in the factors that relate to drug poisoning should point to spatial differences in drug source or intent. Drug poisoning is hypothetically related to a range of environmental, economic and demographic factors and these relationships will be reviewed following a description of relevant geographic characteristics of the Semnan province.

To explain rates of drug poisoning in the Semnan province and the apparent paradox of rates in this province, we examined two factors (spatial and time trend) affecting drug poisoning using spatial methods. Although recent studies have linked high elevations with increased death related drug use (10, 11), intentional or unintentional Drug intoxication all of which influence drug poisoning rates.

2. Objectives

The present study was designed and conducted to identify the time trend and investigate high-risk and low-risk clusters, time trend and spatial analysis of drug poisoning in 2011 - 2017 in Semnan Province. The results may contribute to the application of public health strategies to reduce preventable deaths.

3. Methods

3.1. Design of Study

A descriptive-analytical and retrospective longitudinal study was conducted using the Exploratory Mixed Design method, which simultaneously evaluated the effects of place (multiple-group study) and time (time-trend study) (12-14). It evaluated and compared the incidence of drug poisoning in different places of Semnan over different periods.

3.2. Target Population

The target population consisted of all people in Semnan Province from 2011 to 2018, including people who visited the emergency department of the hospitals of province during the period with a diagnosis of drug poisoning. Inclusion criteria were as follows: (1) access to patient spatial information; and (2) all patients of both genders and all ages with a poisoning code based on the ICD-10 book (Appendix 1) (15) from 2011 to 2018 and being registered in the medical records unit. The population statistics of the counties of Semnan Province were also received from the health deputies of Semnan and Shahroud universities of medical sciences. All samples were included using the census method.

3.3. Location of the Study Area

Semnan Province is located at a longitude of 53º 23' and latitude of 35º 34', and its average altitude is 1130 m. The area of this province is 97,491 km2 that constitutes 5.9% of the total area of Iran. The population of the province was 702,000 in the 2016 census.

3.4. Data Collection Sources

Data related to the patients in the mentioned period were received from the medical records unit of hospitals (Velayat, Kowsar, Amiralmomenin, Emam Khomeini, Shafa, Motamedi, Bahar, Emam Hosein, Khtamol’anbia). Demographic data of Semnan Province in the study period were received from the health deputies of Semnan and Shahroud universities of medical sciences. The patients' geographic data were obtained from the latest information available in the Google Map and the reference land layer of Iranian provinces (Georeferenced).

3.5. Procedure

The incidence of the drug poisoning indicates the rate of infection and is the result of dividing the number of new cases of the drug poisoning during a certain period by the average number of people at risk of developing the drug poisoning in that period. The incidence of the drug poisoning (per 100000 people) was considered as a dependent variable and the year as an independent variable. In addition to identifying the exact time of Joinpoints, the values of AAPC-APC and slope changes of the regression line were estimated at a 95% confidence interval for the annual drug poisoning trend from 2011 to 2011 according to the location.

The Global Moran's index was used to analyze the spatial pattern and evaluate the degree of spatial autocorrelation of the drug poisoning.

To this end, first the cumulative incidence of drug poisoning in the province from 2011 to 2018 was calculated and then the zoning of the cumulative incidence in the province was drawn using the ArcGIS 10.8; Finally, the spatial autocorrelation of drug poisoning was measured using Global Moran's I. The statistics varied from +1 to -1. The closer the value of the statistic is to +1, the more it is clustered, and the closer it is to -1, the more it is dispersed, Furthermore, if its value is zero, the spatial distribution of the phenomenon is random. The null hypothesis indicated that the distribution of drug poisoning was randomly dispersed in Semnan Province (13, 16).

The Kuldorff's purely spatial scan statistic test and the Gini Optimized Cluster Collection method were used to detect the high-risk and low-risk clusters and estimate the relative risk of drug poisoning in any county of Semnan Province. This new method separates clusters on a more local scale and optimizes the radius of the circular scanning window for cluster detection, thereby demonstrating the best results in cluster detection (13, 16, 17).

3.6. Statistical Analysis

The statistical analysis consisted of three main sections: (1) descriptive analysis of epidemiological and demographic indices the patients; (2) analysis of the evaluation of the time trend of theincidence of drug poisoning; (3) analysis of the spatial pattern of the incidence and clustering ofdrug poisoning. All statistical tests were in α of two domains at a significance level of 0.05.

SPSS 26 was used to conduct the descriptive analysis. To perform the time trend analysis, the annual incidence rate of drug poisoning per hundred thousand people per year from 2011 to 2018 wasfirst measured for the entire province and calculated by place and time. Then, changes in the time trend of the incidence of drug poisoning were evaluated using the Joinpoint Regression Program version 4.9.0.0 in March 2021. This application is used to detect the time trend and points when the incidence of the drug poisoning changes significantly.

4. Results

4.1. Descriptive Section

In the present study, 3360 patients hospitalized in Semnan Province were diagnosed with drug poisoning from 2011 to 2018 after applying the inclusion and exclusion criteria. The average annual number of cases per hundred thousand people was 598.4. The mean age of the patients was 26.7 years. The total number of the patients, who visited hospitals in Semnan Province and the annual incidence rate per 100000 people show in Table 1. As seen in this table, the highest rate of patients is in 2017 and the lowest in 2011 (Table 1).

Table 1.

Number of Patients and the Annual Incidence Rate per 100000 People (2011 - 2017)

Year2011201220132014201520162017
Number of patients219254290357576494600
Annual incidence rate per 100000 people45.152.457.466.580.58098.2

Figure 1 depicts the number of the incidence of drug poisoning in the counties of in Semnan Province during 8 years. The highest incidence rate was related to Garmsar County in 2011, 2012, 2013, 2014, and 2016, and the highest incidence rate belonged to Shahroud County in 2015 and 2018, as well as the highest incidence rate was related to Damghan County in 2017.

The incidence of drug poisoning as well as counties of Semnan Province during 8 years
The incidence of drug poisoning as well as counties of Semnan Province during 8 years

4.2. Spatial Distribution Results

Figure 2 shows the spatial distribution zoning of the cumulative incidence of drug poisoning.

Spatial distribution map of the cumulative incidence of drug poisoning in Semnan province
Spatial distribution map of the cumulative incidence of drug poisoning in Semnan province

Figure 3 presents the cumulative incidence of each county of the province. As the figure shows, the highest and lowest rates of cumulative incidence were in Garmsar and Semnan, respectively.

Spatial autocorrelation

Degree of the spatial autocorrelation of drug poisoning in Semnan Province using the Global Moran's I index
Degree of the spatial autocorrelation of drug poisoning in Semnan Province using the Global Moran's I index

Figure 3 shows the result of the patient's autocorrelation and analysis of the Global Moran's I index. The value of the Global Moran's I index was equal to 0.506, and P-value = 0.74, indicating that the degree of the spatial autocorrelation of drug poisoning was non-significant (Figure 3).

4.3. Estimation of the Relative Risk and Identification of Clusters by the Gini Optimized Cluster Collection Method in Each County of Semnan Province

The relative risk value in the low-risk cluster (Semnan County) was 0.26 and significant, indicating that the risk of drug poisoning in regions inside the cluster was 0.26 times higher that in regions outside the cluster (74% less). The relative risk in the high-risk cluster (Shahroud County) was 2.07 and significant, revealing that the risk of drug poisoning in regions inside the cluster was 2.07 times higher than the regions outside the cluster (107% more). Shahroud and Garmsar were two high-risk regions for drug poisoning so that the relative risk in the regions was higher than 1 (Table 2); in other words, the people living in these regions were at higher risk of developing drug poisoning. Semnan and Damghan are low-risk and protective regions, since they have a relative risk lower than 1.

Table 2.

Geographical Coordinates of City Centers Identified as High Risk and Low Risk in All Clusters Using the Gini Optimized Cluster with Gini Optimized Cluster Collection Method

No.CityCoordinate/RadiusPopulation SizeNumber of CasesExpected CaseAnnual Cases/100000Relative RiskObserved/ExpectedLog Likelihood RatioP-Value
1Semnan35.34 N & 53.23 E2331425131397.92219.60.260.37533.0652< 0.000
2Shahroud36.4167 N & 55.0167 E1647211573987.679532.071.59219.94< 0.000
3Garmsar35.12 N & 52.20 E91659954549.591038.72.011.74150.47< 0.000
4Damghan36.168 N & 54.348 E99370491595.82493.10.80.8211.664630.001

Among low-risk clusters, Semnan with a Log-likelihood ratio (LLR) of 533 was the most important low-risk cluster, i.e. the cluster was more likely to be a low-risk cluster. The geographical coordinates of the center of the cluster were 53.23 E and 35.34 N and had 513 poisoned patients. Damghan County with a relative risk of 0.8 and LLR = 11.66 was a secondary cluster. The geographical coordinates of the center of the cluster were 36.167N and 54.348E and had 491 poisoned patients.

Moreover, Shahroud was the most important high-risk cluster with LLR = 219.94, demonstrating that the cluster was more likely to be a high-risk cluster. The geographical coordinates of the center of the cluster were 36.4167N and 55.0167E and had 1573 patients. Garmsar County with a relative risk of 2.01 and LLR = 150.47 was a low-risk secondary cluster. The geographical coordinates of the center of the cluster were 36.188N and 52.20E and had 954 poisoned patients (Table 2).

4.4. Time Distribution Results

Table 3 and Figure 4 present the results of the time trend of the incidence of drug poisoning based on the residence county. The trend of poisoning was increasing for Semnan County (Figure 4A) from 2011 to 2014, and the value of APC was 71.1 for the period; in other words, the incidence of drug poisoning increased by 71.1% per year on average in this period, and then it had a non-significant reduction (APC = -13.0). The incidence of poisoning had a significant increase in Shahroud County (Figure 4B) from 2011 to 2015 (APC = 35.9), however,the rate of poisoning decreased after 2015. The trend of poisoning also increased in Damghan County from 2011 to 2014 and increased by 8% (APC) on average, however, the incidence of poisoning increased rapidly after 2015.

Temporal analysis with join-point regression models fitted to crude incidence rate drug related poisoning in counties of Semnan Province (2011 - 2018).
Temporal analysis with join-point regression models fitted to crude incidence rate drug related poisoning in counties of Semnan Province (2011 - 2018).

Garmsar County also had an increasing and non-significant trend from 2011 to 2014 (APC = 7.8), however, the incidence of poisoning was declining after 2015.

The poisoning time trend in the whole province was increasing by 14.1% per year on average, which was statistically significant, from 2011 to 2015, but a joinpoint was found in 2015, and the incidence of poisoning decreased by an average of 0.2% annually from 2015 to 2018. The value of AAPC was 9.8 and significant from 2011 to 2018, indicating that the incidence of poisoning increased by an average of 9.8 % in the whole province per year from 2011 to 2018 (Table 3 and Figure 4).

Table 3.

More Details of the Temporal Analysis with Join-Point Regression Models Fitted to the Incidence Rate of Drug-Related Poisoning in Semnan Province (2011 - 2018)

The Graph Related with Figure 4CityTime Frame of the TrendAPCP-ValueAAPCP-Value
Point Estimate95% CIPoint Estimate95% CI
ASemnan2011 - 201471.1-16.2 - 249.40.0916.3-9.1 - 48.8< 0.01
2014 - 2018-13.0-44.6 - 36.70.4
BShahroud2011 - 201535.9 a15.8 - 59.60.00924.1 a13.7 - 35.4< 0.01
2015 - 20189.9-14.7 - 41.60.4
CDamghan2011 - 20158.0-35.3 - 80.20.614.0-13.8 - 50.7< 0.01
2015 - 201822.5-45.5 - 175.20.4
DGarmsar2011 - 20137.8-46.1 - 115.80.71.3-11.9 - 16.51
2013 - 2018-1.2-15.4 - 15.40.8
ETotal of Province2011 - 201614.1*10.8 - 17.60.0019.8*6.9 - 12.8< 0.1
2016 - 2018-0.2-12.7 - 14.10.9

5. Discussion

Knowledge about the poisoning pattern in a particular region plays an important role in identifying risk factors and early diagnosis of poisoning. Understanding the poisoning pattern also helps to design suicide prevention strategies and reduce the risks of accidental poisoning (18). Semnan Province is a region with a high risk of deaths due to substance abuse (19); hence, the present study set to identify high-risk clusters and the time trend of drug poisoning in the province. Therefore, spatial analyses at this point are valuable for identifying risk factors for informing resource allocation.

The annual number of cases per 100,000 people was 598.4 on average and in young people. The AAPC of drug poisoning incidence was 9.8 and significant in the whole province, indicating that poisoning increased by an average of 9.8% per year in the whole province from 2011 to 2018. The results of cluster identification by the Gini Optimized Cluster Collection method demonstrating that the relative risk value was 0.26 and significant in the low-risk cluster (Semnan County), showing that the risk of drug poisoning was 0.26 times higher in regions inside the cluster than regions outside the cluster (74% less). Furthermore, the relative risk value was 2.07 and significant in the high-risk cluster (Shahroud County), revealing that the risk of drug poisoning was 2.07 times higher in the regions inside the cluster than regions outside the cluster (107% more). Semnan City with LLR = 533.07 was the most likely low-risk cluster, and Shahroud was the most likely high-risk cluster with LLR = 219.94.

Semnan Province is in the transit route of drugs from Afghanistan to Tehran owing to its proximity to Khorasan Province of Iran (20). The drug trafficking route in this province has provided easy access of people to drugs and alcohol in this region. Since international evidence shows that easy access to drugs plays key role in intentional drug poisoning (4), the trend of visits of drug addicts to hospitals in Semnan Province was increasing during the 2020s (8). Additionally, its mortality rate due to substance abuse is higher than the national average (7). In case of a lack of proper planning and intervention, this trend is predictedto increase (8).

Owing to the geographical location of Semnan Province, it has the highest number of Afghan illegal asylum seekers compared to other provinces (21), leading to a high number of addicts and suicide attempts according to the large population. Furthermore, positive attitude toward addiction in families, high migration in rich strata, more than a quarter of the county's population being below the poverty line, and unfair distribution of income in families have made caused Garmsar County to have the highest cumulative incidence from drug poisoning in the province (22, 23). Regions with low-income populations are more at risk of deadly poisoning however, the mechanism of this association is unknown (24). The history of addiction in families and positive attitude toward addiction in families in Shahroud County can increase the relative risk of addiction in this county (23). This results coincides with previous research that drugusers may be use drug more available and less costly and used by peoples who misuse opioids (25).

The data indicated that most poisoned individuals were in the second decade of their lives. A similar study indicated that the highest age of suicide attempt was from 15 to 29 years in Southeast Asian countries (26). These findings correspond with the results from studies on heroin user in Seattle where it was found that younger heroin users and opioid injectors were more at risk than older heroin users (27). The major causes of this upward include increased access to opioids, injection behaviors and sharing of syringes (28).

Studies in western countries indicated that there was nosignificant relationship between gender, age group, and addiction (29). Cultural factors and regional differences play crucial roles in this field as change of cultural norms and different paths of economic development over time change this index (26).

The results indicated that there were 598.4 annual cases/100000 people on average in the province. In the United States, the figure was 232 per 100000 in 2007 (30); it was 17 per 100,000 in 2013 in Japan (31), and 115 per 100,000 in New Zealand (32). It is important to note that these statistics provide contrary definitions of the number of people poisoned by drugs and other chemicals (e.g. alcohol, carbon monoxide, and pesticides) in various studies, and the statistics are not exclusively for drug poisoning (31).

The present article reveals a more representative picture of drug poisoning in this particular region. Probably risk factors for acute poisoning in Semnan Province include regions with lower population density and a higher poverty scale, lower urban development, and young people. The interrelated nature of the published data related to illicit opioid epidemics indicates that they cannot be ignored and the people in this province are potentially at risk of the dissemination of the data on control, problematization, stigma, class and/or regional tension; however, the numbers will probably increase in the subsequent years. Irrespective of the policies that are likely to be implemented, death rates will rise and quality of life will considerably decrease for many people in the coming years. It is better to design planning for naloxone distribution that is considered worldwide as an effective antidote for intervening deaths from drug overdoses (25).

To support the current approaches to reducing the harm, there are important research gaps to be filled in future studies. Among the gaps are surveillance, ethnographic studies on drug use behaviors; epidemiologic studies on exposure, natural histories to describe transitions in routes of administration and use, and risks associated with new synthetic opioids that are illegally manufactured; evolution of treatment trajectories; opioid markets changes; and measurement of the social and economic effects of the use of opioids, particularly heroin.

Further studies are recommended to report the causes of poisoning as well as the exact status of death or recovery of individuals according to the type of drugs and using forensic reports so that the officials may know how to better plan the quality and quantity of substance available to the public in Semnan Province and to prevent more serious accidents with timely interventions.

5.1. Conclusion

Semnan Province is a region with a high risk of deaths due to substance abuse. Shahroud and Garmsar as two border towns of Semnan province had the highest incidence and relative risk of substance abuse. After clustering regions in poisoning incidents, drug policies should be designed according to geographical differences in the region. Training packages should be provided for the most at-risk populations, and differences in demographic characteristics should be considered inside and outside the province. It is necessary to upgrade preventive and therapeutic services in high -risk substance abuse areas.

5.2. Limitations

Lack of access to some social, economic, and climatic data to measure their effects on substance abuse is one of the most important limitations of the present study that we suggest that other researchers can measure their effects on the Substance Abuse by resolving this limitation.

A limitation of the study was the lack of cooperation of the province's hospitals in providing accurate statistics on the type of toxic substance and the individuals' gender. To better understand the determinants of the spatial distribution of drug poisoning, more data should be collected and analyzed at the personal level, including the history of prescribing drugs and sources of prescribed drugs (legal vs. diversion; by a physician vs. via the Internet) Since there may be confounding factors that affect the results of spatial regression. The information system and records of patients are not without errors in the hospitals; and there is no specific national registration system in Iran to determine the number of suicides or poisoning cases due to addiction. Different organizations and sources report different statistics; thus, there are many under-enumerations in this field (7).

Acknowledgements

References

  • 1.

    Soltaninejad K, Shadnia S. Lead Poisoning in Opium Abuser in Iran: A Systematic Review. Int J Prev Med. 2018;9:3. [PubMed ID: 29416839]. [PubMed Central ID: PMC5787876]. https://doi.org/10.4103/ijpvm.IJPVM_22_17.

  • 2.

    Hatcher S, Sharon C, Collins N. Epidemiology of intentional self-harm presenting to four district health boards in New Zealand over 12 months, and comparison with official data. Aust N Z J Psychiatry. 2009;43(7):659-65. [PubMed ID: 19530023]. https://doi.org/10.1080/00048670902970833.

  • 3.

    Spiller HA, Ackerman JP, Smith GA, Kistamgari S, Funk AR, McDermott MR, et al. Suicide attempts by self-poisoning in the United States among 10-25 year olds from 2000 to 2018: substances used, temporal changes and demographics. Clin Toxicol (Phila). 2020;58(7):676-87. [PubMed ID: 31587583]. https://doi.org/10.1080/15563650.2019.1665182.

  • 4.

    Kumpula EK, Lambie B, Quigley P, Nada-Raja S, Norris P. Prescribers aware: a cross-sectional study from New Zealand emergency departments on the substances used in intentional self-poisoning and their sources. J Prim Health Care. 2020;12(3):235-43. [PubMed ID: 32988445]. https://doi.org/10.1071/HC20017.

  • 5.

    Rahimian Boogar E, Alavi K, Esfahani M. [The psychological, demographical and socio-economic factors related with suicidal ideation among general population of Semnan city in 1391]. J Arak Uni Med Sci. 2014;17(1):0. Persian.

  • 6.

    Galea S, Vlahov D. Social determinants and the health of drug users: socioeconomic status, homelessness, and incarceration. Public Health Rep. 2002;117 Suppl 1(Suppl 1):S135-45. [PubMed ID: 12435837]. [PubMed Central ID: PMC1913691].

  • 7.

    Babakhanian M, Zarghami M, Alipour A, Khosravi A, Hashemi-Nazari SS, Saberi M, et al. An Estimation of Drug-Related Deaths in Iran, Using the Capture-Recapture Method (2014-2016). Addict Health. 2020;12(2):87-97. [PubMed ID: 32782731]. [PubMed Central ID: PMC7395932]. https://doi.org/10.22122/ahj.v12i2.266.

  • 8.

    Moosazadeh M, Afshari M, Rezaie S, Sahraie Z, Ghazanfarpour M, Babakhanian M. Forecasting trends in hospitalisations due to drug poisoning in Semnan, Iran up to 2021: time series modelling. Alcoholism and Drug Addiction. 2020;33(2):151-60. https://doi.org/10.5114/ain.2020.99871.

  • 9.

    Taghaddosinejad F, Arefi M, Fayaz AF, Tanhaeivash R. Determination of substance overdose in two Iranian centers: comparison between opioids and non-opioids. J Forensic Leg Med. 2013;20(3):155-7. [PubMed ID: 23472794]. https://doi.org/10.1016/j.jflm.2012.06.012.

  • 10.

    Alipour A, Babakhanian M, Zarghami M, Khosravi A, Saberi M. Years of Life Lost (YLLs) due to drug-related deaths in the Islamic Republic Of Iran (2014-2017): a temporal and geographic pattern. Cien Saude Colet. 2022;27(7):2843-54. [PubMed ID: 35730851]. https://doi.org/10.1590/1413-81232022277.23252021.

  • 11.

    Shahbazi F, Mirtorabi D, Ghadirzadeh MR, Shojaei A, Hashemi Nazari SS. Years of Life Lost (YLL) Due to Substance Abuse in Iran, in 2014-2017: Global Burden of Disease 2010 Method. Iran J Public Health. 2020;49(11):2170-8. [PubMed ID: 33708738]. [PubMed Central ID: PMC7917516]. https://doi.org/10.18502/ijph.v49i11.4735.

  • 12.

    Morgenstern H. Ecologic studies in epidemiology: concepts, principles, and methods. Annu Rev Public Health. 1995;16:61-81. [PubMed ID: 7639884]. https://doi.org/10.1146/annurev.pu.16.050195.000425.

  • 13.

    Rahimi Pordanjani S, Kavousi A, Mirbagheri B, Shahsavani A, Etemad K. Geographical Pathology of Acute Lymphoblastic Leukemia in Iran with Evaluation of Incidence Trends of This Disease Using Joinpoint Regression Analysis. Arch Iran Med. 2021;24(3):224-32. [PubMed ID: 33878881]. https://doi.org/10.34172/aim.2021.34.

  • 14.

    Rothman KJ, Greenland S, Lash TL. Modern epidemiology. Philadelphia, USA: Lippincott Williams & Wilkins; 2008.

  • 15.

    World Health Organization. International Statistical Classification of Diseases and Related Health Problems: Alphabetical index. Geneva, Switzerland: World Health Organization; 2004.

  • 16.

    Rahimi Pordanjani S, Kavousi A, Mirbagheri B, Shahsavani A, Etemad K. Temporal trend and spatial distribution of acute lymphoblastic leukemia in Iranian children during 2006-2014: a mixed ecological study. Epidemiol Health. 2020;42. e2020057. [PubMed ID: 32777882]. [PubMed Central ID: PMC7871159]. https://doi.org/10.4178/epih.e2020057.

  • 17.

    Pordanjani SR, Kavousi A, Mirbagheri B, Shahsavani A, Etemad K. Spatial analysis and geoclimatic factors associated with the incidence of acute lymphoblastic leukemia in Iran during 2006-2014: An environmental epidemiological study. Environ Res. 2021;202:111662. [PubMed ID: 34273372]. https://doi.org/10.1016/j.envres.2021.111662.

  • 18.

    Abd-Elhaleem ZAE, Al Muqhem BA. Pattern of acute poisoning in Al Majmaah region, Saudi Arabia. Am J Clin Exp Med. 2014;2(4):79-85.

  • 19.

    Shahbazi F, Mirtorabi SD, Ghadirzadeh MR, Hashemi-Nazari SS, Barzegar A. Characterizing Mortality from Substance Abuse in Iran: An Epidemiological Study during March 2014 to February 2015. Addict Health. 2017;9(3):166-74. [PubMed ID: 29657697]. [PubMed Central ID: PMC5894796].

  • 20.

    UNODC. Drug Supply Reduction: An overview of drug supply and trafficking in Iran. Dakar, Africa: United Nations Office on Drugs and Crime; 2006, [cited 2022]. Available from: https://www.unodc.org/pdf/iran/drug_crime_situation/dsr/Supply_Reduction_trends_and_trafficking.pdf.

  • 21.

    Pishgahifard Z, Azarkamand R, Jorabloo M, Yourdkhani MS. Study of socio-economic and political effects caused by the immigration of Afghans to the City of Semnan. World Appl Sci J. 2012;19(4):497-503.

  • 22.

    Azkia M. Poverty, Vulnerability and Development: Case Study of the Garmsar and Dasht-e Azadegan Villages. Critique: Critical Middle Eastern Studies. 2005;14(1):109-25. https://doi.org/10.1080/10669920500057187.

  • 23.

    Jorablo S. [Factors affecting the change of drug consumption tastes from traditional to industrial in Garmsar]. Semnan Police Knowledge Quarterly. 2020;10(3):118-50. Persian.

  • 24.

    Hester L, Shi X, Morden N. Characterizing the geographic variation and risk factors of fatal prescription opioid poisoning in New Hampshire, 2003–2007. Annals of GIS. 2012;18(2):99-108. https://doi.org/10.1080/19475683.2012.668558.

  • 25.

    Stewart K, Cao Y, Hsu MH, Artigiani E, Wish E. Geospatial Analysis of Drug Poisoning Deaths Involving Heroin in the USA, 2000-2014. J Urban Health. 2017;94(4):572-86. [PubMed ID: 28639058]. [PubMed Central ID: PMC5533669]. https://doi.org/10.1007/s11524-017-0177-7.

  • 26.

    Värnik P. Suicide in the world. Int J Environ Res Public Health. 2012;9(3):760-71. [PubMed ID: 22690161]. [PubMed Central ID: PMC3367275]. https://doi.org/10.3390/ijerph9030760.

  • 27.

    Banta-Green C, Jackson TR, Freng S, Hanrahan M, Miller G, Reid S, et al., editors. Drug Abuse Trends in the Seattle/King County Area: 2010. Epidemiologic Trends in Drug Abuse. 2011; Bethesda, Maryland. National Institute on Drug Abuse; 2011. 1912 p.

  • 28.

    Peavy KM, Banta-Green CJ, Kingston S, Hanrahan M, Merrill JO, Coffin PO. "Hooked on" prescription-type opiates prior to using heroin: results from a survey of syringe exchange clients. J Psychoactive Drugs. 2012;44(3):259-65. [PubMed ID: 23061326]. https://doi.org/10.1080/02791072.2012.704591.

  • 29.

    Maddux JF, Hoppe SK, Costello RM. Psychoactive substance use among medical students. Am J Psychiatry. 1986;143(2):187-91. [PubMed ID: 3946652]. https://doi.org/10.1176/ajp.143.2.187.

  • 30.

    Xiang Y, Zhao W, Xiang H, Smith GA. ED visits for drug-related poisoning in the United States, 2007. Am J Emerg Med. 2012;30(2):293-301. [PubMed ID: 21367556]. https://doi.org/10.1016/j.ajem.2010.11.031.

  • 31.

    Okumura Y, Sakata N, Takahashi K, Nishi D, Tachimori H. Epidemiology of overdose episodes from the period prior to hospitalization for drug poisoning until discharge in Japan: An exploratory descriptive study using a nationwide claims database. J Epidemiol. 2017;27(8):373-80. [PubMed ID: 28242045]. [PubMed Central ID: PMC5549249]. https://doi.org/10.1016/j.je.2016.08.010.

  • 32.

    Peiris-John R, Kool B, Ameratunga S. Fatalities and hospitalisations due to acute poisoning among New Zealand adults. Intern Med J. 2014;44(3):273-81. [PubMed ID: 24372720]. https://doi.org/10.1111/imj.12364.