Chiang Mai Med J 2010;49(4):161-166.

Original article

AIDS Incidence Rates in Province,

Wattanavadee Sriwattanapongse,M.S., Sukon Prasitwattanaseree, D.E.A., and Surin Khanabsakdi, M.S. Department of Statistics, Faculty of Science, Chiang Mai University

Abstract The objective of this study was to investigate the epidemic patterns of hospital-diagnosed HIV/ AIDS incidence by year, and age group for Chiang Mai; a northern province of Thailand. This study was based on individual hospital case records of HIV/AIDS patients routinely reported from 1997 to 2006. Linear regression to forecast HIV/AIDS incidence rates was used by and age groups in order to prevent disease epidemics that are likely to occur in these areas in the near future. The model contained additive effects associated with the season of the year, district, age group, and the HIV/AIDS incidence rates of previous years, and it can be used to provide useful short-term forecasts. This study found that all coeffi cients of age group factor were signifi cant, as were all those of district factor, except for Samoeng, Hot, Doi Tao, Chai Prakan and . All coeffi cients of year time factor were signifi cant, except for the 1998 factor. The linear regression model obtained an r-square of 74.47%.This study also found that the average incidence rate of AIDS appears to be highest at 2.10 per 10,000 people per year in . Having a model that provides such forecasts of disease outbreaks can provide a useful basis for allocation of resources for disease prevention. Chiang Mai Medical Journal 2010;49(4):161-166. Keywords: AIDS incidence, , linear regression model

Acquired immune defi ciency syndrome Asia new HIV infections declined by nearly 25% (AIDS) is a set of symptoms and infections result- and in South and South East Asia by 10% in the ing from damage to the human immune system same time period. In Eastern Europe, after a caused by the human immunodefi ciency virus dramatic increase in new infections among inject- (HIV).(1) HIV is transmitted through direct con- ing drug users, the epidemic has leveled off con- tact of a mucous membrane or the bloodstream siderably. However, in some countries there are with a bodily fl uid containing HIV, such as blood, signs that new HIV infections are rising again.(3) semen, vaginal fl uid, preseminal fl uid, and breast Statistics regarding AIDS cases occurring in milk.(2) According to new data in the 2009 AIDS Thailand between 1985 and 2006 show that an epidemic update, new HIV infections have been estimated 1,102,628 people (adults and children) reduced by 17% over the past eight years. Since were infected with HIV, and that 558,578 died 2001, the number of new infections in sub-Saha- of AIDS-related complications (data from the ran Africa is approximately 15% lower, which is Thailand A2 Team cited in AIDS Thailand, about 400,000 fewer infections in 2008. In East 2007).(4) Thailand experienced its fi rst case of

Address requests for reprints: Wattanavadee Sriwattanapongse, M.S., Department of Statistics, Faculty of Science, Chiang Mai University, Chiang Mai 50200, Thailand: E-mail:[email protected] Received 26 July 2010, and in revised form 6 August 2010. 162 Sriwattanapongse W, Prasitwattanaseree S, and Khanabsakdi S.

AIDS in 1984. By 1995, approximately 800,000 Doi Saket, Mae Taeng, Mae Rim, Samoeng, Fang, Thais were infected with HIV, and 1 million Mae Ai, Phrao, San Pa Tong, San Kamphaeng, Thais became infected by the year 2000. There San Sai, Hang Dong, Hot, Doi Tao, Omkoi, have been fi ve major epidemic waves: among Saraphi, Wiang Haeng, Chai Prakan, Mae Wang, male homosexuals (beginning 1984-85); intrave- Mae On and Doi Lo. Chiang Mai province has nous drug users (1988); female commercial sex been one of Thailand’s major tourist destinations. workers (1989); male clients of commercial sex It is considered one of the most scenic provinces workers (1990); and housewives and the newborn in the country, due to its mountain ranges, valleys, (1991).(5) fl ora and fauna. Unlike most of Thailand, the The fi rst case of an AIDS patient in the north- climate in the north, especially Chiang Mai, is ern part of Thailand was reported in 1987,(6) cool, fresh and misty. marking the outbreak of the epidemic. In a recent cross-sectional, qualitative study, Angelina MATERIALS AND METHODS Arbisi and Panpanich(7) examined the current atti- tudes toward and motives for acupuncture use and Study Design disuse among people with HIV/AIDS (PHA) in Data used in this study were taken from the northern Thailand. AIDS and STD Prevention and Control Program, Sun et al.(8) discussed the regression analysis in which case information was collected routinely of such studies, and a simple estimating equation in the province by the Chiang Mai Provincial Pub- approach was proposed under the proportional lic Health Offi ce. Data for the years after 1997- hazards model. This method can be easily imple- 2006 are available on computer fi les. There is a mented and does not involve any iteration among record for each case, with that fi elds consist of unknown parameters, as do full likelihood ap- characteristics of each subject including age, gen- proaches proposed in the literature. The asymptot- der, and location (district). ic properties of the proposed regression coeffi cient After cleaning to correct or impute data-entry estimates are derived and an AIDS cohort study is errors, the records for the province were stored analyzed to illustrate the proposed approach. in a Structured Query Language (SQL) database. SQL programs were used to create AIDS case OBJECTIVES counts by year (1997-2006), age group (0-9, 10- 19, 20-29, 30-39, 40-49, 50-59 and 60+ years), This study aimed to fi nd a suitable statistical and district. Incidence rates were computed as the model for predicting the yearly incidence rates of number of cases per 10,000 residents in the dis- reported HIV/AIDS cases by age group, districts trict, according to the 2000 Population and Hous- and year in Chiang Mai province, northern Thai- ing Census of Thailand. land, where a high risk of the disease exists, based on routinely collected data available from Chiang Linear regression model Mai provincial health offi ces. More specifi cally, This method describes data structures with the study aimed to ascertain whether any district multi-categorical determinants (i.e., a multi- had experienced an unusually high incidence of categorical determinant regression model). If

HIV/AIDS, so that preventive disease-control variables are denoted by xx 12 , ,...... x c the extended measures could be put in place. model for Chiang Mai uses the form: Chiang Mai occupies an area of 20,107 square c kilometers and has a total population of 1,500,127 (1) ya=+∑ bxjj people. Neighboring provinces comprise (from j=1 the northeast, clockwise) Chiang Rai, Lampang,

Lamphun, Tak, and Mae Hong Son, most of The constants b 1 ,b 2 ,..... b c are called regression which border Myanmar’s Shan state. Chiang Mai coeffi cients. These measure the province consists of 24 districts: Muang Chiang associations between the determinants and the Mai , Chom Thong, Mae Chaem, Chiang Dao, outcome variable.(9) For this paper, the outcome is AIDS incidence rate in Chiang Mai 163

an incidence rate, as the determinants are district, this model the parameters were constrained in this β age group, and year. study so that α1= 0, 1= 0 and η1= 0. While linear The simplest model is based on linear regres- time trends could be included in the model, they sion, with the outcome variable defi ned as the are less useful for short-term forecasting purposes incidence rate in a cell indexed by district, age in the presence of high serial correlations, and not group, and calendar month (allowing for a sea- considered in this study. sonal effect) as categorical determinants. Such incidence rates generally have positively skewed RESULTS distributions, so it is conventional to transform them by taking logarithms. Since monthly dis- There was a total of 17,535 HIV/AIDS cases ease counts based on small regions are often zero, reported in Chiang Mai province during the ten it is necessary to make some adjustment to avoid year period of 1997-2006. The number of cases taking logarithms of zero. for a particular age group and district varied from The method used to defi ne the outcome is: 0 to 143; and the corresponding maximum disease rates were 11.04 cases per 1,000. The time series of average yearly age-specifi c HIV/AIDS rates ⎛⎞nijt YK ijt =+ ln ⎜⎟ 1 (2) per 1,000 population at risk in all districts of Chi- P ⎝⎠ ang Mai are plotted in Figure 1. where n is the number of disease cases in the cell, In Figure 1, the incidence rates show a marked P is the population at risk, and K is a specifi ed seasonal periodicity, and decrease from very high constant. Then levels in age groups 20-29, 30-39 and 40-49 over the 10-year period; whereas in age groups 0-9, 10- Y P  DE  K  H ijt i j t ijt (3) 19, 50-59 and 60+ years, the reported HIV/AIDS incidence rate shows a slightly less decreas over the same period.

where Nijt is a random variable denoting the Figure 2 shows how the average HIV/AIDS reported number of disease cases in age group incidence rate varied by district in Chiang Mai i, district j and year t for the region of interest, province. The lowest rates occurred in Om Koi and nijt is the corresponding number observed; district (0.15 cases/1,000/year), and the highest

Yijt is the outcome variable specifi ed in Equation were in Hang Dong district (2.10 cases/1,000/

(2) and yijt the corresponding number observed; year). Table 1 shows the results of fi tting a lin- and εijt comprises a set of independent normally ear model to the total number of HIV/AIDS for distributed random variables with a mean of 0. In Chiang Mai province. Equation (3) was used to

AIDScaseincidence:ChiangMai Case/1000 1997Ͳ2006 5 4.5 4 0-9 3.5 10-19 3 20-29 2.5 30-39 2 40-49 1.5 50-59 1 60+ 0.5 Figure 1. Time series of yearly 0 disease rates for each age group 1997 1998 1999 2000 2001 2002 2003 2004 2005 2006 in Chiang Mai province from 1997 to 2006. 164 Sriwattanapongse W, Prasitwattanaseree S, and Khanabsakdi S.

Figure 2. Average annual incidence rates for all AIDS cases in Chiang Mai province from 1997 to 2006.

calculate the log-transformed incidence rates. In and Mae Wang district. Finally, all coeffi cients of each case, K = 10,000, was taken giving reason- year time factor were signifi cant, except for the ably linear normal score plots. The number of year 1998 factor. zero counts is 350 (20.8%). While all compo- nents in the model for each province are statisti- CONCLUSIONS cally signifi cant, and the coeffi cients incorporat- ing further correlations between age groups and This study has shown that AIDS is a problem districts are also quite substantial, the largest in Chiang Mai. The multiple linear regression residual obtained is 2.264, which corresponds to models provide a fi t to age group, districts and 4 cases reported in the age group of more than 60 year. The results from demographic variables years old in during 1997 (incidence show that 61.2% of AIDS patients are male and rate 2.3 per 1,000). However, this residual does 38.8% female; also, that 43.6% of AIDS patients not appear as an outlier on the normal score plot, are aged 30-39 years. This study found that the and the numbers of cases reported in the same dis- average prevalence of AIDS appears to be highest trict and age group were small (0, 1, 1, 0, 0, 0, 1, in 2.1 per 1,000 people per year in Hang Dong 0 and 0, respectively) during the following nine district and the lowest in 0.3 per 1,000 people per years. year in . This study found the Table 1 gives estimated coeffi cients and their highest age-specifi c AIDS incidence rates to be in standard errors for this linear model. When the the group of 30-39 years (28.2 per 1,000 per year). coeffi cient of model was tested by statistics-Z, it The lowest age-specifi c AIDS incidence rate is 1.0 was found that linear regression models and all per 1,000 per year in the group of 10-19 years. coeffi cients of age group factor were signifi cant. During the study period from 1997 to 2006, All coeffi cients of district factor were signifi cant, 17,535 cases of HIV/AIDS were reported in except for Samoeng, Hot, Doi Tao, Chai Prakan Chiang Mai province, with 10,732 (61.2%) and AIDS incidence rate in Chiang Mai 165

Table 1. Results of fi tting a linear model for HIV/AIDs 6,803 (38.8%) being male and female and male, in Chiang Mai province respectively. Most of the AIDS patients were Chiang Mai resident in Muang Chiang Mai, followed by San Coeffi cient Std.Error Pa Tong district and the least number of patients was in . The trend in AIDS patients Constant 1.4428 0.0983 has decreased in females, but are approximately Age group: 0-9 1 1 double that number in males. AIDS patients 10-19 -0.6174 0.0589 decreased mostly at adult age (30-39 years), 20-29 1.6381 0.0589 followed by childhood , old age and teenage, 30-39 1.9904 0.0589 40-49 1.2173 0.0589 consecutively, respectively. 50-59 0.3185 0.0589 The method of back-calculation was used 60+ -0.6011 0.0589 for data on the AIDS incidence rates in Thailand District: 1 Mueang Chiang Mai 1 1 from January 1984 until December 2000. The 2 Chom Thong 0.2324 0.109 data was taken from that released in the UN- AIDS/WHO 2004 Epidemiological Unit, which 3 Mae Chaem -0.6334 0.109 predicts the future trend of the AIDS epidemic 4 Chiang Dao 0.7473 0.109 in Thailand. The simulation indicates a possible 5 Doi Saket 0.3442 0.109 error in the recorded AIDS incidence cases dur- 6 Mae Taeng 0.4507 0.109 ing the early stages of the epidemic.(10) Therefore, 7 Mae Rim 0.4115 0.109 recorded AIDS patients should be recognized by 8 Samoeng 0.0101 0.109 the Department of Public health. 9 Fang 0.4945 0.109 Based on these fi ndings, these researchers pro- 10 Mae Ai 0.6045 0.109 pose that all levels of Thai society join in fi nding 11 Phrao 0.4239 0.109 solutions to promote the prevention of AIDS, and develop strategies to encourage protective behav- 12 San Pa Tong 0.6231 0.109 ior from HIV/AIDS infection. As a result of these 13 San Kamphaeng 0.1496 0.109 and other HIV prevention efforts, and an increase 14 San Sai 0.4347 0.109 in societal awareness of and response to the AIDS 15 Hang Dong 0.7906 0.109 epidemic, new infections in Thailand should be 16 Hot 0.0908 0.109 drastically reduced. 17 Doi Tao 0.1678 0.109 18 Omkoi -0.946 0.109 ACKNOWLEDGEMENTS 19 Saraphi 0.1945 0.109 20 Wiang Haeng 0.3259 0.109 We would like to express our gratitude to the Chiang Mai University and Chiang Mai Provin- 21 Chai Prakan 0.0538 0.109 cial Public Health Offi ce for allowing us to use 22 Mae Wang 0.1191 0.109 the disease data from their offi ce. 23 Mae On 0.3141 0.109 24 Doi Lo 0.4300 0.109 REFERENCES Year: 1997 1 1 1998 -0.1362 0.0704 1. Weiss RA. How does HIV cause AIDS? Science, 1999 -0.2175 0.0704 1993;260:1273–9. 2000 -0.5108 0.0704 2. Centers for Disease Control and Prevention HIV 2001 -0.5766 0.0704 and Its Transmission Atlanta, GA: Divisions of 2002 -0.6649 0.0704 HIV/AIDS Prevention; 1999. 2003 -0.4534 0.0704 3. World Health Organization AIDS Epidemic 2004 -0.6288 0.0704 Update, 2009. 2005 -0.8149 0.0704 2006 -0.8892 0.0704 4. Roberts J. Aids epidemic in Thailand: good news, bad news, and a warning, J Soc Sci 2008; 16:79- R-squared statistic 0.7447 166 Sriwattanapongse W, Prasitwattanaseree S, and Khanabsakdi S.

80. 8. Sun J, Liao Q, Pagano M. Regression analysis of 5. Ruxrungtham K, Phanuphak P. Update on HIV/ doubly censored failure time data with applica- AIDS in Thailand. J Med Assoc Thai 2001; tions to AIDS studies. Biometrics 1999;55:909- 84(Suppl 1):1-17. 14. 6. Subsai K, Kanoksri S, Siwaporn C, Ling H. Neu- 9. McNeil DR. Epidemiological Research Methods, rological complications in AIDS patients: the 1- John Wiley & Sons: New York, 1995. year retrospective study in Chiang Mai Univer- 10. Choon OH, Chye PY, Wah EC. Modelling the sity, Thailand. Eur J Neurol 2004;11:7559. AIDS Epidemic in Thailand: June 13-15, 2nd 7. Arbisi A, Panpanich R. Acupuncture use among IMT-GT Regional Conference on Mathematics, people living with HIV/AIDS in northern Thai- Statistics and Applications, University Sains Ma- land: motives, barriers, and attitudes. J Med As- laysia, Penang: Malaysia, 2006. soc Thai 2008;91:533-4.

อัตราการเกิดอุบัติการณโรคเอดสในจังหวัดเชียงใหม ประเทศไทย

วัฒนาวดี ศรีวัฒนพงศ, ปร.ด. สุคนธ ประสิทธิ์วัฒนเสรี, D.E.A., และ สุรินทร ขนาบศักดิ์, พบ.ม. ภาควิชาสถิติ คณะวิทยาศาสตร มหาวิทยาลัยเชียงใหม

บทคดยั อ วตถั ประสงคุ  เพอศื่ กษารึ ปแบบการแพรู กระจายของอ บุ ตั การณิ การเก ดโรคเอดสิ  ตามป  อำเภอ และกลมอายุ ของผุ ปู วยท อยี่ ในเชู ยงใหมี ซ งเปึ่ นจ งหวั ดหนั งทางภาคเหนึ่ อของประเทศไทยื การ ศกษานึ เปี้ นการศ กษาผึ ปู วยท ไดี่ จากการบ นทั กึ ในระหวางป  พ.ศ. 2540 ถงึ พ.ศ. 2549 โดยการใช ตวแบบอนั กรมเวลาเชุ งเสิ น เพอพยากรณื่ อ บุ ตั การณิ การระบาดของโรคเอดส ในอนาคตตาม อำเภอ และกลมอายุ ุ และเพอปื่ องก นการเกั ดโรคิ ตวแบบเปั นป จจ ยบวกทั เกี่ ยวขี่ องก บฤดั กาลของปู  อำเภอ กลมอายุ ุ และอตราการเกั ดอิ บุ ตั การณิ ของโรคในป ก อนหน าน ี้ และสามารถใชก บการพยากรณั  ระยะสนั้ จากการศกษาพบวึ า สมประสั ทธิ ของปิ์ จจ ยอายั มุ นี ยสำคั ญั สมประสั ทธิ ปิ์ จจ ยดั านอำเภอม ี นยสำคั ญยกเวั นเส ยแตี อำเภอสะเม งิ ฮอด ดอยเตา ไชยปราการ และแมวาง และ สมประสั ทธิ ของิ์ เวลาปท งหมดมั้ นี ยสำคั ญยกเวั นป จจ ยปั  พ.ศ. 2541 คา R-square ของตวแบบการถดถอยเชั งเสิ น เทาก บรั อยละ 74.47 การศกษาในครึ งนั้ ยี้ งพบวั า อตราอั บุ ตั การณิ เฉล ยของโรคเอดสี่ ส งสู ดซุ งึ่ เทาก บั 2.10 ตอ 10,000 ตอป  ซงอยึ่ ในอำเภอหางดงู ตวแบบไดั ให การพยากรณ การระบาดของโรค เพอใชื่ ประโยชน สำหร บกำหนดทรั พยากรเพั อการปื่ องก นโรคเอดสั ต อไป เชยงใหมี เวชสาร 2553; 49(4):161-166. คำสำคญั : อบุ ตั การณิ การเก ดโรคเอดสิ  จงหวั ดเชั ยงใหมี  ตวแบบการถดถอยเชั งเสิ น