两样本T检验,和终于来到T检验的最后一个章节,两样本T检验。
两样本T检验的应用条件为:1.独立的随机样本 2.资料应当服从正态分布 3.方差齐性
即我们要在前面两种T检验的前提下,做多一个方差齐性检验,且需要注意数据的结构
废话不多说,do it
例
例3-3为研究国产四类新药阿卡波糖胶囊的降血糖效果,某医院用40名Ⅱ型糖尿病病人进行同期随机对照试验。试验者将这些病人随机等分到试验组(阿卡波糖胶囊组)和对照组(拜唐苹胶囊组),分别测得试验开始前和8周后的空腹血糖,算得空腹血糖下降值如表3-2,能否认为该国产四类新药阿卡波糖胶囊与拜唐苹胶囊对空腹血糖的降糖效果不同?
1.SAS
data EX3_5;
do c = 1 to 2;
do i = 1 to 20;
input x@@;
output;
end;
end;
cards;
-0.7 -5.6 2 2.8 0.7 3.5 4 5.8 7.1 -0.5 2.5 -1.6 1.7 3 0.4 4.5 4.6 2.5 6 -1.4
3.7 6.5 5 5.2 0.8 0.2 0.6 3.4 6.6 -1.1 6 3.8 2 1.6 2 2.2 1.2 3.1 1.7 -2
;
proc univariate normal;
var x;
by c;
proc ttest cocharn;
var x;
class c;
run;
sas语法就不累述了,大部分人都不会选择sas吧!!!
来吧,我们看结果。
C1的正态性检验,即实验组的正态性检验
C2的正态性检验,即对照组的正态性检验
认为都服从正态分布
我们看ttest下的结果,先看方差等价,这个例子中的方差认为的齐的。也就是完美的符合T检验的标准。要是不齐,怎么办呢我们看到上面有‘方差’等于则认为是齐性的。若方差不齐则选择方差不等于的结果。
在这道题中,我们看汇总的结果。认为该国产四类新药阿卡波糖胶囊与拜唐苹胶囊对空腹血糖的降糖效果没有不同。
2.SPSS
sas好像太难了,那我们用spss试试看吧!
录入数据
同意我们先做一个正态性检验。如不知如何点进这一步可以参考我之前写过的文章。正态性检验
查看结果
认为都是服从正态分布的。这里要注意我们不是计算两组合在一起的正态性检验。我们是分别求出两组的正态性检验。
我们进行分析,选到独立样本T检验。
选入变量
点击确定,查看结果
先看方差齐性检验,这里是方差齐的,所以我们选择第一行。SIG=0.525>0.05。该国产四类新药阿卡波糖胶囊与拜唐苹胶囊对空腹血糖的降糖效果没有不同。
3.R语言
搞了好久,才发现自己的版本过低,然后又搞错了数据结构。所以个人觉得想好R语言,必须要有耐心。数据结构很重要。
导入car包为了后面做方差齐性检验。个人习惯开开头导入所有需要的包。
结果和上述的方法一样,不一一赘述。如果有R方面的问题可以直接百度解决。
个人感觉用起来R给我的感觉还是很好的。
4.python
ttest_ind(),equal_car参数默认为TRUE。此例中可以不填写。