質問: マイナスのt値 マイナスのt値 実験の結果を見ていると、t値をマイナスで表記しているものを見ることがあります。 実際に検定をする際には絶対値として使用するから実質的には問題ないのでしょうが、そもそもマイナスがでるような計算をすること自体に強烈な違和感があります。 結果の部分では差があるかどうかだけを問題にすればよいので、必ず平均値が大きい方から小さい方を引くものだと考えていました。 統計ソフトの出力がそうなっているといった結果論ではなく、もう少し理論的に、必ず正の値になるようにするべきだ、いや両方使っていい、いやいやそこは絶対値を書くべきだ、といった説明が知りたいです。 あるいは、このように習ったという例でも結構です。 よろしくお願いします。
回答: > もう少し理論的に、必ず正の値になるようにするべきだ、いや両方使っていい、いやいやそこは絶対値を書くべきだ、といった説明が知りたいです。 数学的にt分布において、いわゆる検定統計量であるt値はマイナスの値をとりうるので、マイナスの値がt値として得られることに何ら問題はありませんし、そもそもそれが"理論的"でしょう。 私は逆に > マイナスがでるような計算をすること自体に強烈な違和感があります。 というのが不思議だなぁと思いました(^_^;) > 統計ソフトの出力がそうなっている これはプログラマーがどのような式を用いているか(どのように計算させるか)によるものです。ごく一般的な公式としては分子に (x1_bar - x2_bar) という計算式がでてきますから、プログラム中にこの式を記せば、関数の引数にx1のデータベクトル、x2のデータベクトルを指定したとき、x1がx2より大きければプラスのt値が計算されるでしょうし、逆にx1がx2よりも大きければマイナスのt値が得られるでしょう。 rのt.test()でもそのようになります。 > a <- rnorm(10, mean=100) # 平均が100の正規乱数を10個用意する > b <- rnorm(10, mean=50) # 平均が50の正規乱数を10個用意する > t.test(a, b) # aを先に指定するx1_barは変数aの平均値になる welch two sample t-test data: a and b t = 112.9208, df = 13.99, p-value < 2.2e-16 # t値は正の値 > t.test(b, a) # bを先に指定するx1_barは変数bの平均値になる welch two sample t-test data: b and a t = -112.9208, df = 13.99, p-value < 2.2e-16 # t値は負の値 もちろん、質問者さんがご指摘されているように、実質的には符号の正負は関係ないので abs(x1_bar - x2_bar) といったように、計算式に絶対値として書いてしまう場合だってありえるわけです。