반응형 9.비율비교/카이제곱검정3 [R] 카이제곱검정의 사후검정 (+피셔정확검정의 사후검정) 먼저 rcompanion 패키지를 설치합니다. 해당 라이브러리 안에 있는 데이터와 함수를 이용할 것이기 때문입니다. install.packages("rcompanion") 설치한 패키지를 불러옵니다. library(rcompanion) 사용할 데이터는 아래와 같습니다. > data(Anderson) > Anderson Pass Fail Bloom 21 5 Cobblestone 6 11 Dougal 7 8 Heimlich 27 5 사후검정은 pairwiseNominalIndependence 함수를 사용합니다. > pairwiseNominalIndependence(Anderson,chisq = TRUE,method="bonferroni") Comparison p.Fisher p.adj.Fisher p.Gte.. 2021. 9. 10. R 카이제곱검정,피셔정확검정,오즈비,오즈비신뢰구간 한번에 구하기 R 카이제곱검정,피셔정확검정,오즈비,오즈비신뢰구간 한번에 구하기 환자-대조군 연구 : 폐암환자 100명과 폐암이 아닌 사람 100명을 대상으로 지난 10년 동안 흡연 여부를 조사함. 폐암 비폐암 합계 흡연 30 10 40 비흡연 70 90 160 합계 100 100 200 위 예제를 이용할 것입니다. 굉장히 유용한 패키지가 있어서 사용하겠습니다. 아래 결과들을 한번에 출력해줍니다. 카이제곱검정 피셔정확검정 오즈비 오즈비신뢰구간 먼저 패키지를 설치합시다. install.packages("epitools") 아래와 같이 코딩합니다. oddsratio.wald 함수를 사용합니다. riskratio.wald를 사용하면 오즈비 대신 위험비를 계산해줍니다. library(epitools) r1=c(30,10) r2.. 2020. 5. 7. R에서 카이제곱검정 하는 방법 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) 아래는 결과입니다. > chisq.test(data,correct=FALSE) Pearson's Chi-squared test data: data X-squared = 12.5, df = 1, .. 2020. 5. 5. 이전 1 다음 반응형