HOWTO · R
R での Levene テスト
このチュートリアルでは、R で Levene のテストを実行する方法を示します。
Levene の検定は、複数のグループに対して決定された変数の分散の等号を評価します。 このチュートリアルでは、R で Levene のテストを実行する方法を示します。
R での Levene テスト
Levene の検定は、複数のグループに対して決定された変数の分散の等号を評価します。 この検定では、母分散が等質性または等分散性と呼ばれる帰無仮説を調べ、k 個のサンプルの分散を比較します。ここで、k は複数のサンプルである可能性があります。
Levene の検定は、感度の低い Barlett の検定に代わるものです。 R 言語には、Levene のテストを実行するメソッド leveneTest() があります。 このメソッドは、R 言語の care パッケージからのものです。
まず、車のパッケージがまだインストールされていない場合は、インストールする必要があります。
install.packages('car')
カーパッケージがインストールされたら、Levene テストを実行できます。
コード例:
library(car)
person_weight <- data.frame(weight_program = rep(c("Program1", "Program2", "Program3"), each = 40),
weight_loss = c(runif(40, 0, 4),
runif(40, 0, 6),
runif(40, 1, 8)))
#first six rows of data frame
head(person_weight)
#conduct Levene's Test to check equality of variances
leveneTest(weight_loss ~ weight_program, data = person_weight)
上記のスニペットは、3つの異なる減量プログラムを含むデータ フレームを作成します。これは、さまざまなプログラムでどれだけの体重が減ったかを示します。 まず、コードはデータ フレームの先頭を表示し、次に Levene のテストを実行します。
出力：
weight_program weight_loss
1 Program1 1.027946
2 Program1 2.147631
3 Program1 3.181947
4 Program1 3.933521
5 Program1 3.093434
6 Program1 3.652795
Levene's Test for Homogeneity of Variance (center = median)
Df F value Pr(>F)
group 2 5.2128 0.006782 **
117
---
Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
Levene の検定の p 値は 0.006782 で、有意水準 0.05 よりも小さいです。 Levene の検定概念に基づいて帰無仮説を棄却するため、3つのグループ間の分散は等しくありません。
ボックス グラフをプロットして、重みの分布を示しましょう。
boxplot(weight_loss ~ weight_program,
data = person_weight,
main = "Weight Loss Distribution by the Weight Program",
xlab = "Weight Program",
ylab = "Weight Loss",
col = "grey90",
border = "black")
箱ひげ図の出力:
