Skip to menu

mytable에서 하는 통계법 질문

수수 2019.07.11 13:47 Views : 1643

안녕하세요.

몇년 전 교수님 책으로 R에 입문하여 현재까지 열심히 사용 중입니다.

의료보건분야 R 관련 책은 (지금은 안 찾아보아 모르겠지만) 당시에는 이 책밖에 없었던 것 같고,

입문자에게도 쉽게 잘 읽혀서 R을 사용하게 된 확실한 계기가 되었습니다.

책을 쓰시는 것 외에도 다른 많은 활동을 하고 계신 걸 근래에 알게되었네요. 질문을 드리기에 앞서 늦었지만 매우 감사드립니다.^^


논문 작성시 및 변수들을 스크리닝 목적으로 살펴볼 때에도 moonBook 패키지가 매우 유용한 것 같은데요, mytable 함수와 관련하여 궁금한 점이 있습니다.


mytable 함수를 사용하면, 변수 특징에 따라서, 자동으로 적절한 통계법을 사용하여 p값을 구해주는데요,

어떤 변수에 구체적으로 어떤 통계법을 사용하여 나온 p값인지 궁금한데, 그에 대해서는 매뉴얼을 읽어봐도 제가 잘 못 찾은 것인지 설명을 찾지 못했습니다.


제가 이번에 하고자 하는 것은 '그룹별 연속형 변수의 차이 비교' 입니다.

제가 가지고 있는 기본 통계책을 보면 그룹별 연속형 변수 비교시에 아래와 같은 순서로 판단하여 통계방법을 선택하고 적용하는 것 같습니다.

(1) 정규성 검정[만족] -> 등분산 가정 [만족] -> Student t-test

(2) 정규성 검정[만족] -> 등분산 가정 [불만족] -> Welch t-test

(3) 정규성 검정[불만족] -> Mann Whitney U test


mytable 함수를 통해 제가 가지고 있는 데이터를 돌려서 나온 p값이, 어떤 통계방법을 통해 나온 것인지 궁금해서, 

제가 별도로 계산한 p값과 비교해보았습니다. 그랬더니,

(1) bartlett.test 등분산 가정 [만족] -> Student t-test

(2) bartlett.test 등분산 가정 [불만족] -> Welch t-test

를 적용하는 것 같았습니다. 다만, 매뉴얼에서 설명하고 있는 것처럼 정규성을 만족하지 않는 변수에 대하여 비모수적 검정방법(Mann Whitney U test)을 사용하고 있는 것 같지는 않았습니다.


mytable에서 show.all=TRUE로 설정하면 ptest라는 항목으로 변수마다 'non-normal', 'normal'을 구분해주고 있는데,

여기에서 'non-normal'이라고 표시된 변수도 Mann Whitney U test가 아니라 Student t-test에서 나온 p값과 일치하였습니다.


제가 사용해본 정규성 검정 방법과 달라서 제가 잘못 판단한 것일 수 있습니다.

그렇다면 mytable함수에서 정규성 검정은 어떤 방법으로 하는지(Lilliefors? Shapiro-Wilk? 혹은 다른 방법?)도 궁금하고,

정규성 검정 결과에 따라 Student (or Welch) t-test 혹은 Mann Whitney U test를 구분해서 하고 있는지도 궁금합니다.


평상시에는 특별히 관계가 없더라도, 특히 논문에서 결과 제시할 때에 어떤 통계법을 사용했는지에 대한 설명이 필요하기 때문에 이렇게 문의를 드립니다.

더불어, 그룹별 연속형변수 차이 분석 뿐만 아니라, 범주형 등의 경우에도 어떤 경우에 통계법을 사용하고 있는지 등도 알려주시면 더욱더 잘 활용 할 수 있을 것 같습니다^^


감사합니다.!!

No. Subject Author Date
1093 Box plot 그리프 설정에 관해 문의합니다. papahong 2022.06.20
1092 cluster data 의 risk difference 에 대해 여쭙고자 합니다. lagrima 2022.06.19
1091 Pair plot을 할 경우 [1] PipeDragoN 2022.06.13
1090 PSM에서 matching 버튼 위치 문의 통계어려워 2022.06.13
1089 PSM (4그룹) 관련 문의드립니다. emksh77 2022.06.08
1088 정회원후원하였습니다. [1] secret rosemary 2022.06.08
1087 PSM 문의 [2] secret sfumatoi 2022.06.07
1086 multipleROC 패키지에서 multivariable의 cutoff 구하기 file 찬쓰 2022.06.07
1085 survival analysis관련 문의 드립니다 [1] cebukim 2022.06.04
1084 autoReg 설치 관련 오류 [5] 닥터강 2022.05.24
1083 survival tree (RPA) 분석 관련 문의 드립니다. [2] 곰곰11 2022.05.23
1082 통계 질문입니다. [1] file 봄봄봄 2022.05.23
1081 데이터 불러오기 문제 [3] strat 2022.04.25
1080 r프로그래밍 문자는 factor로 읽을 수 없나요? [1] file 채채채추ㅐㅊ채야채 2022.04.21
1079 R 코드와 관련하여 문의 드려도 될까요..? [5] file 꾸이뭉 2022.04.19
1078 R 업데이트 후 mytable 함수 에러가 있습니다. [1] 카이토르 2022.04.18
1077 Survival curve 에서 cumulative incidence % 로 표현 가능할까요 [1] nyknicks 2022.04.16
1076 multiple logistic regression test에서 설명변수 선택에 대해 여쭙습니다. [2] file kms9057 2022.04.16
1075 PSM 문의 드립니다 [1] caesar 2022.04.07
1074 Survival regression에서 Survival probability 그래프에서 회색음영이 무엇인지요? [2] file 케로스 2022.03.24