扩展功能
文章信息
- 夏小雪, 杨陆一, 滕蓉, 宁磊, 王守东, 牡琦丽
- XIA Xiaoxue, YANG Luyi, TENG Rong, NING Lei, WANG Shoudong, MU Qilli
- 基于遗传算法的骨性Ⅱ类错牙合患者垂直向颅面部骨关系的预测
- Predicti on onvertical craniofacial bone relationship in patients with skeletal class Ⅱ malocclusion based on genetic algorithms method
- 吉林大学学报(医学版), 2017, 43(04): 794-799
- Journal of Jilin University (Medicine Edition), 2017, 43(04): 794-799
- 10.13481/j.1671-587x.20170425
-
文章历史
- 收稿日期: 2016-12-05
骨性Ⅱ类错牙合是临床上常见的错牙合畸形,患者表现出明显的凸面型,其主要原因是颅面部骨的矢状向不调,表现为上颌前突或正常,下颌后缩,或者两者兼有[1]。研究[2]表明:与正常牙合比较,骨性Ⅱ类错牙合患者上颌骨处于回缩位置,而此时的下颌位置更加后缩,发育不良。不同的垂直向骨型会影响骨性Ⅱ类错牙合患者临床表型的严重程度,尤其是对于具有生长发育能力的青少年来说,垂直向高度会随生长发育发生改变。无论是矫治方案的制订,还是治疗过程中垂直向的控制,其垂直向的特征都是临床医师要考虑的因素。因此,有关骨性Ⅱ类错牙合患者垂直向方面的研究受到了国内外学者的广泛关注。
遗传算法是一种源于自然选择原理和自然遗传机制的全局搜索算法[3]。其实质就是在有限的、离散的数学结构上,寻找一个满足给定约束条件并使目标函数值达到最大或最小的解[3]。遗传算法作为具有全局搜索功能的优化算法,其中组合优化问题是其主要的应用领域之一,并得到了广泛的应用和研究[4-6]。本实验采用遗传算法,对骨性Ⅱ类错牙合患者不同垂直类型的颅面骨建立预测方程,目前国内外尚未见相关报道。本文作者通过分析骨性Ⅱ类错牙合畸形患者不同垂直类型各自形态特点及相关性,为其临床的诊断和治疗提供参考, 同时为颅面部生长发育方面的研究提供新思路。
1 资料与方法 1.1 初步筛选选取2012—2015年在吉林大学口腔医院正畸科就诊的安氏Ⅱ类即骨性Ⅱ类错牙合畸形患者300例,均为初诊病例,按以下条件进行初步筛选:汉族,年龄10~18岁;两侧尖牙、磨牙为Ⅱ类关系,上下齿槽座角(ANB角)大于5°;牙列完整(第三磨牙无要求),牙齿形态、数目无异常,无系统性疾病及遗传病史,无正畸治疗史,无颞下颌关节疾病及颅颌面部外伤史;病历完整,头颅侧位片清晰。
1.2 分组根据投影测量结果,将初步筛选后符合要求的200例患者以下颌平面角(FH-MP和SN-MP)值分为3组。① 高角组:FH-MP>32°,SN-MP>40°;② 均角组:22° < FH-MP < 32°,29° < SN-MP < 40°;③ 低角组:FH-MP < 22°,SN-MP < 29°。去除不符合标准的样本,最终获得样本155例。其中高角组50例,男性20例,女性30例,平均年龄(13.7±2.4) 岁;均角组58例,男性32例,女性26例,平均年龄(13.2± 2.5) 岁;低角组47例,男性22例,女性25例,平均年龄(13.4±2.8) 岁。从3组患者中分别随机取出5例为检验样本,余下为实验样本。实验样本用于遗传算法优化公式的获得,检验样本用于误差分析。
1.3 实验方法头颅侧位片经同一名医师严格按照要求由同一台设备拍摄,并存储于Winceph8.0头影测量软件的图像管理系统(日本Rise公司)中进行测量分析(主要分析方法为Coben分析法和Steiner & Tweed北医大分析法),由研究者本人在计算机上进行定点测量,每次测量间隔1周,共测3次,取平均值,记为最终结果。投影测量参数为Winceph8.0头影测量软件的图像管理系统中包含项目。
1.4 遗传算法参数优化方法遗传算法主要参数为:种群规模80 000,选择率0.05,交叉率0.1,变异率0.2,进化代数30。需要优化的方程结构为:y=xx(1)+xx(2)×p(1, i),y=xx(1)+xx(2)×p(1, i)+xx(3)×p(2, i),y= xx(1)+xx(2)×p(1, i)+xx(3)× p(2, i)+xx(4)×p(3, i)。式中,xx(i)为需要优化的参数,p(1, i)、p(2, i)、p(3, i)为影响因子。根据自变量个数选择不同方程进行优化。使用遗传算法对方程进行参数优化(数据来源于实验样本),并得出预测方程。对检验样本进行遗传算法预测值与实测值误差比较,并对结果进行统计学分析。
1.5 统计学分析采用SPSS22.0统计软件进行统计学分析。所有测量参数均以x±s表示,对同一类型不同性别患者颅面部垂直向的相关项参数[前面高(N-Me)、上面高(N-Ans)、前鼻棘至上中切牙切缘垂直距离(Ans-U1)、下中切牙切缘至颏下点的垂直距离(L1-Me)、下面高(Ans-Me)、后面高(N-Go)、鼻根点至蝶鞍点的垂直距离(N-S)、蝶鞍点至关节点的垂直距离(S-Ar)、下颌升支高度(Ar-Go)、蝶鞍点至下颌角点的垂直距离(S-Go)、下颌角(∠Go)、面角(Facial angle)和后面高前面高比(S-Go/N-Me) ]比较采用独立样本t检验,如不同性别患者各测量参数比较差异无统计学意义,则可将不同性别组合并进行分析。上述测量参数3组间比较采用单因素方差分析,组间两两比较采用SNK-q检验。分别对3组上述测量参数进行逐步回归分析,确定影响因子。对检验样本进行遗传算法预测值与实测值误差比较,采用独立样本t检验。以α=0.05为检验水准。
2 结果 2.1 骨性Ⅱ类错牙合畸形患者颅面部垂直向头影测量结果对同一类型不同性别患者颅面部垂直向的相关参数进行独立样本t检验,结果显示:不同性别患者间各测量参数比较差异无统计学意义(P > 0.05),因此将不同性别组参数合并进行分析。3组间S-Go/N-Me比较差异有统计学意义(P < 0.01),表明以FH-MP和SN-MP值作为垂直向分组条件具有可行性并且准确度高。高角组、均角组和低角组3组间比较,前面高、上面高、下面高、后面高和下颌升支高度差异均有统计学意义(P < 0.01),前面高、上面高、下面高逐渐减小,后面高、下颌升支高度逐渐增大。高角组前面高明显高于均角组和低角组(P < 0.01),前上面高和前下面高明显高于低角组(P < 0.01),前下面高高于均角组(P < 0.05),下颌升支高度较低角组明显缩短(P < 0.01)。低角组前面高、上面高和下面高明显低于高角组(P < 0.01),下颌升支高和后面高明显高于均角组和高角组(P < 0.01)。见表 1。
(x±s) | |||||||||||||||
Group | n | N-Me(l/mm) | N-Ans(l/mm) | Ans-U1(l/mm) | L1-Me(l/mm) | Ans-Me(l/mm) | N-Go(l/mm) | n | N-S(l/mm) | S-Ar(l/mm) | Ar-Go(l/mm) | S-Go(l/mm) | ∠Go(θ/°) | Facial angle(θ/°) | S-Go/N-Me(η/%) |
High-angle | 50 | 120.56±8.60 | 52.86±5.78 | 29.58±4.34 | 36.79±5.52 | 62.40±8.61 | 80.38±7.37 | 50 | 7.96±3.63 | 30.46±4.73 | 38.96±6.30 | 72.59±6.18 | 126.26±6.11 | 82.33±2.81 | 60.03±3.19 |
Average-angle | 48 | 110.16±9.01** | 51.03±4.27 | 27.52±3.07** | 35.74±3.31 | 59.13±5.89* | 78.46±6.58** | 58 | 7.96±3.28 | 30.61±3.45 | 39.89±4.39* | 70.68±6.50** | 121.38±5.48** | 84.76±2.89** | 64.51±3.27** |
Low-angle | 47 | 108.66±7.95** | 50.41±6.65** | 26.31±2.49** | 36.23±3.42 | 57.61±4.68** | 83.06±5.96** | 47 | 8.99±3.85 | 30.89±3.14 | 43.18±5.03**△ | 74.05±6.79** | 113.59±5.46** | 87.71±3.08** | 73.50±3.04** |
P < 0.01 vs average-angle group.. |
采用逐步回归分析法确定与颅面结构相关的影响因子。高角组:年龄与Ans-U1呈正相关关系(r=0.470, P < 0.05),N-Me与Ans-Me呈正相关关系(r=0.964, P < 0.05),Ans-U1与Ans-Me呈正相关关系(r=0.805, P < 0.05)、与面角呈负相关关系(r=-0.023, P < 0.05),N-Go与N-Me呈正相关关系(r=0.926, P < 0.05), L1-Me与Ans-Me呈正相关关系(r=0.898, P < 0.05)、与∠Go呈负相关关系(r=-0.468, P < 0.05);均角组:年龄与N-Me呈正相关关系(r=0.531, P < 0.05),Ans-U1与Ans-Me呈正相关关系(r=0.878, P < 0.05)、与面角成负相关关系,Ans-Me与N-Me呈正相关关系(r=0.920, P < 0.05)、与N-Ans呈负相关关系(r=-0.560, P < 0.05)、与Ar-Go呈负相关关系(r=-0.652, P < 0.05),N-Go与S-Go呈正相关关系(r=0.867, P < 0.05)、与N-Ans呈正相关关系(r=0.252, P < 0.05)、与L1-Me呈正相关关系(r=0.754, P < 0.05), S-Ar与S-Go呈正相关关系(r=0.671, P < 0.05)、与Ar-Go呈负相关关系(r=-0.250, P < 0.05)、与L1-Me呈正相关关系(r=0.552, P < 0.05);低角组:年龄与S-Go呈正相关关系(r=0.602, P < 0.05)、与∠Go呈负相关关系(r=-0.346, P < 0.05)、与L1-Me呈正相关关系(r=0.576, P < 0.05),N-Me与Ans-Me呈正相关关系(r=0.869, P < 0.05)、与N-Go呈正相关关系(r=0.859, P < 0.05)、与面角呈负相关关系(r=-0.177, P < 0.05),N-Ans与N-Me呈正相关关系(r=0.605, P < 0.05)、与Ans-U1呈负相关关系(r=-0.113, P < 0.05),Ans-Me与N-Me呈正相关关系(r=0.869, P < 0.05)、与面角呈正相关关系(r=0.070, P < 0.05)、与Ans-U1呈正相关关系(r=0.785, P < 0.05),N-Go与N-Me呈正相关关系(r=0.859, P < 0.05)、与S-Go呈正相关关系(r=0.829, P < 0.05)。
2.3 相关方程的建立采用遗传算法建立了不同垂直类型骨性Ⅱ类患者颅面部垂直向各标志点的定量关系方程。高角组:年龄=5.8836+0.2690×Ans-U1, N-M=22.0266+1.4945×Ans-Me, Ans-U1= 34.9594+0.4545×Ans-Me-0.4097× Facial angle,N-Go=-4.5882+0.4724×N-Me, L1-Me=-12.5905+0.5322×Ans-Me+0.1243×∠Go。均角组:年龄=-2.9441+0.1468×N-Me, Ans-U1=18.917 0+0.4764×Ans-Me-0.2302×Facial angle, Ans-Me=-0.6205+1.0145×N-Me-0.9741×N-Ans-0.0576×Ar-Go, N-Go=1.6311+0.8978×S-Go+ 0.9197×N-S+0.1688×L1-Me, S-Ar=-1.8231+0.8453×S-Go-0.8670×Ar-Go+0.2024×L1-Me。低角组:年龄=11.7406+0.1527×S-Go-0.1699×∠Go+0.2525×L1-Me, N-Me=61.1530+0.9643×Ans-Me+0.6286×N-Go-0.6892×Facial angle, N-Ans=-4.9492+1.0658×N-Me-2.3165×Ans-U1, Ans-Me=-25.1800+0.4184×N-Me+0.2803× Facial angle+0.4776×Ans-U1, N-Go=8.6842+0.4099×N-Me+0.4033×S-Go。
2.4 逐步回归与遗传算法预测结果采用遗传算法建立的方程预测值与实测数据进行比较差异均无统计学意义(P>0.05),表明遗传算法建立的方程能很好地显示其定量规律性,精度高。见表 2。
(n=5, x±s) | |||
Group | Measured value | Predicted value | P |
High-angle | |||
Age(year) | 13.40±2.97 | 13.55±1.18 | 0.184 |
N-Me(l/mm) | 107.30±10.70 | 110.25±8.20 | 0.401 |
Ans-U1(l/mm) | 28.49±4.15 | 29.07±2.32 | 0.262 |
N-Go(l/mm) | 69.60±10.69 | 70.26±10.31 | 0.926 |
LI-Me(l/mm) | 35.59±3.52 | 34.27±2.54 | 0.926 |
Average-angle | |||
Age(year) | 13.00±3.46 | 12.88±2.99 | 0.220 |
Ans-U1(l/mm) | 26.14±1.15 | 26.76±1.07 | 0.758 |
Ans-Me(l/mm) | 56.92±2.06 | 57.06±1.95 | 0.885 |
N-Go(l/mm) | 76.21±4.47 | 76.03±4.34 | 0.970 |
S-Ar(l/mm) | 30.37±3.11 | 29.98±2.68 | 0.804 |
Low-angle | |||
Age(year) | 15.2±2.77 | 12.39±1.37 | 0.136 |
N-Me(l/mm) | 108.48±6.72 | 107.28±6.74 | 1.000 |
N-Ans(l/mm) | 50.24±4.71 | 54.01±4.41 | 0.284 |
Ans-Me(l/mm) | 58.24±5.32 | 57.83±4.24 | 0.643 |
N-Go(l/mm) | 82.27±4.42 | 81.76±4.34 | 0.979 |
本研究选取10~18岁就诊患者,该年龄段患者存在一个完整的颌骨生长迸发期[7],利于观察在整个生长迸发期内骨性Ⅱ类错牙合患者的垂直向变化。本研究结果显示:从颅面部垂直向骨骼特征分析,不同骨面型大部分颌面部硬组织测量参数比较差异有统计学意义。高角型病例和低角型病例的形成主要与前、后面部高度的生长发育失调有关。高角组患者下颌角大,下颌平面顺时针旋转,使其前面高增大并伴有开牙合趋势,前面高的增大主要表现在前下面高的增大;而低角组患者下颌升支高度明显高于高角和均角组,同时下颌的逆时针旋转使其后面高明显高于其他组,使这类患者有深覆牙合趋势,并且下颌后缩明显,与有关骨性Ⅱ类错牙合畸形患者垂直向方面的研究[2, 8-10]比较可以看出:无论是横向还是纵向研究,无论是锥形束CT(CBCT)还是头颅侧位片,均可得到相似结果。但Karlsen等[11]通过对6~15岁不同垂直面型的受试者颅面生长变化的研究发现:无论何种垂直骨面型,下颌均向前上旋转,只有几个病例发生后下旋转,但这种现象并不始终存在,故认为高角是受试者的固有特性,而不是下颌骨旋转所导致的。但上述研究并未对受试者进行矢状向分类(安氏Ⅰ、Ⅱ、Ⅲ类)以排除面部生长不同的影响,而是进行纵向研究,与本研究的方法及样本选择均有不同。本研究结果显示:不同性别骨性Ⅱ类错牙合畸形患者垂直向各测量参数比较差异无统计学意义。但Sinclair等[12]的结论与本研究结论不同,可能是因为其研究为纵向研究,且选择的研究样本的种族不同,年龄段及垂直向标志点与本研究也不完全相同。
近年来遗传算法在医学文献检索[13]、图像分割[14]和牙冠宽度预测[15]等医学领域有很多新的研究进展。本研究采用遗传算法得到骨性Ⅱ类错牙合患者垂直向的预测方程,直观表达了不同垂直类型骨性Ⅱ类错牙合患者颅面部垂直向各标志点的关系。以年龄方程为例,由结果可以看出:高角组患者与年龄最相关的为上颌切牙区的高度,且随年龄增大而增大;均角组患者年龄与前面高呈正相关关系;低角组患者年龄与后面高呈正相关关系,与下颌角呈负相关关系。这使得骨性低角患者随年龄增长,有后面高增加、下颌逆时针旋转和覆牙合增大的趋势。以本研究中高角组的预测方程为例可以看出:高角患者下颌骨顺时针旋转、面角减小、下颌角增大,下颌骨的切牙区和上颌骨的切牙区高度均增大,最终导致前面高增大,而导致其增大的主要原因是前下面高的增大;从均角和低角组中也可看出均角和低角患者的垂直向生长变化。所以,遗传算法得到的预测方程可以反映各种垂直类型的生长趋势,与Chung等[2]的纵向研究结果相似。
对于有关生长发育的研究来说,纵向研究的资料收集困难,观察时间长;横向研究虽简单、快捷,但研究结果会隐藏生长变异的细节。如果在横向研究中收集大量生长(骨龄)更加平均分布的样本,结合遗传算法进行优化,其结果也可对生长发育进行预测。
[1] | Chung CH, Wong WW. Craniofacial growth in untreated skeletal Class Ⅱ subjects:a longitudinal study[J]. Am J Orthod Dentofacial Orthop, 2002, 122(6): 619–626. DOI:10.1067/mod.2002.129195 |
[2] | Stahl F, Bacccetti T, Franchi L, et al. Longitudinal growth changes in untreated subjects with Class ⅡDivision 1 malocclusion[J]. Am J Orthod Dentofacial Orthop, 2008, 134(1): 125–137. DOI:10.1016/j.ajodo.2006.06.028 |
[3] | 苑希民, 李鸿雁, 刘树坤, 等. 神经网络和遗传算法在水科学领域的应用[M]. 北京: 中国水利水电出版社,2002: 29-74. |
[4] | 李鸿雁, 赵娟, 王玉新, 等. 扩域搜索遗传算法优化马斯京根参数及其应用[J]. 吉林大学学报:地球科学版, 2011, 41(3): 861–865. |
[5] | 陈顺立. 农业机器人路径优化及轨迹跟踪研究-基于遗传算法[J]. 农机化研究, 2017, 39(8): 17–21. |
[6] | 李纯思, 衣淑娟. 水稻浸种催芽箱温度传感器优化配置——基于遗传算法[J]. 农机化研究, 2017, 39(6): 27–32. |
[7] | Proffit WR, Fields Jr HW, Sarver DW. Contemporary orthodontics[M]. 5th.ed. St.Louis:Moshy Year Book Inc., 2012:90-111. |
[8] | 叶庆, 赵志河, 赵美英. 安氏Ⅱ1类错(牙合)不同骨面型的牙颌垂直向形态特征分析[J]. 上海口腔医学, 2006, 15(4): 370–374. |
[9] | 赵营, 王建国, 魏志强. 骨性Ⅱ类错(牙合)畸形患者不同垂直骨面型下颌骨三维形态特征研究[J]. 天津医科大学学报, 2015, 21(4): 342–344. |
[10] | Karlsen AT. Craniofacial growth differences between low and high MP-SN angle males:a longitudinal study[J]. Angle Orthod, 1995, 65(5): 341–350. |
[11] | Karlsen AT. Association between facial height development and mandibular growth rotation in low and high MP -SN angle face:a longitudinal study[J]. Angle Orthod, 1997, 67(2): 103–110. |
[12] | Sinclair PM, Little RM. Dentofacial maturationof untreated normals[J]. Am J Orthod Dentofacial Orthop, 1985, 88(2): 146–156. DOI:10.1016/0002-9416(85)90239-8 |
[13] | 张晗, 韩爽, 白星, 等. 利用遗传算法确定医学文献的研究热点[J]. 情报分析与研究, 2011, 3: 57–61. |
[14] | 许良凤, 林辉, 罗珣, 等. 基于模拟退火并行遗传算法的Otsu双阈值医学图像分割[J]. 工程图学学报, 2011, 32(2): 88–92. |
[15] | 董妍, 杨陆一, 朱惠芳, 等. 遗传算法预测不同错(牙合)类型牙冠宽度[J]. 口腔医学研究, 2015, 31(7): 728–731. |