统计第12次作业

1、WORD格式可编辑统计学第12次作业 郭晓兰 微生物学122136412.将18名原发性血小板减少症患者按病情和年龄都相近的原则配为6个单位组，每个单位组中的3名患者随机分配到 A、B C三个治疗组中，治疗后患者的血小板升高，结果见教材表16-16，问3种治疗方法的疗效有无差别？教材表16-16不同人用鹿茸草后血小板的升高值/ (X 1012 L1)年龄组ABC24.66.311.937.69.214.766.26.913.4分析：该题为随机区组设计，应用随机区组方差分析。总变异可分解为处理变异、区组变异和误差。解：（1 ）统计

2、描述表1 18名原发性血小板减少症患者用不同治疗方法后血小板升高值情况TreatBlockNMea nStd.Deviatio nA1666.2001.793B1667.8331.619C16612.5172.411Total18188.8503.319（2 ）统计推断建立检验假设，确定检验水准：对于处理组：H 0 :三种治疗方法的平均血小板升高值相等。H1 :三种治疗方法的平均血小板升高值不全相等。对于区组：Ho :各区组的总体效应相等。H i :各区组的总体效应不全相等。均取:-=0.05;假设资料符合方差分析的前提条件，用SPS软件进行随机区组方差分析，步骤如下：An alyze 宀 G

3、en eral Lin ear ModelUni variate Depe ndents variable :血小板升高值；Fixed factors : treat ; Block ; Model: Custom： Build terms : main effects ; Treat : Model; Block : Model; Continue 点击 Post Hoc , Treat: Post Hoc Test for ; EqualVariances Assumed: LSD 、 SNK Bonferroni ; Continue 点击 Save, Unstandardized (P

4、redicted Values) , Standardized (Residuals)OK将数据进行方差分析后，计算结果如表2:表2方差分析表变异来源SSVMSFP总变异187.26517处理间129.003264.50279.3380.000 0区组间50.132510.02612.3330.000 1误差8.13100.813由表可知：对处理因素和区组因素统计检验后的P均小于0.05，按=0.05水准拒绝H 0接受H1，可认为三种治疗方法的平均血小板升高值不全相等；可认为：对于同一种治疗方法而言，各区组血小板升高值的差异有统计学意义。Bonferroni法对各处理组进行两两比较建立检验假设

5、，确定检验水准：H ° :叫二丄2，即相互比较的两种治疗方法血小板升高值相等。H 1 :叫=J2，即相互比较的两种治疗方法血小板升高情况不相等。 -0.0167。用SPSS软件进行多重比较后得出表3:表3三种治疗方法血小板升高值的多重比较（ Bonferroni 法）对比组1与2£ - X2t值P值A与B治疗组1.6333.3170.032A与C治疗组6.31712.1340.000B与C治疗组4.6838.9960.000由于表3的P值是由SPSS软件计算而来，该 P值与=0.05检验水准比较，各个对比组的P值均小于0.05，拒绝H °，接受H i，可认为B治疗

6、方法血小板升高值高于A治疗方法，两者平均血小板升高值差值的95喏信区间为(0.139 , 3.127 ) (x g12 l-1 );可认为C治疗方法血小板升高值高于A治疗方法，两者平均血小板升高值差值的95酒信区间为(4.823 ,7.811 ) (x 1012 L1);可认为C治疗方法血小板升高值高于B治疗方法，两者平均血小板升高值差值的 95%置信区间为(3.189 , 6.177 ) (x 1012 l-1 )(3 )残差图判断资料是否满足方差分析前提条件绘制残差图：Graphs Scatter Simple 宀 Define: Y Axis: ZRE_1, X Axis:treat O

7、KGraphs Scatter Simple Define: Y Axis: ZRE_1, X Axis:block OKGraphs Scatter Simple Define: Y Axis: ZRE_1, X Axis: pre_1OK_ J. _ m 佃 脑In sks!1na.udi!z:niQxlrlLIia匚日-吋口曲載 n 盲第申曾金丘d IRa*c.u*右5 水«卄百协ac>2 juTawta帖血Predicted Value fur k4KAtreat与残差的残差图图2：预测值与残差的残差图sun&niirg.?«.dua 二 or 血小戦

8、斤m 们t.DD0 DQ-1 0D-£410blackz.oo图3 :区组与标准残差的残差图三个残差图显示无特殊值，基本满足方差分析要求的方差齐性和服从正态分布的条件，可以用方差分析对资料进行分析。4.某研究人员以0.3 ml/kg剂量纯苯给大鼠皮下注射染毒，每周3次，经45天后，使实验动物白细胞总数下降至染毒前的50%左右，同时设置未染毒组。两组大鼠均按照是否给予升高白细胞药物分为给药组和不给药组，实验结果见教材表16-18，试作统计分析。教材表16-18 实验效应指标（吞噬指数）数据未染毒组染毒组不给药给药不给药给药3.803.881.851.943.903.842.012.25

9、4.063.962.102.033.853.921.922.103.843.802.042.08解：本资料为析因设计资料，应用析因设计资料的方差分析。（1 ）统计描述表42 2析因设计四种组合的均数和标准差升高白细胞药物纯苯皮下注射染毒X.j未染毒染毒不给药3.890 ±0.1011.984 ±0.0992.937给药3.880 ±0.0632.080 ±0.1132.980Xi.3.8852.0322.959（2 ）统计推断建立检验假设，确定检验水准；H10 : （-："'） =0即染毒组和未染毒组不存在交互效应。H11 : （-：）

10、-'） = 0即染毒组和未染毒组存在交互效应。H20:? =0即染毒因素的主效应为 0。H 21 ： =0即染毒因素的主效应不为0。H30 :1 =0即给药因素的主效应为 0。H31 :卩=0即给药因素的主效应不为0。均取=0.05 ;假设资料符合方差分析的前提条件，用SPS软件进行析因设计资料的方差分析，步骤如下：An alyze Gen eral Lin ear Model Uni variate Depe ndents variable :吞噬指 数；Fixed factors :药物；纯苯染毒；Model: Full , Continue点击 Plots ,Horizontal

11、 Axias:纯苯染毒；Separate Lines:药物；Add , Continue 点击 Save, Unstandardized (Predicted Values) , Standardized (Residuals) , Continue 丿点击 Options , Display means for: 纯苯染毒、药物； Compare main effects ； Display: descriptive statistics Continue OK结果如表5:表5 方差分析表变异来源SSVMSFP总变异17.33919染毒与否17.168117.1681857.0090.000

12、给药与否0.00910.0091.0000.332交互作用0.01410.0141.5190.236误差0.148160.009由表可见，F交互=1.519 , P =0.236, P>0.05 ，接受H ° ,可认为染毒和给药无交互效应。F染毒=1857.009, P<0.05，拒绝H。，接受H i，可认为染毒与不染毒对应的总体均 数不等，即染毒因素存在主效应，主效应为1.853，其95酒信区间为(1.762 ,1.944 )。F给药=1.000 , P = 0.332, P>0.05，接受H。，可认为给药与不给药对应的总体均数相等，即给药因素没有存在主效应。(3 )残差图判断资料是否满足方差分析前提条件绘制残差图：Graphs Scatter Simple Define: Y Axis: ZRE_1, X Axis:给药OKGraphs Scatter Simple Define: Y Axis: ZRE_1, X Axis:染毒OK专业技术分享Graphs S catter Simple 宀Define: Y Axis: ZRE_1, X Axis: pre_1OKQ:Land£J<hMi1ilHll-eL匚里 rGi矗好析竄slEJnabrEL円ftPKniSJa 匚：1二口4 布 離DaDaCi图4给药与标准残差的