铜仁建设厅官方网站,南京的网站建设,现在网站开发哪个语言好,网上房地产上海一、案例介绍
由单因素方差分析案例中#xff0c;为研究郁金对低张性缺氧小鼠存活时间的影响#xff0c;将36只小鼠随机生成A、B以及 C 三组#xff0c;每组12个#xff0c;雌雄各半#xff0c;分别以10g/kg、20g/kg、40g/kg三种不同剂量的郁金灌胃#xff0c;各组小鼠均…一、案例介绍
由单因素方差分析案例中为研究郁金对低张性缺氧小鼠存活时间的影响将36只小鼠随机生成A、B以及 C 三组每组12个雌雄各半分别以10g/kg、20g/kg、40g/kg三种不同剂量的郁金灌胃各组小鼠均同时置于放有钠石灰的250ml密闭广口瓶中观察并记录小鼠存活时间。想要研究不同剂量的郁金下的小鼠的存活时间是否不同经过分析后发现模型是有显著性差异说明对个总体均数不全相等但是究竟是A、B、C三组总体均数全部不等还是其中两个总体均数不等需要进一步研究案例数据来源于颜虹, 徐勇勇. 医学统计学[J]. 人民卫生出版社, 2015.。
二、问题分析
经过单因素方差分析后发现模型是有显著性差异说明对个总体均数不全相等但是究竟是甲、乙、丙三组总体均数全部不等还是其中两个总体均数不等需要进一步对多个均数做两两比较用多重比较方法。
三、软件操作及结果解读
(一) 数据导入
1.数据格式
首先将数据整理成正确的格式一般X一列Y为一例并且分析的数据带有数据标签的需要另添加一个表格进行说明数据格式如下 2.导入数据
将整理好的数据上传至SPSSAU系统内如下 上传结果如下 (二)事后多重比较分析
对于单因素方差分析过程可以参考前面所发的文章这里对于其相关过程不再赘述。
软件操作 单因素方差分析后进行两两比较进一步研究分析路径为点击【通用方法】→【方差分析】然后进行分析 事后多重比较方法说明
多重比较有很多种方法一般根据他的发展所采用的概率分布和发展过程可以发现一般包括t检验和q检验衍生的方法目前这两种方法的应用较多另外还有一些基于F分布的方法具体说明如下
1LSD方法
LSD是最早的多重比较方法LSD与独立样本t检验非常相近主要差别在于LSD法在首先满足F检验达到显著的基础上将F检验的误差均方作为合并方差LSD法因其计算简单检验效能高所以被广泛应用有研究者模拟试验比较了不同的事后多重比较方法发现在方差中F检验显著时LSD方法是检验效率最高的多重比较方法但是LSD法也存在明显的不足比如涉及过多的要比较均数对犯I型错误的概率较高。其计算公式如下 tα/2为t分布的临界值通过查t分布表得到其自由度为n-kn为样本总数k为因素中不同水平的水平个数MSE为组内方差ni和nj分别为第i个样本和j个样本的样本量。 4.根据显著性水平α做作出决策如果均值之差的绝对值大于LSD则拒绝H0否则不拒绝H0。
2Scheffe
Scheffe检验设定所有可能线性组合Scheffe检验的临界值是用组间自由度乘以F检验的临界值该方法对于非均衡设计更有效但是Scheffe计算相较其他方法更为复杂计算公式如下 3Bonferroni校正
Bonferroni不等式是指一个或多个时间发生的总概率不高于这些事件各自发生概率的加和。但一般认为该方法过于保守更难得到拒绝无效假设的结论Bonferroni法适用于比较次数不是太多的多重比较比较次数越多越保守。计算公式如下 4sidak
Sidak基本思路与Bonferroni法接近一般sidak优于Bonferroni法它使用Bonferroni修正中的一个特例来估算所有假设测试之间存在的共同显著水平。
5Tamhane T2方差不齐
如果方差不齐,但希望进行多重比较则使用此方法其使用前提是需要数据满足正态性但是不需要满足方差齐性。
6SNK Q检验
SNK法的基本目标是把处理均数分成不同处理的子集。这些子集内部都是同质的但是不同子集之间却存在差异因此SNK又被称为分层检验就是按照范围的不同采用不同的检验标准来判定显著性。因为其统计量为q又称为q检验。其统计量的计算如下 式中分子为任意两个对比组A、B的样本均数之差值分母是差值的标准误式 中和中和 分别为A和B两个样本的例数MS_误差为前述方差分析中算得的误差均方。
7Duncan检验
Duncan氏新复极差法目前更多的统计学文献不推荐使用该方法原因在于Duncan多重比较检验误差率的推导是基于单调性条件的单调上升或单调下降且没有震荡这也是很多动物试验结果应用Duncan氏新复极差法比较多的原因。
3.进行两两比较
本例子以SNK Q检验进行两两比较。结果如下
方差分析表格 SNK Q检验如下: 平均值 A组和B组比较时平均值就分别是A组数据和B组数据的平均值。其余同理。差值绝对值
I平均值-J平均值的绝对值差值绝对值比如|40.083-52.985|12.875,其余以此类推。
3K
Q检验的界值不但要考虑自由度还需要考虑秩次R的差别这里用K进行表示 而A组、B组C组的秩如下 比如A组和B组K|1-2|12以此类推。
4)df
df为自由度这里的自由度为误差自由度一共有3组数据每组12个样本所以总的自由度为12×3-125自变量的自由度为3-13所以误差自由度为35-233
55和1
根据K以及自由度进行查看q检验界值表确定p值。
6se
Se为差值的标准误计算如下 7LSR值
LSR值Q临界值*se比如2.887*1.6654.792以此类推
四、结论
经过单因素方差分析后发现模型是有显著性差异说明对个总体均数不全相等经过事后多重比较的SNK Q检验的方法发现A组和B组在0.01水平上呈现显著性说明两者具有差异并且A组和C组B组和C组都是在0.01水平上呈现显著性说明三者之间两两比较均具有显著性差异。
五、知识小贴士
1、如何做Dunnett法事后多重比较
如果需要进行Dunnett事后检验法可使用SPSSAU通用方法里面的非参数检验并且选择多重比较方法为Dunn’s t法即可。
2、提示‘数据质量异常’如何解决
如果是X的某个组别下Y的个数小于2个此时会出现异常情况。建议使用分类汇总进行检查确认此种异常情况后使用筛选样本功能处理后再次分析即可。
