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")
箱ひげ図の出力:
Sheeraz is a Doctorate fellow in Computer Science at Northwestern Polytechnical University, Xian, China. He has 7 years of Software Development experience in AI, Web, Database, and Desktop technologies. He writes tutorials in Java, PHP, Python, GoLang, R, etc., to help beginners learn the field of Computer Science.
LinkedIn Facebook