1、统计学第 12 次作业 郭晓兰 微生物学 122136412将 18 名原发性血小板减少症患者按病情和年龄都相近的原则配为 6 个单位组，每个单位组中的 3 名患者随机分配到 A、B、C 三个治疗组中，治疗后患者的血小板升高，结果见教材表 16-16，问 3 种治疗方法的疗效有无差别？教材表 16-16 不同人用鹿茸草后血小板的升高值/（ L-1）120年龄组 A B C1 3.8 6.3 8.02 4.6 6.3 11.93 7.6 10.2 14.14 8.6 9.2 14.75 6.4 8.1 13.06 6.2 6.9 13.4分析：该题为随机区组设计，应用随机区组方差分析。总变异可分
2、解为处理变异、区组变异和误差。解：（1）统计描述表 1 18 名原发性血小板减少症患者用不同治疗方法后血小板升高值情况Treat Block N Mean Std.DeviationA 16 6 6.200 1.793B 16 6 7.833 1.619C 16 6 12.517 2.411Total 18 18 8.850 3.319（2）统计推断建立检验假设，确定检验水准：对于处理组：三种治疗方法的平均血小板升高值相等。0H：三种治疗方法的平均血小板升高值不全相等。1对于区组：各区组的总体效应相等。0：各区组的总体效应不全相等。均取 ；1H05.假设资料符合方差分析的前提条件，用SPSS软
3、件进行随机区组方差分析，步骤如下：AnalyzeGeneral Linear ModelUnivariateDependents variable：血小板升高值；Fixed factors：treat；Block；Model：Custom：Build terms：main effects；Treat : Model；Block : Model；Continue点击Post Hoc，Treat: Post Hoc Test for；Equal Variances Assumed: LSD、SNK、Bonferroni；Continue点击Save，Unstandardized (Predicte
4、d Values)，Standardized (Residuals) OK将数据进行方差分析后，计算结果如表2：表2 方差分析表变异来源 SMSFP总变异 187.265 17处理间 129.003 2 64.502 79.338 0.000 0区组间 50.132 5 10.026 12.333 0.000 1误差 8.13 10 0.813由表可知：对处理因素和区组因素统计检验后的 均小于0.05，按 水准拒绝 ，P5.0H接受 ，可认为三种治疗方法的平均血小板升高值不全相等；可认为：对于同一种治疗1H方法而言，各区组血小板升高值的差异有统计学意义。Bonferroni法对各处理组进行两两
5、比较建立检验假设，确定检验水准： ，即相互比较的两种治疗方法血小板升高值相等。0H21： ，即相互比较的两种治疗方法血小板升高情况不相等。67.用 SPSS 软件进行多重比较后得出表 3：表 3 三种治疗方法血小板升高值的多重比较（Bonferroni 法）对比组 1 与 2 21Xt 值 值PA 与 B 治疗组 1.633 3.317 0.032A 与 C 治疗组 6.317 12.134 0.000B 与 C 治疗组 4.683 8.996 0.000由于表 3 的 值是由 SPSS 软件计算而来，该 值与 检验水准比较，各个对比组PP05.的 值均小于 0.05，拒绝 ，接受 ，可认为
6、B 治疗方法血小板升高值高于 A 治疗方0H1法，两者平均血小板升高值差值的 95%置信区间为（0.139，3.127） （ L-1）；可认为 C12治疗方法血小板升高值高于 A 治疗方法，两者平均血小板升高值差值的 95%置信区间为（4.823，7.811） （ L-1）；可认为 C 治疗方法血小板升高值高于 B 治疗方法，两者平均120血小板升高值差值的 95%置信区间为（ 3.189，6.177） （ L-1）。120（3）残差图判断资料是否满足方差分析前提条件绘制残差图：Graphs Scatter SimpleDefine: Y Axis: ZRE_1, X Axis: treatO
7、KGraphs Scatter SimpleDefine: Y Axis: ZRE_1, X Axis: blockOKGraphs Scatter SimpleDefine: Y Axis: ZRE_1, X Axis: pre_1OK图1：treat 与残差的残差图 图2：预测值与残差的残差图图3：区组与标准残差的残差图三个残差图显示无特殊值，基本满足方差分析要求的方差齐性和服从正态分布的条件，可以用方差分析对资料进行分析。4某研究人员以 0.3 ml/kg 剂量纯苯给大鼠皮下注射染毒，每周 3 次，经 45 天后，使实验动物白细胞总数下降至染毒前的 50%左右，同时设置未染毒组。两组大鼠
8、均按照是否给予升高白细胞药物分为给药组和不给药组，实验结果见教材表 16-18，试作统计分析。教材表 16-18 实验效应指标（吞噬指数）数据未染毒组 染毒组不给药 给药 不给药 给药3.80 3.88 1.85 1.943.90 3.84 2.01 2.254.06 3.96 2.10 2.033.85 3.92 1.92 2.103.84 3.80 2.04 2.08解：本资料为析因设计资料，应用析因设计资料的方差分析。（1）统计描述表 4 析因设计四种组合的均数和标准差2纯苯皮下注射染毒升高白细胞药物未染毒 染毒j.X不给药 3.890 0.1011.984 0.0992.937给药 3
9、.880 0.063 2.080 0.113 2.980.iX3.885 2.032 2.959（2）统计推断建立检验假设，确定检验水准； ： 即染毒组和未染毒组不存在交互效应。 10H)(： 即染毒组和未染毒组存在交互效应。 ： 即染毒因素的主效应为0。 20： 即染毒因素的主效应不为0。 1H： 即给药因素的主效应为0。 30： 即给药因素的主效应不为0。 31H均取 ；5.假设资料符合方差分析的前提条件，用SPSS软件进行析因设计资料的方差分析，步骤如下：AnalyzeGeneral Linear ModelUnivariateDependents variable：吞噬指数；Fixed
10、 factors：药物；纯苯染毒；Model：Full ，Continue点击Plots， Horizontal Axias: 纯苯染毒； Separate Lines: 药物；Add ，Continue 点击Save，Unstandardized (Predicted Values) ，Standardized (Residuals)，Continue 点击Options，Display means for: 纯苯染毒、药物； Compare main effects ；Display: descriptive statistics ContinueOK结果如表5：表5 方差分析表变异来源
11、SMSFP总变异 17.339 19染毒与否 17.168 1 17.168 1857.009 0.000 给药与否 0.009 1 0.009 1.000 0.332 交互作用 0.014 1 0.014 1.519 0.236误差 0.148 16 0.009 由表可见， ， , P0.05 ，接受 ，可认为染毒和给药无交互效519.交 互F236.0P0H应。 P0.05 ，接受 ，可认为给药与不给药对应的总体均数相.1给 药F32.P0等，即给药因素没有存在主效应。（3）残差图判断资料是否满足方差分析前提条件绘制残差图：Graphs Scatter SimpleDefine: Y Axis: ZRE_1, X Axis: 给药OKGraphs Scatter SimpleDefine: Y Axis: ZRE_1, X Axis: 染毒 OKGraphs Scatter SimpleDefine: Y Axis: ZRE_1, X Axis: pre_1OK图4 给药与标准残差的残差图 图5 染毒与标准残差的残差图图6 预测值与标准残差的残差图三个残差图显示无特殊值，基本满足方差分析要求的方差齐性和服从正态分布的条件，可以用方差分析对资料进行分析。
