圓形分布法在胎兒畸形篩查中的應(yīng)用

楊蕾

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

劉小燕

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

孫紅

（長(zhǎng)江大學(xué)計(jì)算機(jī)科學(xué)學(xué)院，湖北 荊州434023）

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

1 資料來(lái)源與方法

1.1 資料來(lái)源

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

1.2 方法

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

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

1.2.3 統(tǒng)計(jì)公式

1）集中趨勢(shì)的測(cè)度指標(biāo)r：

例外情況：

3）圓標(biāo)準(zhǔn)差S：

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

2 結(jié)果

2.1 胎兒畸形基本情況

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

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

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

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

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

2.3 計(jì)算r值

2.4 查圓分布r值表

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

2.5 求平均角α

169.3693°×365／360＝171.7216°，對(duì)應(yīng)時(shí)間即高峰時(shí)點(diǎn)為6月23日。

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

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

3 討論

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

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

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

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

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

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

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

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

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

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

［6］楊樹(shù)勤.中國(guó)醫(yī)學(xué)百科全書(shū).醫(yī)學(xué)統(tǒng)計(jì)學(xué) ［M］.上海：上海科學(xué)技術(shù)出版社，1985：217－222.

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

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

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

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