AstroBlog

宇宙系大学院生の戯言

Rでのt検定(片側検定、両側検定)のやり方

SPONSORED LINK

今回は自分のためのメモ的な投稿です。

たとえば、
x = -1.22, -1.17, 0.93, -0.58, -1.14
y = 1.03 , -1.59, -0.41, 0.71, 2.10
という2つのデータがあったとします。まずはデータセットをします。

> x <- c(-1.22, -1.17, 0.93, -0.58, -1.14)
> y <- c(1.03 , -1.59, -0.41, 0.71, 2.10)

この2つの標本についてt検定を行うと、

片側検定は

> t.test(x, y, var.equal = FALSE)

    Welch Two Sample t-test

# data:  x and y
# t = -1.3339, df = 6.8431, p-value = 0.2249
# alternative hypothesis: true difference in means is not equal to 0
# 95 percent confidence interval:
#  -2.79205  0.78405
# sample estimates:
# mean of x mean of y 
#   -0.636     0.368 

となり、続いて両側検定は

> t.test(x, y, var.equal = FALSE, alternative="less")

    Welch Two Sample t-test

# data:  x and y
# t = -1.3339, df = 6.8431, p-value = 0.1125
# alternative hypothesis: true difference in means is less than 0
# 95 percent confidence interval:
#       -Inf 0.4268973
# sample estimates:
# mean of x mean of y 
#    -0.636     0.368 

となります。