메뉴 건너뛰기

웹에서 하는 R 통계

mytable에서 하는 통계법 질문

2019.07.11 13:47

수수 조회 수:63

안녕하세요.

몇년 전 교수님 책으로 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를 구분해서 하고 있는지도 궁금합니다.


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

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


감사합니다.!!

번호 제목 글쓴이 날짜 조회 수
공지 묻고답하기 게시판을 운영합니다 [1] cardiomoon 2015.04.20 716
» mytable에서 하는 통계법 질문 [1] 수수 2019.07.11 63
476 여러 군의 평균 값 비교에 대해 질문드립니다. [2] jera0131 2019.07.10 42
475 mytable에서의 P value [2] file RR 2019.07.08 52
474 mytable 및 ztable 관련 질문드립니다 [1] file jera0131 2019.07.05 41
473 정회원 서버접속이 안됩니다. [1] 세희맘 2019.07.05 159
472 정회원 로그인이 안됩니다 [1] jaljary 2019.07.05 23
471 정회원 서버와 web r studio 접속이 안됩니다. [1] file jera0131 2019.07.05 37
470 정회원서버 [1] 박소이 2019.07.02 48
469 survival analysis diddbwls9 2019.07.01 52
468 ggcorrplot 에서 이분형 변수를 그냥 넣어도 될까요? Dr.Yun 2019.06.19 33
467 범주형 독립변수를 이용한 로지스틱 회귀분석 문의 드립니다. [1] 키아리 2019.06.17 149
466 Propensity score matching analysis로 질문드립니다. [1] 아누누 2019.06.07 96
465 콘솔창 실행결과가 출력되지않음. (RStudio) [1] wjdEoddl 2019.06.04 126
464 안녕하세요 서버질문드립니다 [1] 집가고싶다 2019.06.01 42
463 반복문 질문드립니다. [1] lineplay 2019.05.31 39
462 subgroup analysis 에서 forest plot을 구할때 김날 2019.05.31 55
461 paired t-test 적용유무 [2] wjdEoddl 2019.05.22 59
460 프로그램별 p-value [1] wjdEoddl 2019.05.21 61
459 rstudio-server 패키지 설치오류 file 하리보는맛있어 2019.05.14 331
458 ggPredict 그래프 theme을 바꿀 수 있을까요? [4] eyedoc광현 2019.05.13 79