采用对数线性模型分析靳三针疗法治疗儿童自闭症的临床疗效资料

郎建英，袁 青，徐 谦，陈丽霞，柯丽萍，俞裕天，吴至凤

(1.广州中医药大学预防医学教研室，广州 510405；2.广州中医药大学针灸康复临床医学院，广州 510405；3.广州中医药大学教务处，广州 510405；4.第三军医大学新桥医院儿科康复中心，重庆 400037)

采用对数线性模型分析靳三针疗法治疗儿童自闭症的临床疗效资料

郎建英1，袁 青2△，徐 谦1，陈丽霞1，柯丽萍3，俞裕天2，吴至凤4

(1.广州中医药大学预防医学教研室，广州 510405；2.广州中医药大学针灸康复临床医学院，广州 510405；3.广州中医药大学教务处，广州 510405；4.第三军医大学新桥医院儿科康复中心，重庆 400037)

目的 探讨对数线性模型在分析靳三针疗法治疗儿童自闭症临床资料中的应用。方法采用对数线性模型分析靳三针疗法与对照疗法治疗218例儿童自闭症的临床资料，在考虑相关因素的情况下比较靳三针针刺组(治疗组)和行为干预组(对照组)疗效。结果当病情程度为重度时，组别与疗效的正相关有统计学意义(P<0.05)，治疗组疗效优于对照组。结论用对数线性模型分析临床资料，可获得更丰富的信息，为指导临床实践提供参考。

针刺疗法；对数线性模型；列联表；儿童自闭症；靳三针

儿童自闭症(children autism)又称为孤独症，是以不同程度的社会交往能力缺陷、语言发育障碍、刻板的行为方式，以及狭隘的兴趣等为主要症状的一种疾病。靳三针疗法治疗儿童自闭症取得了满意的效果[1-3]。以往临床资料均采用常规χ2检验等统计方法比较不同疗法的疗效，较少用对数线性模型分析临床资料。本文对218例儿童自闭症临床资料，用对数线性模型研究相关因素的主效应及其交互效应，从而合理评价比较靳三针疗法和行为干预疗法治疗儿童自闭症的疗效，为指导临床治疗提供参考，研究结果如下。

1 资料与方法

1.1 一般资料 选择2003年9月至2008年12月某三个医疗机构收治的符合入选标准的儿童自闭症患者218例，按家长自愿的原则，分为靳三针针刺组(治疗组)和行为干预组(对照组)。其中治疗组118例，对照组100例。两组患儿性别、年龄、病情程度构成差异均无统计学意义(P>0.05)。

1.2 方法 治疗组采用靳三针疗法，以头部穴组为主；对照组采用行为干预，由儿童特殊训练学校的专业人员对其进行训练。两组均采用每日治疗1次，每周6次，2个月为1个疗程，共治疗2个疗程[4-5]。采用儿童自闭症评定量表(childhood autism rating scale,CARS)进行治疗前后疗效评价，由施术者以外的评价者进行盲法测评，评价者不知道患儿具体分组情况。CARS总分降低5分以上者为有效，降低少于5分者为无效。

1.3 统计学处理 采用SPSS17.0统计软件进行对数线性模型分析，以P<0.05为差异有统计学意义。

2 结 果

2.1 最优对数线性模型的选择及拟合优度检验结果 纳入分析的变量有4个。X：病情程度(1-轻中度；2-重度)；Y：就诊年龄(1-<4岁；2-≥4岁)；Z：组别(1-治疗组；2-对照组)；U：疗效(1-有效；2-无效)。

根据拟合优度检验的G2差值，用后退法选择模型变量，最终获得包含3个变量二阶交互效应的最优模型{XZU}

最优模型的拟合优度检验结果：G2=4.202，df=7，P=0.756。P>0.05，表明模型对资料有很好的拟合效果，可以认为该模型为最优模型。即变量X、Z和U的二阶交互效应存在，且有统计意义，可以认为病情程度不同时，组别与疗效的交互效应不同。

表1 {XZU}模型单元计数和残差(n)

2.2 残差分析 由表1可见，三维列联表所有格子的标准化残差的绝对值均小于1.96，表明模型很好地拟合实际数据。

表2 {XZU}模型参数估计

3 讨 论

本研究数据来自国家中医药管理局中医临床诊疗技术整理与研究项目“靳三针治疗儿童自闭症临床规范化研究”，研究结果已发表多篇文章[4,6-9]。研究结果显示治疗组疗效优于对照组。但该研究设计不属于实验研究，符合入选标准的患儿不是随机分组而是按家长意愿分组。虽然两组儿童在性别、就诊年龄和病情轻重构成上差异无统计学意义，但仍不排除这些混杂因素存在影响。

两分类变量列联表资料的分析，通常采用χ2检验，如果有第3个分类变量为混杂因素，可以采用分层χ2检验。但若有多个混杂分类变量时,χ2检验无法对多个分类变量间的关系作出综合的评价。对数线性模型可以分析各变量的主效应及变量间的交互效应，可在控制混杂因素的同时对变量的效应做出估计。

对数线性模型将列联表的总效应取自然对数后分解成主效应和各因素间的交互效应[10]，通过最优模型来分析因素间的关系[11]。统计学家提出许多选择最优模型的准则和筛选策略，其中以“后退法”逐步筛选变量最常用。

选出最优模型后，根据模型拟合优度检验的P值评价模型对实际数据的拟合效果，P>0.05表示模型对实际资料有好的拟合效果。模型的拟合优度检验从整体上反映资料的拟合情况，而残差分析则可更细致地考察每个格子的拟合情况。一般标准化残差在(±1.96)范围内表明拟合效果好，模型选择合适[12]。根据最优模型中各效应参数估计值的正负和P值，结合实际专业知识对计算结果进行综合解释。当模型中主效应和交互效应的参数估计值大于0时为正效应，小于0时为负效应。

本研究应用对数线性模型分析，探讨4个主要因素(病情程度X、就诊年龄Y、组别Z及疗效U)的主效应及其交互效应，通过“后退化”选择模型变量，结果3个变量(病情程度X、组别Z及疗效U)进入模型,得到包含X 、Z和U二阶交互效应的最优模型{XZU}，拟合优度检验P>0.05，表明病情程度、组别和疗效三变量间存在二阶交互作用。结合专业知识，可以认为病情程度不同，组别与疗效的关系不同。病情程度为轻中度时，组别与疗效为正效应，但无统计学意义(P>0.05)；病情程度为重度时，组别与疗效为正效应，且有统计学意义(P<0.05)，即组别的值越小(1-治疗组；2-对照组)，疗效的值也越小(1-有效；2-无效)。可以认为，对于重型患者，靳三针疗法的疗效优于行为干预疗法。这一结论具有一定的临床指导意义。

本研究结果表明，对数线性模型分析处理临床资料可获得比传统χ2检验更丰富的信息，是挖掘高维列联表临床资料丰富信息的有用工具。

[1]马瑞玲，袁青，靳瑞．针刺配合行为干预疗法对儿童自闭症行为的影响[J]．中国中西医结合杂志，2006，26(5)：419-422．

[2]袁青，柴铁劬，郎建英，等．针刺治疗儿童自闭症40例疗效观察[J]．广州中医药大学学报，2007，24(3)：208-210,214．

[3]Yuan Q,Ma RL,Jin R.Effect of acupuncture on cerebral images in autism children[J].World J Acupu-Mox,2004,14(3):3.

[4]袁青，吴至凤，汪睿超，等.针刺对儿童自闭症不同临床症状影响的临床研究[J].新中医，2009，41(3)：83-85.

[5]袁青，刘祎思，俞裕天，等．头穴留针配合行为训练治疗儿童自闭症疗效观察[J]．中国针灸，2013，33(7)：609-613．

[6]吴至凤，袁青，汪睿超，等．靳三针治疗不同年龄段自闭症儿童疗效观察[J]．重庆医学，2009，38(21)：2685-2687．

[7]袁青，吴至凤，汪睿超，等．靳三针治疗儿童自闭症不同中医证型疗效分析[J]．广州中医药大学学报，2009，26(3)：241-245．

[8]袁青，汪睿超，吴至凤，等．靳三针治疗重度自闭症疗效对照观察[J]．中国针灸，2009，29(3)：177-180．

[9]袁青，吴至凤，汪睿超，等．针刺对不同病情程度儿童自闭症疗效分析[J]．针刺研究，2009，34(3)：183-187．

[10]方积乾．医学统计学与电脑实验[M]．上海：上海科学技术出版社，1997：521-530．

[11]柳青．多元统计分册——中国医学统计百科全书[M]．北京：人民卫生出版社，2004：259-261．

[12]王保进．多变量分析-统计软件与数据分析[M]．北京：北京大学出版社，2007：411-450．

验交流·

10.3969/j.issn.1671-8348.2015.23.033

国家中医药管理局中医临床诊疗技术整理与研究项目(2003ZL 第49号)。

郎建英(1960-)，副教授，博士，主要从事医学统计学在中医药研究中的应用研究。

，E-meil:yuanqing1005@126.com。

R246.4

B

1671-8348(2015)23-3260-03

2015-02-10

2015-07-15)

△