昨天讲过了两个样本均数比较的t检验的SPSS教程,今天我们再和大家分享一下SAS操作,文末还有SAS与SPSS的一个简单比较,感兴趣的不要错过哦。
一、问题与数据
随机抽取20只小鼠分配到A、B两个不同饲料组。每组 10 只,喂养一定时间后,测得鼠肝中铁的含量(μg/g),数据如下。试问不同饲料对鼠肝中铁的含量有无影响?
二、对数据结构的分析
整个数据资料涉及2组,每组随机抽取10个数据,测量指标为鼠肝中铁的含量,因此属于两组设计的定量资料。
要知道不同饲料对鼠肝中铁的含量有无影响,则要比较2组间的总体均数差值是否有统计学差异。若数据满足独立性、样本均数服从正态分布、两个总体方差齐,可选用两个样本均数比较的t检验。
三、SAS分析方法
1. 数据录入SAS
2. 使用proc ttest命令→将Group定义为分组变量(class),Fe定义为反应变量(var)→run。
四、结果解读
1. 前2个表格给出了两组鼠肝中铁含量的部分统计信息,包括组别(Group)、数量(N)、均值及其95%置信区间、标准差及其95%置信区间、样本均数的标准误、最小值和最大值,以及两组差值的上述指标。
2. 方差等价表格给出了两组数据方差齐性检验的结果,F值=3.74,P(Pr > F)=0.0627 < 0.10(方差齐性检验一般α选0.10,为了减小犯II型错误的概率),说明两组数据方差不齐。
3. 由于两组数据方差不齐,因此看方差“不等于”的结果(即方法为Satterthwaite一行),t值=2.11,P(Pr > |t|)=0.0541。(若方差齐,则使用方差“等于”一行的结果)
五、撰写结论
A饲料组鼠肝中铁含量为(2.40 ± 1.09)μg/g,B饲料组鼠肝中铁含量为(1.58 ± 0.56)μg/g。尚不能认为两组饲料对鼠肝中铁含量影响的差异有统计学意义(t=2.11,P=0.054)。
六、延伸阅读
SAS软件与SPSS软件关于方差齐性检验的计算:
SAS的计算公式为:
S1和S2分别为两样本的标准差,如果H0(即两总体方差相同)成立,则S12/ S22服从F分布。
本例用SAS进行方差齐性检验,得出的F值=3.74,P=0.0627;如果用SPSS进行分析,F值=8.246,P=0.010,二者结果略有差别。SPSS进行方差齐性检验时,采用的方法为Levene’s检验,可以检验多组方差是否相等。实际上Levene’s检验,就是用测量值减去各组测量值的均数取绝对值后做单因素方差分析的结果。
(如果你想使用文中数据进行练习,请随时给小咖(微信:xys2016ykf)发消息,小咖将原始数据发给你。)
研究进展精彩回顾
研究方法精彩回顾
医咖会微信:medieco-ykh
我们建了一个微信群,有临床研究设计或统计学方面的难题?快加小咖个人微信(xys2016ykf),拉你进群,和其他小伙伴们一起交流学习吧。
点击左下角“阅读原文”,看看医咖会既往推送了哪些研究方法。还可以到医咖会公众号下方的自定义菜单,点击“直接搜索”,查找你想了解的内容。