Grubbs 검정으로 극단값 테스트하기

R로 시작하는 이상치 탐지 입문

Alastair Rushworth

Data Scientist

시각적 판단은 항상 신뢰할 수 없습니다!

boxplot(temperature, ylab = "Celsius")

이상값을 나타내는 상자 그림

Grubbs 검정

특정 점이 이상값인지 판단하는 통계 검정
데이터가 정규 분포를 따른다고 가정
정규성 가정을 먼저 확인해야 함

히스토그램으로 정규성 확인하기

hist(temperature, breaks = 6)

데이터 분포를 나타내는 히스토그램

대칭적이고 종 모양인가요?

Grubbs 검정 실행하기

grubbs.test() 함수 사용:

grubbs.test(temperature)

    Grubbs test for one outlier
data:  temp
G = 3.07610, U = 0.41065, p-value = 0.001796
alternative hypothesis: highest value 30 is an outlier

p-값 해석하기

grubbs.test(temperature)

    Grubbs test for one outlier

data:  temperature
G = 3.07610, U = 0.41065, p-value = 0.001796
alternative hypothesis: highest value 30 is an outlier

p-value

0에 가까울수록 - 이상값에 대한 _강한_ 근거
1에 가까울수록 - 이상값에 대한 _약한_ 근거

이상값의 행 인덱스 구하기

최댓값의 위치

which.max(weights)

최솟값의 위치

which.min(temperature)

연습해 봅시다!

R로 시작하는 이상치 탐지 입문