[R 통계분석] Brown-Forsythe test (등분산검정)

[R 통계분석] Brown-Forsythe test (등분산검정)

#1. 방법

levene.test(y, group, location = c("median", "mean", "trim.mean"),
  trim.alpha = 0.25, bootstrap = FALSE, num.bootstrap = 1000,
  kruskal.test = FALSE, correction.method = c("none",
  "correction.factor", "zero.removal", "zero.correction"))

y자리에 전체 데이터, group 자리에 분류기준을 입력한다. 예를들어 남녀 그룹의 키를 비교하는 경우, y자리에는 전체 키 데이터를 입력하고 group 자리에는 각 데이터의 성별을 입력한다. levene's test에서 location 을 median으로 선택한 경우가 Brown-Forsythe test 입니다.

#2. 예제

> 코드

#Levene's test를 하기 위해 lawstat 패캐지를 먼저 설치합니다.
install.packages("lawstat")

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

#표본생성
set.seed(1)
group1=rnorm(50,0,1)
group2=rnorm(50,0,1)

#형식에 맞게 수정
y=c(group1,group2)
group=c(rep(1,50),rep(2,50))

#Brown-Forsythe test 수행. levene 검정에서 median옵션을 적용하면 됩니다.
levene.test(y,group,location="median")

> 실행결과

> levene.test(y,group,location="median")

    Modified robust Brown-Forsythe Levene-type test based on the absolute
    deviations from the median

data:  y
Test Statistic = 1.1141, p-value = 0.2938