黒バスと弱ペダキャラの身長平均値の差の検定

医学部の友人(@_surrealisme_)に統計学を黒バスに例えて教えてたら調子に乗って実際にやってみた。 概要だけ感覚的に掴めるように説明したので端折った部分、正確でない部分があります。また全力でふざけた結果なのでお見苦しい点もありますがご了承ください。 私自身まだまだ統計学を理解しきれてない部分ばかりですが、もし大学で統計学をかじってて分からないという人がいましたら少しでも親しみを持てるきっかけになれば幸いです。
2
ちえみ⛩ @_surrealisme_

@rikka327 @gojiharuwanwan 片側検定と両側検定の違いがわからないっす! あと有意水準がなんなのかも微妙

2014-07-24 22:43:36
律歌(りっか) @rikka327

@_surrealisme_ 有意水準。例えばある事象が偶然起こる確率が0.01であったとする。これは100回に一回は起こる、珍しい事象だと言える。確率が低ければ低いほど偶然では起こらない。これを意味のある確率、有意確率という。ではどこからが有意確率なのか。(続

2014-07-24 22:56:27
律歌(りっか) @rikka327

@_surrealisme_ 大抵心理学では0.05未満や0.01未満を有意であるとしている。医学系なら0.001もざら。確率がどこから有意であるかの水準、これを有意水準と呼ぶ。

2014-07-24 22:59:19
ちえみ⛩ @_surrealisme_

すげぇ!!!!! 統計が日本語だ!!!! わかる!!!!! わかりやすすぎて泣きそう

2014-07-24 23:00:35
ちえみ⛩ @_surrealisme_

@rikka327 おけーい! すげえ、αとpの違いがよくわかった りっか先生だ……

2014-07-24 23:05:42
律歌(りっか) @rikka327

ちえみに説明するために一切使わなかった統計の教科書開いてんだけど、絶対わからんしょこれ。私がわからんもん(え

2014-07-25 00:02:37
律歌(りっか) @rikka327

統計は実際に使ってみると割としっくりくること多いから一回統計ソフトをどうやって使うか教えて、そのあとに基礎を教えた方がいい気がしてくる

2014-07-25 00:14:16
律歌(りっか) @rikka327

ついでに思いついたので黒バスでt検定の説明でもしとく。

2014-07-25 00:32:09
律歌(りっか) @rikka327

あるオタク「黒バスのキャラの身長は他の漫画に比べて絶対高い!」 んまぁそりゃそうだろうけどこれを統計的に説明しようと思う。具体的にゆうと、黒バスキャラの身長の平均とよわペダのキャラの平均を比較しようと思う。これがよくいうt検定の中でも『対応のない平均値の差の検定』ってやつ。

2014-07-25 00:35:33
律歌(りっか) @rikka327

対応のない差の検定でも色々順序を踏むんだけど。とりあえず黒バス勢の平均値が180でよわペダ勢の平均値が160だとして、あと平方和出したり自由度出してt値出して正規分布のt値と比較して正規分布のそれよりデカかったら有意。とか多分そんなん。もう統計ソフトでやってるからわからんちん

2014-07-25 00:42:58
律歌(りっか) @rikka327

んでだ。これが黒バスとよわペダならt検定でええんやけど同じ身長が有利なHQも比較しようぜって話になった。するとこれは三群やからt検定(二群のみ有効)は出来ん。ので三群以上は分散分析を使うん。

2014-07-25 00:46:53
律歌(りっか) @rikka327

またHQをいれるんじゃなくて、「いや待て、黒バスの中でも女子キャラは普通やろ。男女と作品別比較しようぜ」という話になり男女×作品別の四群となる。これも分散分析。更に要因が二つになる(さっきのは作品が増えただけで要因は一つ)ので二要因の分散分析となる。

2014-07-25 00:49:36
律歌(りっか) @rikka327

二要因にすると何が起きるかと言えば、「黒バスの中でも男子キャラはどの群よりも平均値が高い」という性差×作品別の効果、つまり"二つの要因が互いに及ぼす影響"が検討できる。この影響を専門用語で交互作用という。

2014-07-25 00:52:27
律歌(りっか) @rikka327

どちらの要因のみの効果(これは専門用語で主効果)があっても、交互作用がないこともある。逆に片方だけ主効果があって、交互作用があるって場合もある。交互作用は医学系の例だと「年齢が高い人だけこの薬が効く」とかそんなんに使う。確かねー

2014-07-25 00:54:45
律歌(りっか) @rikka327

まあそんな感じで要因数とか水準数(さっき行った所の作品数)によって使う検定って違うねん。他にも対応ある平均値の差の比較とかもあるで。これはさっきみたいに黒バスとよわペダという二つの作品間を比較するんじゃなくて、中学生時代のキセキと現在のキセキの身長差を比較するときに使う。

2014-07-25 00:57:39
律歌(りっか) @rikka327

つまり同じ個体での時間的な比較ね。医学系なら薬を常用して一週間後と一ヵ月後の比較とかね。作品間での比較を対応のない平均差の検定、作品内での比較を対応のある平均差の検定と呼ぶ。t検定見てて群間とか郡内とか出てきたらこれやで。

2014-07-25 01:00:17

実際にやってみた

律歌(りっか) @rikka327

みんな喜べ黒バスと弱ペダ間で有意差出たで pic.twitter.com/u3NGJA0jQ7

2014-07-25 01:30:37
拡大

画像には出てないが等分散性の検定を前もって行ってある。

黒バス群:N=15, Mean=185.73, SD=10.16
弱ペダ群:N=15, Mean=176.40, SD=9.28
t(28)=2.63,p<.05
二群の平均値に有意差が認められたため、黒バスのキャラは弱ペダのキャラより身長が高いと言う仮説が支持された。

律歌(りっか) @rikka327

とりあえず保存しとくんでエクセルデータ欲しいリア友は言ってください。ちなみに男女差が出ないように男子のみ、メンバーの抽出は黒バス三巻までしか読んだことない私が聞いたことある人を選んだので偏ってる場合がございます。これをサンプル抽出の偏りといいます。

2014-07-25 01:34:37
律歌(りっか) @rikka327

よわペダ12巻まで読んでないけど有意差出るなら葦何とかさんでなく御堂筋入れれば良かった……

2014-07-25 01:35:24

おまけ

律歌(りっか) @rikka327

知らない葦木場さんを御堂筋に変えたら有意確率0.01下回りましたわクッソワロ pic.twitter.com/icyKPefzwB

2014-07-25 01:45:06
拡大

黒バス群:N=15, Mean=185.73, SD=10.16
弱ペダ群:N=15, Mean=176.27, SD=6.57
t(28)=3.35,p<.01
二群の平均値に有意差が認められたため、黒バスのキャラは弱ペダのキャラより身長が高いと言う仮説が支持された。