圆形分布法在胎儿畸形筛查中的应用

杨蕾

（公安县中医医院妇产科，湖北 公安434300）

刘小燕

（公安县妇幼保健院超声科，湖北 公安434300）

孙红

（长江大学计算机科学学院，湖北 荆州434023）

据文献报道，我国缺陷儿出生率约为4%～6%，每年有80万～100万缺陷儿出生，畸形儿已成为新生儿死亡的主要原因，也给家庭和社会带来巨大的精神痛苦和沉重的经济负担［1］。如何保障我国的优生优育政策，最大限度地减少缺陷胎儿的出生率，提高整体国民素质，是医疗机构妇产科工作者亟待解决的课题。分析研究胎儿畸形的发生规律及变动模式，对探讨病因，提出相应干预措施，进行超前预防和控制，降低畸形胎儿的出生率，保护育龄妇女身心健康，无疑具有重要意义。笔者应用圆形分布统计方法［2］，分析公安县中医医院2007年至2013年三维彩超筛查出的胎儿畸形资料，探讨胎儿畸形发生的季节性规律，旨在为进一步研究胎儿畸形的影响因素提供有用的信息。

1 资料来源与方法

1.1 资料来源

收集整理公安县中医医院2007年至2013年三维彩超筛查出的胎儿畸形病案资料。

1.2 方法

1.2.1 资料整理与基础数据计算 首先将胎儿畸形病案资料进行分类、归纳整理，再将胎儿畸形发生数按月合并，用圆形分布法将发生时间换算成角度。胎儿畸形发生月份为LMP月份。为了计算方便，以元旦零时为零角度，通过三角函数代换原理，以各月发病数乘以各月的度数来计算基础数据。

1.2.2 月中位角计算 先求出各月的月中点，以元旦零时为零角度，则1月份有31d，月中点距零为15.5d，2月份有28d，月中点距零为31＋14＝45d，依次类推。再将时间天数换算成角度，全年365d，相当于圆周360°，则1d相当于360／365＝0.9863°，1月份月中点为15.5d，则月中位角为：

1.2.3 统计公式

1）集中趋势的测度指标r：

例外情况：

3）圆标准差S：

1.2.4 实际计算 本文在具体计算时，应用Excel软件的函数功能进行数据处理。首先将上述公式及相应的一些关系式按Excel软件的函数格式［3］输入Excel表，再将某病各月的病例数输入，系统即可自动给出所求的X值、Y值、值、r值、S值等结果数据。

2 结果

2.1 胎儿畸形基本情况

本研究收集整理公安县中医医院2007年至2013年筛查的胎儿畸形病案资料见表1。公安县中医医院2007年至2013年共筛查孕妇12177例，检出并经引产证实为畸形胎儿259例，畸形率2.13%，以内脏畸形、脑脊畸形和肛脐畸形为主，其中内脏畸形占24.71%，脑脊畸形占23.16%，肛脐畸形占21.62%。胎儿畸形总的构成比之间，经统计处理，差异有统计学意义（P＜0.01）。各年份筛查的胎儿畸形率经统计处理，差异无统计学意义（P＞0.05），即各年份的畸形率无统计学差异。

表1 胎儿畸形基本情况统计

2.2 各月发生频数及圆分布情况

将表1资料中的胎儿畸形发生数按月合并，输入Excel表，同时在相应的单元格输入计算公式，将发生时间及其与发生频数换算成角度，计算结果见表2。

表2 各月发生频数及圆分布情况统计

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 求高峰时区α±S

即：169.3693°±S＝169.3693±75.9339＝93.4353°～245.3032°，胎儿畸形发生的高峰时区为每年3月25日至9月7日，表明3月下旬至9月上旬受孕增加了发生胎儿畸形的风险。

3 讨论

众所周知，Gregg于1941年证实先天性白内障的季节性与孕早期母体感染风疹病毒有关，揭示了环境因素的致畸效应。1961年Edwands将圆形统计思想用于先天畸形的季节性分析［4］，使得推算风险度最大的受孕时间，采取干预措施进行控制成为可能。大量的文献报道，近代三维超声技术的发展，在早孕期胎儿畸形筛查中发挥了巨大作用［5］。

文献记载，反映疾病季节性规律的研究的方法，有季节性指数、余弦函数、季节性升高所占的月份数、病例占全年的百分比、线状图等［6］，但这些方法都不能用一个简单的指标反映全年的变化，而圆形分布统计方法则综合了全年各月疾病发生情况，比较准确地反映了这种变化，是描述疾病季节性特征较完善的方法之一［7］。

笔者应用圆形统计方法，分析公安县中医医院2007年至2013年筛查的胎儿畸形资料，探讨胎儿畸形发生的季节性规律，结果令人满意。按平均角估计，胎儿畸形发生的高峰时点为6月23日，高峰区间为3月下旬至9月上旬，提示这个区间受孕增加了发生胎儿畸形的风险。

本研究结果表明，离散程度指标r＝0.4155，P＜0.01，表明胎儿畸形的发生存在集中趋势［8］。提示育龄妇女在3月25日至9月7日期间受孕，要倍加注意防控发生胎儿畸形的风险和规避因素。

本组7年资料259例胎儿畸形的季节变动周期是否符合余弦曲线模式，发生率与受孕月份是否呈指数曲线关系［9］，尚未深入研究。一般来说，确定型模型在探讨疾病的发病规律、进行预测方面，比通常的统计指标能提供更多的信息［10］。如何建立胎儿畸形发生的季节性规律模式，是我们今后进一步研究的方向。

［1］李胜利.胎儿畸形产前超声诊断学 ［M］.北京：人民军医出版社，2004：86.

［2］郭祖超.医用数理统计方法 ［M］.3版.北京：人民卫生出版社，1988：122－135.

［3］耿贯一.流行病学上册 ［M］.北京：人民卫生出版社，1985：28－31.

［4］王成科.利用数学模型探讨计划免疫前后麻疹发病季节特征 ［J］.现代预防医学，1998，25（3）：301－302.

［5］周艳.三维彩超技术在早中孕胎儿畸形筛查中的临床价值 ［J］.河南医学研究，2013，22（2）：239－240.

［6］杨树勤.中国医学百科全书.医学统计学 ［M］.上海：上海科学技术出版社，1985：217－222.

［7］郑彦玲，郑艳梅，张利萍.基于随机时间序列法的新疆乙肝发病率预测 ［J］.数理医药学杂志，2013，6：631－635.

［8］张利娟，徐志敏，郑浩，等.湖区5省急性血吸虫感染时空聚集性分析 ［J］.中国血吸虫病防治杂志，2014，26（3）：245－249.

［9］任飞林，仝振东，唐安，等.应用圆形分布法分析舟山市手足口病发病高峰变化 ［J］.浙江预防医学，2014，26（7）：683－686.

［10］薛黎坚，王文明，王华，等.应用圆形分布法分析戊肝发病高峰 ［J］.数理医药学杂志，2014，27（3）：262－263.