サイズの異なる2つのデータに対して対応のあるt検定を行いたい

Question

pythonで対応のあるt検定を行いたいのですが、2つの群のデータ数が異なる場合の対処方法がわかりません。

from scipy import stats
stats.ttest_rel(A, B)

ValueError: unequal length arrays

サンプル数の異なる2つの群に対してt検定を行う方法を教えていただきたいです。
ご回答よろしくお願いします。

Comments

[magichan]

サンプル数の異なる段階で、対応はないのでは？

[pdaocmhii]

ご回答ありがとうございます。
勉強不足でした。私が示したかったものがそもそも対応のあるt検定ではなかったです。
A群には、そのまま何もせず、B群には薬を投与して、1ヶ月後のA群、B群の生体情報に差が見られるかというものを示したかったのですが...。A群とB群では被験者の数が異なっているので検定が行えず、悩んでいました。
統計に関する知識が乏しく申し訳ありません。

Answer 1

解決済みですが、コメントします。
対応のある検定は同一群に対して介入をした場合としない場合を比較することで、介入の効果を評価する時に使います。今回の場合は群が異なるので対応のない検定を使うほうが妥当と思います。対応のない検定の場合、分散が同じかどうかでやり方が違うようですが、jeanbiegoさんの紹介する記事の通り、ウェルチのt検定を使っておけば問題ないでしょう。個人的にはt検定のような解析的アプローチよりモンテカルロシミュレーションとかを使った数値計算的アプローチのほうがいまどきぽくていいと思います

Answer 2

なぜ対応のあるt検定なのにデータ数が異なるのでしょうか。欠損がある等であれば、そのデータセットを除くなど前処理をした方がいいかもしれません。

↓こういう検定がしたいということですよね？
対応のある2標本t検定

Comments

[pdaocmhii]

ご回答ありがとうございます。
勉強不足でした。私が示したかったものがそもそも対応のあるt検定ではなかったです。
A群には、そのまま何もせず、B群には薬を投与して、1ヶ月後のA群、B群の生体情報に差が見られるかというものを示したかったのですが...。A群とB群では被験者の数が異なっているので検定が行えず、悩んでいました。
統計に関する知識が乏しく申し訳ありません。

[jeanbiego]

対応のないt検定をscipyで紹介している例がありました。
ご参考まで。
https://qiita.com/ynakayama/items/b9ec31a296de48e62863

[pdaocmhii]

ありがとうございます！