Okumura's Blog

Twitterでの某先生のつぶやきがきっかけで，ググって

• 等分散の検定
• マイクロソフト ExcelのFtestのバグ
• ExcelのF検定
• ひどい話です！のF検定のところ
• エクセル2003でのＦ検定時のＰ値 | OKWave

などをざっと見ながら，等分散の検定にある練習用データをExcel 2010で試してみた（図）。結果はExcel 2003とまったく同じ。古いExcelでは片側検定の結果を間違って両側検定と表示していたが，Excel 2003以降は正しく片側検定と表示するようになった。ただ，この場面では両側検定が普通なので，直す場所が違ったようにも思う。「P(F<=f) 片側」という書き方も微妙である。「F 境界値 片側」もよく考えないと誤解しやすい。

MacのExcel 2008でやってみようとしたら，そうだ，Mac版には分析ツールがなかったのだ：

• Mactopia Japan | Office 2008 Q & A ページの「ソルバーがない！？！　分析ツールもない！？！ VB も使えない！？！」

ついでにRでやってみた結果：

> ftall = c(170, 166, 157, 164, 157, 160, 163, 160, 153, 165, 162,
  150, 165, 151, 165, 169, 158, 153, 172, 162, 160, 175, 165, 155,
  177, 161, 160, 166, 162, 161, 165, 158, 158, 164, 160, 159, 165,
  163, 160, 175, 159, 161, 173, 158, 159, 162)
> mtall = c(175, 178, 185, 163, 165, 175, 168, 180, 180, 170, 167,
  175, 164, 178, 175, 169, 171, 190, 173, 171, 175, 171, 166, 181,
  175, 163, 180, 188, 180, 185, 181, 167, 178, 177, 172, 198, 171,
  169, 190, 164, 170, 188, 177, 169, 186, 179, 176, 190, 171, 171,
  186, 169, 179, 180)
> var.test(ftall, mtall)

	F test to compare two variances

data:  ftall and mtall 
F = 0.5655, num df = 45, denom df = 53, p-value = 0.05203
alternative hypothesis: true ratio of variances is not equal to 1 
95 percent confidence interval:
 0.3223957 1.0050989 
sample estimates:
ratio of variances 
         0.5654643