葉孟良，鐘曉妮

（重慶醫(yī)科大學(xué)公共衛(wèi)生學(xué)院 400016）

在醫(yī)學(xué)領(lǐng)域中，有些數(shù)據(jù)以一年中的月、日或一晝夜中的時、分來表示，可以研究是否有集中于某一時刻的傾向，這一類時間性的資料不屬于線性的正態(tài)分布，可轉(zhuǎn)化成角度資料來處理，采用圓形分布法進(jìn)行統(tǒng)計(jì)分析［1－2］。本文用圓形分布法探討重慶市2009年各系統(tǒng)疾病的住院時間分布特征，為合理配置臨床醫(yī)療資源及疾病的預(yù)防控制提供依據(jù)。

1 資料與方法

1.1 一般資料 重慶市2009年住院患者資料，共882 727例。按ICD－10疾病編碼提取出各系統(tǒng)疾病總例數(shù)［3］，見表1。

1.2 方法 應(yīng)用圓形分布法進(jìn)行統(tǒng)計(jì)分析［4－5］，根據(jù)圓形分布統(tǒng)計(jì)模型，計(jì)算各月天數(shù)的組中值，再轉(zhuǎn)換成角度（α），計(jì)算其正弦值（sinα）、余弦值（cosα），并分別與每月的病例數(shù)（f）相乘，然后應(yīng)用下列公式計(jì)算平均角度和角度標(biāo)準(zhǔn)差（s），用±s推算平均住院高峰期：x＝（Σfcosα）／n，y＝（Σfsinα）／n，r＝最后將所求出的角度均值）及標(biāo)準(zhǔn)差（s）化為相應(yīng)時間，并對r作顯著性檢驗(yàn)，當(dāng)n≤500，查r分布界值表，但當(dāng)n足夠大時，r的分布符合瑞利分布［6－7］。由于r的值在［0，1］區(qū)間內(nèi)，即r又符合Alpha＝2，Beta＝1的韋伯分布，可以利用Excel精確計(jì)算圓形分布角均數(shù)的臨界概率。

1.3 統(tǒng)計(jì)學(xué)處理 數(shù)據(jù)的提取采用SAS9.1，圓形分布的計(jì)算采用Excel2003完成。

表1 2009年重慶市住院患者疾病構(gòu)成（%）及順位

2 結(jié) 果

以傳染病和寄生蟲病為例，詳細(xì)闡述其應(yīng)用圓形分布分析住院患者的季節(jié)性分布。

2.1 重慶市2009年傳染病和寄生蟲病的月份分布 全年共有28 597例住院患者，其中4月份有2 979例住院患者，占10.42%，所占的比重最高，1月份最低，占6.01%，可見住院患者有集中于某個時間點(diǎn)的傾向，則這一傾向可用平均角表示其集中位置［8］，見表2。

表2 重慶市2009年傳染病和寄生蟲病住院患者月份分布

2.2 圓形分布分析結(jié)果 見表3。

表3 重慶市2009年傳染性疾病與寄生蟲病住院時間分布

圓形分布的平均角及標(biāo)準(zhǔn)差的計(jì)算：x＝（Σfcosα）／n＝－0.107，y＝（Σfsinα）／n＝－0.029，r＝＝0.110。由于x＜0，所以平均角＝180°＋acstg（y／x）＝164.95°，轉(zhuǎn)化為時間是第167.24天，即6月16日，此即為該地區(qū)的傳染病與寄生蟲病患者住院的理論高峰值。s＝57.295 8×12hr＝120.30°。轉(zhuǎn)化為天數(shù)是121.97d。

查“圓形分布r界值表”，因?yàn)閞＞r0.05（500）＞r0.05（28 597），所以P＜0.05，平均角有統(tǒng)計(jì)學(xué)意義，說明傳染性疾病和寄生蟲病住院分布存在季節(jié)性集中趨勢。

同理，應(yīng)用圓形分布法分析2009年重慶市住院患者中構(gòu)成比前十位中其余九類疾病的季節(jié)性，得到呼吸系統(tǒng)疾病、消化系統(tǒng)疾病、損傷與中毒存、妊娠、分娩和產(chǎn)褥期、泌尿生殖系統(tǒng)疾病、腫瘤、眼和附器疾病在住院月份上存在集中趨勢（P＜0.05）；影響健康狀態(tài)和與保健機(jī)構(gòu)接觸的因素以及循環(huán)系統(tǒng)不存在集中趨勢（P＞0.05），見表4。

表4 2009年重慶市住院患者圓形分布均角的假設(shè)檢驗(yàn)結(jié)果

3 討 論

疾?。ɑ蛏憩F(xiàn)象）的發(fā)生時間在晝夜24h的分布稱為晝時性資料，在全年12個月的分布為季節(jié)性資料。這類資料特點(diǎn)是：（1）角度、鐘點(diǎn)等周而復(fù)始，沒有準(zhǔn)確零點(diǎn)；（2）數(shù)值大小的說法是任意的［9］。這種資料是不能直接進(jìn)行算術(shù)運(yùn)算的，否則會產(chǎn)生謬誤。傳統(tǒng)的方法是用各月發(fā)病數(shù)或率繪制流行曲線。但流行曲線只能對發(fā)病季節(jié)性進(jìn)行粗略的定性描述，不僅得不到確切的發(fā)病高峰時點(diǎn)，更不便進(jìn)行季節(jié)性強(qiáng)弱、高峰期早晚的比較［10］。住院患者的資料，屬于季節(jié)性資料，圓形分布法能真實(shí)、簡明地反映全年各月各個系統(tǒng)疾病的住院情況，結(jié)果判定不受樣本大小的影響［11］。因此可以應(yīng)用圓形分布分析各個系統(tǒng)住院患者是否存在高峰期以及存在的時點(diǎn)，為臨床合理配置資源與疾病的預(yù)防控制提供依據(jù)。

患者住院具有隨機(jī)性，各系統(tǒng)疾病全年都有，常有不確定性，但進(jìn)行大樣本分析則可揭示其規(guī)律。分析各系統(tǒng)疾病的季節(jié)分布，對揭示疾病的發(fā)生時間分布規(guī)律及合理分配各臨床科室人力、物力、財(cái)力資源以及疾病的預(yù)防控制提供依據(jù)，避免資源分配不合理而浪費(fèi)資源。國內(nèi)尚無利用圓形分布綜合分析各系統(tǒng)疾病的季節(jié)分布的文獻(xiàn)報(bào)道。

本研究應(yīng)用圓形分布法探討了2009年重慶市各系統(tǒng)住院患者的季節(jié)分布特征，傳染病和寄生蟲病、呼吸系統(tǒng)疾病、消化系統(tǒng)疾病、損傷與中毒、妊娠、分娩和產(chǎn)褥期、泌尿生殖系統(tǒng)疾病、腫瘤、眼和附器疾病在住院月份上存在集中趨勢。這些系統(tǒng)或類別疾病存在集中趨勢，但其高峰值各不相同，因此醫(yī)院可以根據(jù)其高峰期合理安排相應(yīng)科室的物力、財(cái)力及人力，同時為疾病的預(yù)防控制一些有益的線索。

應(yīng)用圓形分布分析住院患者的季節(jié)性符合統(tǒng)計(jì)學(xué)分析原理方法，分析結(jié)果能為臨床資源的配置提供依據(jù)，但是圓形分布法本身只是描述了各系統(tǒng)疾病的時間分布特征，不能揭示其影響因素，不能解決資料的偏性問題［12］。因而在此基礎(chǔ)上，需進(jìn)一步進(jìn)行流行病學(xué)調(diào)查，分析其影響因素，為疾病的防治提供依據(jù)。

［1］ 陸守曾.醫(yī)學(xué)統(tǒng)計(jì)學(xué)［M］.北京：中國統(tǒng)計(jì)出版社，2002：264－268.

［2］ 何清波，蘇炳華，錢元.醫(yī)學(xué)統(tǒng)計(jì)學(xué)［M］.上海：上海科學(xué)技術(shù)文獻(xiàn)出版社，2002：231－242.

［3］ 衛(wèi)生部統(tǒng)計(jì)學(xué)信息中心，北京協(xié)和醫(yī)院世界衛(wèi)生組織疾病分類合作中心.國際疾病分類（ICD－10）應(yīng)用指導(dǎo)手冊［M］.北京：中國協(xié)和醫(yī)科大學(xué)出版社，2001.

［4］ 王慶昌，李欣.圓形分布分析的EXCEL實(shí)現(xiàn)［J］.中國衛(wèi)生統(tǒng)計(jì)，2006，23（5）：448－449.

［5］ 朱凱，李悅，劉偉新.圓形分布資料平均角置信區(qū)間的Bootstrap方法實(shí)現(xiàn)［J］.中國衛(wèi)生統(tǒng)計(jì)，2011，28（4）：412－413.

［6］ 朱歲松.圓形分布角均數(shù)的顯著性檢驗(yàn)［J］.數(shù)理醫(yī)藥學(xué)雜志，2009，22（1）：53－55.

［7］ 田考聰.描述性統(tǒng)計(jì)分冊［M］.北京：人民衛(wèi)生出版社，2004：108－110.

［8］ 楊樹勤.中國醫(yī)學(xué)百科全書醫(yī)學(xué)統(tǒng)計(jì)學(xué)［M］.上海：上?？茖W(xué)技術(shù)出版社，1982：217－222.

［9］ 石小河，王玉環(huán).用圓形分布法分析職業(yè)傷害的時間分布特征［J］.疾病控制雜志，2006，3（1）：30－32.

［10］ 陳勇，殷菲.應(yīng)用圓形分布法探討江蘇省流腦的發(fā)病季節(jié)性及其變化趨勢［J］.現(xiàn)代預(yù)防醫(yī)學(xué)，2008，35（14）：2630－2631.

［11］ 甘仰本，廖征，蔡軍.應(yīng)用集中度和圓形分布法分析南昌市1985～2008年流行性腦脊髓膜炎發(fā)病季節(jié)性［J］.中國衛(wèi)生統(tǒng)計(jì)，2010，27（4）：379－380.

［12］ Mittleman MA，Maldonado G，Gerberich SG，et al.Alternative approaches to analytical designs in occupational injury epidemiology［J］.Am J Ind Med，1997，32（2）：129－141.