卫生统计学方差分析.pptxVIP

第八章 多组定量或等级资料平均值的比较（P92） ;先来看一个具体的例题;从本例引伸的问题;第一节 完全随机设计多组 均数比较的方差分析 ;表8-1 某受试物不同浓度组高寿命果蝇生存天数;;表8-1 某受试物不同浓度组高寿命果蝇生存天数;? =N-1;三种变异的关系： SS总= SS组间+SS组内 v总= v组间+v组内 ;成组设计方差分析的计算公式;（1） 建立检验假设 H0：多个样本总体均数相等。 H1：多个样本总体均数不相等或不全等。 检验水准为0.05。 （2） 计算检验统计量F值 （3） 确定P值并作出推断结果;表8-1 某受试物不同浓度组高寿命果蝇生存天数;查附表 （方差分析用）F界值表，F3.81,P0.05,可认为4个处理组总体平均生存天数不全等或全不等;表8-2 例8-1的方差分析表;成组设计方差分析数据满足的条件;基本原理;4. 多个均数间的两两比较 当检验结果P≤0.05后，可进一步对多个均数作两两比较。当进行多个均数间的两两全面比较时，采用SNK(Student-Newman-Keuls)检验；当g-1个处理组分别与一个共同的对照组比较，而各处理组间不作两两比较时，采用dunnett-t检验。;;SAS程序;;某一组或多组不服从正态分布或分布类型未知，或各组总体方差不齐，可采用完全随机设计多组分布比较的秩和检验 （Kruskal-Wallis H）检验。 此种非参数检验方法适用面广，但增加了第Ⅱ类错误的概率β，降低了统计检验的功效1－β 。;一、多组连续变量资料的秩和检验 ; 5组资料经方差齐性Levene检验，F=2.3505,P=0.06840.10,可认为方差不齐。;⑴　建立检验假设 H0 :5总体分布位置相同 H1 :5总体分布位置全不相同或不全相同　 　α=0.05 ⑵ 编秩求秩和 将全部数据由小到大按自然数1，2，3，…的顺序及大小统一编秩。当某若干个数相同且位于同一组时按顺序编秩；当某若干个数相同但位于2个及以上组时，将该若干个数按顺序编秩后对秩取均数，此均数即为该若干个数的共同秩次。计算得i(i=1,2,..,5) 组的秩和Ri分别为：167.5，216.5，256.0，299.5，335.5,见表8-2。;表8-2 不同浓度（％）受试物组高寿命果蝇生存天数的比较;⑶计算检验统计量 ;⑷ 确定P值下结论 当组数等于3且每组例数≤5时，可查有关统计书籍附H界值表确定P值。当组数大于3或虽组数等于3但最小样本例数大于5时，H近似服从自由度为组数－1的卡方分布。 本例H=8.289.49=,P0.05,可认为5总体分布位置相同，即5组的平均生存天数相同。;⑸ 两两比较 当结论为“各组总体分布位置全不相同或不全相同”时，需要进一步作两两比较。可采用完全随机设计多组数据秩转换后作方差分析并两两比较的方法，参见本段“2.SAS 程序及结果”中的“⑶ 完全随机设计多组数据秩转换后作方差分析并两两比较的SAS程序”。;也可采用两样本秩和检验的方法，借助SAS或SPSS软件得到相应的P值（参见第七章），按公式α’=α/k计算各次两两比较的检验水平α’，这里α为总检验水平，k为两两比较的次数。;相关SAS程序;二、多组有序变量资料的秩和检验;表8－3不同文化程度组民工对 性自慰/手淫的认识程度等级比较 ;⑴建立检验假设 H0 :4总体分布位置相同 H1 :4总体分布位置全不相同或不全相同　 　α=0.05 ⑵ 编秩求秩和 先计算各等级的合计，再确定各等级的秩次范围及平均秩次，如等级“同意”的合计为778，则其秩次范围为1－778，其平均秩次为（1＋778）/2＝389.5；;又如等级“无所谓”的合计为957，则其秩次范围为779－（778＋957），即779－1735，其平均秩次 为1257。以各等级的平均秩次近似代替该等级各调查对象的秩次，计算各组的秩和，如小学组秩和＝389.5×114＋1257×133＋2322×185＝641154。;⑶ 计算检验统计量 按公式（8－6）计算检验统计量H值，当各样本相同秩较多时，用H除以C的商对H作校正。这里，其中为第j个相同秩次的个数，N为各组例数之和。;⑸ 两两比较 当结论为“各组总体分布位置全不相同或不全相同”时，需要进一步作两两比较。可采用两样本等级资料秩和检验的方法，借助SAS或SPSS软件得到相应的P值（参见第七章），按公式α’=α/k计算各次两两比较的检验水平α’，这里α为