二群の等分散性の検定 var.test

*機能

正規母集団からの二標本の分散が等しいかどうか検定する。

*使用法

 var.test(x, y, ratio = 1, alternative = c("two.sided", "less", "greater"), conf.level = 0.95, ...)
 var.test(formula, data, subset, na.action, ...)

*引数

 x, y		2 つのデータベクトル,または,lm クラスを持つオブジェクト
 ratio		帰無仮説の下での x と y の分散比
 		デフォルトでは 1(つまり,等分散である)
 alternarive	対立仮説の種類 "two.sided", "less", "greater" のいずれか
 		最初の1文字だけでもよい
 		両側検定なら "two.sided"(デフォルト)
 		片側検定なら "less" または "greater"
 conf.level	信頼率
 		デフォルトは 0.95
 formula	「左辺値 ~ 右辺値」の形をしたモデル式
 		左辺値は母分散を検定する数値変数
 		右辺値は 2 群を表す 2 つの水準を持つ factor 変数
 data		モデル式に出てくる変数がデータフレームに含まれる場合は,それを含むデータフレーム名
 subset		観察データのサブセット
 na.action	データに NA が含まれるときに適用される関数(注)
 		デフォルトは getOption("na.action") により示されるもの(通常は na.omit 関数)
 ...		その他の引数

注:[[超訳:NAの扱い]]

*戻り値

**関数への入力

 $ alternarive: 対立仮説の種類

**関数からの出力
 $ statistic  : F 統計量
 $ parameter  : F 分布の自由度
 $ p.value    : P 値
 $ conf.int   : 母分散比の信頼区間
 $ estimate   : 標本分散比
 $ method     : 検定の種別
 $ data.name  : データの記述

*例1 2 つのデータベクトルで検定

 > x <- c(1,2,3,2,3,4,3,3,2,2,1,2,3)
 > y <- c(2,1,3,4,3,2,2,3,4,5,3,1,2)
 
 > var.test(x, y)
 
 	F test to compare two variances
 
 data:  x and y 
 F = 0.5413, num df = 12, denom df = 12, p-value = 0.3014
 alternative hypothesis: true ratio of variances is not equal to 1 
 95 percent confidence interval:
  0.1651628 1.7739390 
 sample estimates:
 ratio of variances 
          0.5412844 
 
 > var.test(lm(x ~ 1), lm(y ~ 1)) # lm クラスがどうとか言っているのは,このこと
 
 	F test to compare two variances
 
 data:  lm(x ~ 1) and lm(y ~ 1) 
 F = 0.5413, num df = 12, denom df = 12, p-value = 0.3014
 alternative hypothesis: true ratio of variances is not equal to 1 
 95 percent confidence interval:
  0.1651628 1.7739390 
 sample estimates:
 ratio of variances 
          0.5412844 

*例2 formula で検定

例1と同じデータをデータフレームにして formula を用いて検定

 > x <- c(1,2,3,2,3,4,3,3,2,2,1,2,3)
 > y <- c(2,1,3,4,3,2,2,3,4,5,3,1,2)
 > df <- data.frame(data=c(x, y), group=rep(1:2, each=13))
 
 > var.test(data ~ group, df)
 
 	F test to compare two variances
 
 data:  data by group 
 F = 0.5413, num df = 12, denom df = 12, p-value = 0.3014
 alternative hypothesis: true ratio of variances is not equal to 1 
 95 percent confidence interval:
  0.1651628 1.7739390 
 sample estimates:
 ratio of variances 
          0.5412844


トップ   編集 差分 バックアップ 添付 複製 名前変更 リロード   新規 一覧 検索 最終更新   ヘルプ   最終更新のRSS