[R] Levene's test (등분산검정)

Levene's test 는 집단들의 분산이 같은지 여부를 검정할 때 사용합니다. 예제를 통해 이해해봅시다. 

 

패키지 설치 및 불러우기

먼저 패키지를 설치합니다. 

install.packages("lawstat")

패키지를 불러옵니다.

 

library(lawstat)

 

데이터 생성

예제에 사용될 테이터를 만들어봅시다. 두 집단의 키 데이터입니다. 그룹별로 50명을 추출하였습니다. 

set.seed(1)
A=rnorm(50,175,10)
B=rnorm(50,185,10)

stack 형태의 데이터로 만들어줍니다. stack 형태의 데이터는 각 열에 독립변수와 종속변수가 입력된 데이터입니다. 위 데이터의 경우 독립변수는 그룹의 종류이고 종속변수는 키입니다. 

x=c(rep(1,50),rep(2,50))
y=c(A,B)
df=data.frame(X=x,Y=y)

 

생성한 데이터는 아래와 같습니다. 

 

 

Levene's test 수행

levene.test 함수를 이용합니다. 결과는 아래와 같습니다. 

> levene.test(df$Y,df$X,location="mean")

Classical Levene's test based on the absolute deviations from the mean ( none not applied because the location is not set to median )

data:  df$Y
Test Statistic = 1.1157, p-value = 0.2934

 

 

#전체 코드

#패키지 불러오기
library(lawstat)

#데이터 생성
set.seed(1)
A=rnorm(50,175,10)
B=rnorm(50,185,10)

#데이터를 stack 형태로 변형
x=c(rep("A",50),rep("B",50))
y=c(A,B)
df=data.frame(X=x,Y=y)

#levene test 수행
levene.test(df$Y,df$X,location="mean")