반응형
R에서 카이제곱검정 하는 방법
환자-대조군 연구 : 폐암환자 100명과 폐암이 아닌 사람 100명을 대상으로 지난 10년 동안 흡연 여부를 조사함.
|
폐암 |
비폐암 |
합계 |
흡연 |
30 |
10 |
40 |
비흡연 |
70 |
90 |
160 |
합계 |
100 |
100 |
200 |
위 예제로 카이제곱검정을 해봅시다. 흡연유무가 폐암 유무와 연관성이 있는지를 검정하는 것입니다.
아래는 R코드입니다.
r1=c(30,10)
r2=c(70,90)
data=rbind(r1,r2)
chisq.test(data,correct=FALSE)
chisq.test(data)
r2=c(70,90)
data=rbind(r1,r2)
chisq.test(data,correct=FALSE)
chisq.test(data)
아래는 결과입니다.
> chisq.test(data,correct=FALSE)
Pearson's Chi-squared test
data: data
X-squared = 12.5, df = 1, p-value = 0.000407
> chisq.test(data)
Pearson's Chi-squared test with Yates' continuity
correction
data: data
X-squared = 11.281, df = 1, p-value = 0.0007829
일반검정과 연속수정검정 모두 유의차가 있다는 것을 알 수 있습니다.
반응형
'9.비율비교 > 카이제곱검정' 카테고리의 다른 글
[R] 카이제곱검정의 사후검정 (+피셔정확검정의 사후검정) (0) | 2021.09.10 |
---|---|
R 카이제곱검정,피셔정확검정,오즈비,오즈비신뢰구간 한번에 구하기 (0) | 2020.05.07 |
댓글