R언어에서 T-test와 Chi-square 검정을 하는 간단한 스니펫입니다.
간단한 것이지만 자주 쓰지 않으면 너무 기본이라서 오히려 기억이 잘 나지 않기 마련입니다.
t.test와 chisq.test 함수 모두 list(object)를 리턴하게 되고 멤버 변수인 p.value에 접근해서 P값을 알아낼 수 있습니다.
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 |
# T-test two vectors a <- c(1, 2, 3, 4, 5, 6, 7, 8, 9, 10) b <- c(2, 3, 4, 5, 6, 7, 8, 9, 10, 11) t.test(a, b) # 결과에 대한 자세한 설명을 화면에 출력해봐 t.test(a, b)$p.value t.test(a, b)$statistic t.test(a, b)$conf.int t.test(a, b)$estimate t.test(a, b)$null.value t.test(a, b)$alternative t.test(a, b)$method t.test(a, b)$data.name t.test(a, b)$conf.int[1] t.test(a, b)$conf.int[2] # chi-square test for two vectors a <- c(1, 2, 3, 4, 5, 6, 7, 8, 9, 10) b <- c(2, 3, 4, 5, 6, 7, 8, 9, 10, 11) chisq.test(a, b) if (chisq.test(a, b)$p.value < 0.05) { print("귀무가설 기각") } else { print("귀무가설 채택") } |