圓形分布法在胎兒畸形篩查中的應用

楊蕾

（公安縣中醫(yī)醫(yī)院婦產(chǎn)科，湖北 公安434300）

劉小燕

（公安縣婦幼保健院超聲科，湖北 公安434300）

孫紅

（長江大學計算機科學學院，湖北 荊州434023）

據(jù)文獻報道，我國缺陷兒出生率約為4%～6%，每年有80萬～100萬缺陷兒出生，畸形兒已成為新生兒死亡的主要原因，也給家庭和社會帶來巨大的精神痛苦和沉重的經(jīng)濟負擔［1］。如何保障我國的優(yōu)生優(yōu)育政策，最大限度地減少缺陷胎兒的出生率，提高整體國民素質(zhì)，是醫(yī)療機構(gòu)婦產(chǎn)科工作者亟待解決的課題。分析研究胎兒畸形的發(fā)生規(guī)律及變動模式，對探討病因，提出相應干預措施，進行超前預防和控制，降低畸形胎兒的出生率，保護育齡婦女身心健康，無疑具有重要意義。筆者應用圓形分布統(tǒng)計方法［2］，分析公安縣中醫(yī)醫(yī)院2007年至2013年三維彩超篩查出的胎兒畸形資料，探討胎兒畸形發(fā)生的季節(jié)性規(guī)律，旨在為進一步研究胎兒畸形的影響因素提供有用的信息。

1 資料來源與方法

1.1 資料來源

收集整理公安縣中醫(yī)醫(yī)院2007年至2013年三維彩超篩查出的胎兒畸形病案資料。

1.2 方法

1.2.1 資料整理與基礎數(shù)據(jù)計算 首先將胎兒畸形病案資料進行分類、歸納整理，再將胎兒畸形發(fā)生數(shù)按月合并，用圓形分布法將發(fā)生時間換算成角度。胎兒畸形發(fā)生月份為LMP月份。為了計算方便，以元旦零時為零角度，通過三角函數(shù)代換原理，以各月發(fā)病數(shù)乘以各月的度數(shù)來計算基礎數(shù)據(jù)。

1.2.2 月中位角計算 先求出各月的月中點，以元旦零時為零角度，則1月份有31d，月中點距零為15.5d，2月份有28d，月中點距零為31＋14＝45d，依次類推。再將時間天數(shù)換算成角度，全年365d，相當于圓周360°，則1d相當于360／365＝0.9863°，1月份月中點為15.5d，則月中位角為：

1.2.3 統(tǒng)計公式

1）集中趨勢的測度指標r：

例外情況：

3）圓標準差S：

1.2.4 實際計算 本文在具體計算時，應用Excel軟件的函數(shù)功能進行數(shù)據(jù)處理。首先將上述公式及相應的一些關(guān)系式按Excel軟件的函數(shù)格式［3］輸入Excel表，再將某病各月的病例數(shù)輸入，系統(tǒng)即可自動給出所求的X值、Y值、值、r值、S值等結(jié)果數(shù)據(jù)。

2 結(jié)果

2.1 胎兒畸形基本情況

本研究收集整理公安縣中醫(yī)醫(yī)院2007年至2013年篩查的胎兒畸形病案資料見表1。公安縣中醫(yī)醫(yī)院2007年至2013年共篩查孕婦12177例，檢出并經(jīng)引產(chǎn)證實為畸形胎兒259例，畸形率2.13%，以內(nèi)臟畸形、腦脊畸形和肛臍畸形為主，其中內(nèi)臟畸形占24.71%，腦脊畸形占23.16%，肛臍畸形占21.62%。胎兒畸形總的構(gòu)成比之間，經(jīng)統(tǒng)計處理，差異有統(tǒng)計學意義（P＜0.01）。各年份篩查的胎兒畸形率經(jīng)統(tǒng)計處理，差異無統(tǒng)計學意義（P＞0.05），即各年份的畸形率無統(tǒng)計學差異。

表1 胎兒畸形基本情況統(tǒng)計

2.2 各月發(fā)生頻數(shù)及圓分布情況

將表1資料中的胎兒畸形發(fā)生數(shù)按月合并，輸入Excel表，同時在相應的單元格輸入計算公式，將發(fā)生時間及其與發(fā)生頻數(shù)換算成角度，計算結(jié)果見表2。

表2 各月發(fā)生頻數(shù)及圓分布情況統(tǒng)計

2.3 計算r值

2.4 查圓分布r值表

本例n＝259，r＝0.4155 P＜0.01。

2.5 求平均角α

169.3693°×365／360＝171.7216°，對應時間即高峰時點為6月23日。

2.6 求高峰時區(qū)α±S

即：169.3693°±S＝169.3693±75.9339＝93.4353°～245.3032°，胎兒畸形發(fā)生的高峰時區(qū)為每年3月25日至9月7日，表明3月下旬至9月上旬受孕增加了發(fā)生胎兒畸形的風險。

3 討論

眾所周知，Gregg于1941年證實先天性白內(nèi)障的季節(jié)性與孕早期母體感染風疹病毒有關(guān)，揭示了環(huán)境因素的致畸效應。1961年Edwands將圓形統(tǒng)計思想用于先天畸形的季節(jié)性分析［4］，使得推算風險度最大的受孕時間，采取干預措施進行控制成為可能。大量的文獻報道，近代三維超聲技術(shù)的發(fā)展，在早孕期胎兒畸形篩查中發(fā)揮了巨大作用［5］。

文獻記載，反映疾病季節(jié)性規(guī)律的研究的方法，有季節(jié)性指數(shù)、余弦函數(shù)、季節(jié)性升高所占的月份數(shù)、病例占全年的百分比、線狀圖等［6］，但這些方法都不能用一個簡單的指標反映全年的變化，而圓形分布統(tǒng)計方法則綜合了全年各月疾病發(fā)生情況，比較準確地反映了這種變化，是描述疾病季節(jié)性特征較完善的方法之一［7］。

筆者應用圓形統(tǒng)計方法，分析公安縣中醫(yī)醫(yī)院2007年至2013年篩查的胎兒畸形資料，探討胎兒畸形發(fā)生的季節(jié)性規(guī)律，結(jié)果令人滿意。按平均角估計，胎兒畸形發(fā)生的高峰時點為6月23日，高峰區(qū)間為3月下旬至9月上旬，提示這個區(qū)間受孕增加了發(fā)生胎兒畸形的風險。

本研究結(jié)果表明，離散程度指標r＝0.4155，P＜0.01，表明胎兒畸形的發(fā)生存在集中趨勢［8］。提示育齡婦女在3月25日至9月7日期間受孕，要倍加注意防控發(fā)生胎兒畸形的風險和規(guī)避因素。

本組7年資料259例胎兒畸形的季節(jié)變動周期是否符合余弦曲線模式，發(fā)生率與受孕月份是否呈指數(shù)曲線關(guān)系［9］，尚未深入研究。一般來說，確定型模型在探討疾病的發(fā)病規(guī)律、進行預測方面，比通常的統(tǒng)計指標能提供更多的信息［10］。如何建立胎兒畸形發(fā)生的季節(jié)性規(guī)律模式，是我們今后進一步研究的方向。

［1］李勝利.胎兒畸形產(chǎn)前超聲診斷學 ［M］.北京：人民軍醫(yī)出版社，2004：86.

［2］郭祖超.醫(yī)用數(shù)理統(tǒng)計方法 ［M］.3版.北京：人民衛(wèi)生出版社，1988：122－135.

［3］耿貫一.流行病學上冊 ［M］.北京：人民衛(wèi)生出版社，1985：28－31.

［4］王成科.利用數(shù)學模型探討計劃免疫前后麻疹發(fā)病季節(jié)特征 ［J］.現(xiàn)代預防醫(yī)學，1998，25（3）：301－302.

［5］周艷.三維彩超技術(shù)在早中孕胎兒畸形篩查中的臨床價值 ［J］.河南醫(yī)學研究，2013，22（2）：239－240.

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

［7］鄭彥玲，鄭艷梅，張利萍.基于隨機時間序列法的新疆乙肝發(fā)病率預測 ［J］.數(shù)理醫(yī)藥學雜志，2013，6：631－635.

［8］張利娟，徐志敏，鄭浩，等.湖區(qū)5省急性血吸蟲感染時空聚集性分析 ［J］.中國血吸蟲病防治雜志，2014，26（3）：245－249.

［9］任飛林，仝振東，唐安，等.應用圓形分布法分析舟山市手足口病發(fā)病高峰變化 ［J］.浙江預防醫(yī)學，2014，26（7）：683－686.

［10］薛黎堅，王文明，王華，等.應用圓形分布法分析戊肝發(fā)病高峰 ［J］.數(shù)理醫(yī)藥學雜志，2014，27（3）：262－263.