스타의 도서관/통계학

R을 이용한 통계 분석 - 모평균의 t검정

스타(star) 2012. 4. 24. 00:54

작년도 대기업의 대졸 신입 사원 월임금은 평균 160만원 이었다. 금년도 대졸자의 초임을 알아보기 위하여 15명을 단순확률 추출하여 성별과 월임금을 조사하니 다음과 같다.

ex salary
1 163
1 162
2 157
1 172
1 161
2 165
1 162
2 165
1 158
1 159
2 160
2 161
2 165
1 168
1 170
이 자료를 통해 대졸 신입 초임이 작년보다 인상되었는지 유의수준 5%로 가설 검정하라.

우선 이 가설은

H0 : u = 160만원 - 귀무가설

H1 : u > 160만원 - 대립가설이다.

 

위와 같은 텍스트 파일을 만들고, 테이블을 로드해 오겠습니다.

ex <- read.table("ex69.txt", header=T)

header=T 파라메터는 맨 위의 헤더가 변수명이라는 것을 말합니다.

 

다음으로 모평균 t검정을 해봅니다.

 t.test(ex$salary, alternative=c("greater"), mu=160)

ex$salary는 ex데이터의 salary 변수를 의미합니다. 또한 alternative는 대립 가설을 말합니다. "two.sided", "less", "greater"중 하나를 선택할 수 있습니다. mu는 검정할 평균값을 말합니다.

위의 문장을 해석하면, t검정을 해서 대립가설의 평균이 160보다 큰지 검정하는 것입니다.

출력 결과에서 p-value가 유의수준은 0.05 뿐만 아니라 0.01 보다 작은 0.006266이라는 값이 나왔습니다.

p-value는 측정된 표본 평균이 모든 가능한 표본평균 중에서 어디에 위치하는지 알려줍니다. 때문에 매우 높은 확률로 160 이상이라고 할 수 있겠네요.