A群から取ったデータと、B群から取ったデータの比の信頼区間を求めたいんだけど…
という相談を受けたので探してみると、3種類あるらしい。
Rですぐできそう、というか既にあったのはmratiosパッケージ。
そもそも平均で差がなければ比は1、差があれば1でないから、もともとはt検定を使っているのだが、比の信頼区間を求めてほしいとなるとFiellerという人がごちゃごちゃやっているのでこれを使う。
set.seed(12345) #真の平均の差は1、比は2倍違う。 a1 <- rnorm(5, 1, 1) a2 <- rnorm(5, 2, 1)
t.test.ratio(a1, a2) Ratio t-test for unequal variances data: x and y t = -0.4263, df = 6.363, p-value = 0.684 alternative hypothesis: true ratio of means is not equal to 1 95 percent confidence interval: 0.4056763 2.6170008 sample estimates: mean x mean y x/y 1.2676164 1.4664953 0.8643849
今回はたまたまサンプル数が少ないので、真の平均の比が2倍違っても、1倍ではないとは言えない結果になった。
t検定で同じp値が出る。
t.test(a1, a2) Welch Two Sample t-test data: a1 and a2 t = -0.4263, df = 6.363, p-value = 0.684 alternative hypothesis: true difference in means is not equal to 0 95 percent confidence interval: -1.3249365 0.9271785 sample estimates: mean of x mean of y 1.267616 1.466495
