【統計】対応のあるt検定

対応のあるt検定のメモ。


帰無仮説、対立仮設


帰無仮説  H_{0} :前後に差がない( E\left( d\right) = 0
・対立仮設  H_{1} :前後に差がある( E\left( d\right) \neq 0


つまり、関心のある時点での変化量がゼロかどうか。

計算式


 t = \dfrac{E\left( d\right) - 0}{\sqrt{V\left( d\right) /n}} ※d=前後の差


下記は架空データでの計算。前後の差dの平均・分散とn数からt値を求める。


Rでも検算。

> ads <- data.frame(before=c(10,12,8,10,11), after=c(16,15,10,12,16))
> t.test(ads$before, ads$after, paired=T)

    Paired t-test

data:  ads$before and ads$after
t = -4.4313, df = 4, p-value = 0.01141
alternative hypothesis: true difference in means is not equal to 0
95 percent confidence interval:
 -5.855595 -1.344405
sample estimates:
mean of the differences 
                   -3.6 


前後の差を見ているので、引数に変化量とゼロを与えてあげても等価。

> ads$chg <- ads$before - ads$after
> ads$x <- 0
> t.test(ads$chg, ads$x, paired=T)

    Paired t-test

data:  ads$chg and ads$x
t = -4.4313, df = 4, p-value = 0.01141
alternative hypothesis: true difference in means is not equal to 0
95 percent confidence interval:
 -5.855595 -1.344405
sample estimates:
mean of the differences 
                   -3.6 


プログラムコード


R

t.test(before, after, paired=T)


SAS

proc ttest data=ADS;
    paired before * after;
run;


参考


Paired Sample T-Test - Statistics Solutions

本ブログは個人メモです。 本ブログの内容によって生じた損害等の一切の責任を負いかねますのでご了承ください。