如果两个样本具有方差齐性,那么做独立样本t检验时,直接套用t检验的公式,计算t值,,查表的自由度为n1+n2-2,然后用函数pt( t value, n1+n2-2)给出p值,小于0.05即为显著。
如果方差不齐,需要计算校正后的自由度,
方法如下:
设样本量为 n,假定为 30 , R 里面提供了 var 函数来求样本方差
var 函数的定义是:
但样本方差的定义是
n <- 30x <- rnorm(30)# 样本标准差为print(var(x)*(n/(n-1)))
如果两个样本具有方差齐性,那么做独立样本t检验时,直接套用t检验的公式,计算t值,,查表的自由度为n1+n2-2,然后用函数pt( t value, n1+n2-2)给出p值,小于0.05即为显著。
如果方差不齐,需要计算校正后的自由度,
方法如下:
设样本量为 n,假定为 30 , R 里面提供了 var 函数来求样本方差
var 函数的定义是:
但样本方差的定义是
n <- 30x <- rnorm(30)# 样本标准差为print(var(x)*(n/(n-1)))