★ 統計学的に**座は死亡事故を起こしやすいか? ★

1004. 統計学的に**座は死亡事故を起こしやすいか? うお座 2003/11/12 (水) 13:27
├1008. 自己レス 有意とはいえないでしょう 結論 うお座 2003/11/12 (水) 14:17
├1007. Re: 統計学的に**座は死亡事故を起こしやすいか? 青木繁伸 2003/11/12 (水) 14:03
└1006. Re: 統計学的に**座は死亡事故を起こしやすいか? 青木繁伸 2003/11/12 (水) 13:54
 └1009. Re^2: 統計学的に**座は死亡事故を起こしやすいか? うお座 2003/11/12 (水) 14:19


1004. 統計学的に**座は死亡事故を起こしやすいか? うお座  2003/11/12 (水) 13:27
【社会】死亡事故最多は「うお座」−北海道警が星座別に分析
http://www.sanspo.com/sokuho/1111sokuho042.html
のデータを使ってカイ二乗検定をしてみました。
(星座別走行距離は等しいと仮定$)

『1433件の交通死亡事故で,「うお座」(149件),「やぎ座」(138件),「かに座」(131件)の順。最も少ない「さそり座」は99件』との情報しかないので,残りの8星座での発生数が等しいと仮定するとP=0.27
残り8星座が非常に偏っていると仮定するとP=0.003
となりました。

ポアソン事象だと思いますが,総計が決まっており,トップ3と最下位の発生数が分かっているときに他の発生数をどのように仮定するのがいいのでしょうか? 

     [このページのトップへ]


1008. 自己レス 有意とはいえないでしょう 結論 うお座  2003/11/12 (水) 14:17
> ポアソン事象だと思いますが,総計が決まっており,トップ3と最下位の発生数が分かっているときに他の発生数をどのように
>仮定するのがいいのでしょうか? 

少し考えて,lamda=1433/12でトップ3とボトム1つを除いた間の8つと仮定するのがリーズナブルだとおもいました。
Rだと
a<-rpois(12,1433/12)
で
for (I in 1:100) {
  a<-rpois(12,1433/12)
  x<-c(143,138,131,sort(a)[4:11],99)
  print(chisq.test(x)$p.value)
}
で100回ほどポアソンで残り8星座の発生数を仮定して
カイ二乗検定すると全てp>.05となりました。
今回の星座によるばらつきは統計学的には有意とは
いえないであろうと思えるようになりました。

     [このページのトップへ]


1007. Re: 統計学的に**座は死亡事故を起こしやすいか? 青木繁伸  2003/11/12 (水) 14:03
> 【社会】死亡事故最多は「うお座」−北海道警が星座別に分析

作家の星座についても,過去ログにあります。
http://aoki2.si.gunma-u.ac.jp/lecture/mb-arc/arc020/364.html

「星座と○○」ということで,○○のところをいろいろ変えてα=0.05でやってみれば,そのような試み20件あたり1件の大発見があるのでしょう。例数が多ければもっと「あたり」の確率は高くなるでしょう。---週刊誌などが喜びそう。
自民党議員は○○座が多い。民主党議員は○○座が多い。
首相経験者は○○座が多い。高級官僚には○○座が多い。
有名タレントには○○座が多い。医者には○○座が多い。
金持ちは○○座が多い。。。。。

     [このページのトップへ]


1006. Re: 統計学的に**座は死亡事故を起こしやすいか? 青木繁伸  2003/11/12 (水) 13:54
> 総計が決まっており,トップ3と最下位の発生数が分かっているときに他の発生数をどのように仮定するのがいいのでしょうか? 

似たようなもの(?)が,過去ログにあります(私は意外な感じがしました)

ジップの法則
http://aoki2.si.gunma-u.ac.jp/lecture/mb-arc/arc025/086.html

     [このページのトップへ]


1009. Re^2: 統計学的に**座は死亡事故を起こしやすいか? うお座  2003/11/12 (水) 14:19
青木先生いつも有難うございます。

> ジップの法則
> http://aoki2.si.gunma-u.ac.jp/lecture/mb-arc/arc025/086.html

はこの場合には当てはまらないのではないでしょうか?

     [このページのトップへ]


● 「統計学関連なんでもあり」の過去ログ--- 026 の目次へジャンプ
● 「統計学関連なんでもあり」の目次へジャンプ
● 直前のページへ戻る