안녕하십니까?
선생님의 책을 통해 공부하면서 로지스틱 회귀분석을 해보고 있는데요, 궁금한 것이 있습니다.
- '웹에서 클릭만으로 하는 R통계분석' 59쪽의 결과물에서도 나온 것처럼, backward regression 최종모형에서 p-value가 0.05보다 큰데도 불구하고 모형에 남아있는 변수들이 있습니다. 이것이 함수에서 p-value를 0.05보다 큰 값(예를 들면 0.1)으로 잡아서 그런 것인지, 아니면 그냥 glm 함수에서 완벽한 모델에 이르지 못해서 그런 것인지 궁금합니다. 예전에 최종모형에서 0.05보다 큰 변수를 제거한 모델과 anova 명령어를 사용해서 비교해보니 유의한 차이가 없는 것으로 나오는 것을 본 적이 있는데, 그렇다면 더 작은 모형을 제시하는게 맞는 것인지, 아니면 최종모형을 그대로 보고해도 되는 것인지요?
- 다른 SPSS를 사용하는 통계학 교재에서 모델 검정과정에서 Hosmer Lemeshow Goodness-of-fit test 같은 것을 사용해서 모형의 적합성을 판단하는 내용을 본 적이 있는데요, R을 사용한 교재들에서는 이러한 검정에 대한 내용이 없더군요. 굳이 이러한 검정을 하지 않더라도 결과 보고에 큰 문제는 없는 건지요?
아무래도 통계분석에는 초보이고, 모르는 상황에서 혼자 힘으로 해결하려고 하다보니 통계결과보고에서 지적을 받지 않을까 하는 걱정부터 드는군요^^:
너무 초보적인 질문일지 모르겠습니다만, 부족한 부분을 깨우쳐 주시면 감사하겠습니다.
Comment 2
-
cardiomoon
2015.12.09 18:12
-
tympa
2015.12.09 22:00
아...그렇군요~~
선생님 덕분에 내용이 정리되는 것 같습니다.너무나 감사합니다^^
No. | Subject | Author | Date |
---|---|---|---|
1194 | PSM 시 범주형 변수가 총 10가지 일 때 어떻게 돌려야 하는지요 [1] | 긍정적인마음가짐 | 2023.05.06 |
1193 | 첫번째 행을 라벨로 사용 버튼 클릭시 에러 문의 [3] | 수담이 | 2023.05.04 |
1192 | 안녕하세요? web-R table 그리고 관련하여 문의드립니다. [4] | mkw7654 | 2023.05.02 |
1191 | 메타분석 U3 그래프 관련 문의드립니다 | 화이팅 | 2023.05.02 |
1190 | 교수님 안녕하세요~~ | 메시12 | 2023.04.29 |
1189 | 연속형 변수인데 범주형으로 인지됩니다. [2] | shararaku | 2023.04.28 |
1188 | 통계마당 [2] | ejo | 2023.04.27 |
1187 | moonBook 오류 관련 재문의드립니다. [1] | holyhyun | 2023.04.24 |
1186 | moonBook 사용시 오류메시지로 문의드립니다. [1] | holyhyun | 2023.04.24 |
1185 | 에러 문의 [2] | 화이팅 | 2023.04.19 |
1184 | Cox regression 에 대한 문의입니다. [2] | saga | 2023.04.18 |
1183 | PSM 문의 [3] | starwars | 2023.04.13 |
1182 | ggPS, ggPSMSummary function [6] | 로렌쵸메디치 | 2023.04.12 |
1181 | Survival 분석에 대해서 문의 드립니다. | 정훈신 | 2023.04.05 |
1180 | 안녕하세요. 분석 툴에 대해 문의드립니다. [1] | 스매직맨 | 2023.04.04 |
1179 | 드디어 publish 되었네요. [2] | 정지윤 | 2023.04.04 |
1178 | 연구동향분석 시 논문 데이터 정제 관련 [5] | vtlbtl | 2023.03.31 |
1177 | 안녕하세요 교수님 에러 문의드립니다 [2] | 소소 | 2023.03.26 |
1176 | 코딩값이 여러개인 변수에 대한 전체 p값을 알 수 있는 방법이 궁금합니다 | 쥬쥬즈주 | 2023.03.22 |
1175 | PSM 매칭 버튼 보이지 않음. | swpapa | 2023.03.19 |
1. 다중회귀모형에서 "최선의" 회귀모형을 고르는 가장 좋은 하나의 기준은 없습니다. 크게 R2를 사용하는 방법, adjusted-R2를 사용하는 방법, AIC를 사용하는 방법 등이 있는데 웹R에서 사용하는 방법은 step()함수를 사용하며 이 함수는 AIC 값을 기준으로 AIC값이 적은 모형을 선택하게 됩니다. 즉 p값을 기준으로 모형을 선택하는 것이 아니라 AIC(AKAIKE'S AN INFORMATION CRITERION)을 기준으로 모형을 선택합니다. 자세한 내용은 제가 쓴 "의학논문작성을 위한 R통계와 그래프" 책 300-308페이지 최선의 회귀모형 고르기 부분을 참조해주세요.
2. 59페이지에 있는 것은 대장암의 치료가 control, Lev,Lev+5FU 세군으로 나누어져 있습니다. rxLev로 되어 있는 항목은 control과 LEV간의 p값이며 rxLEV+5FU는 control과 LEV+5FU 군 을 비교한 p값입니다. 이러한 경우 rxLEV+5FU이 통계적으로 의미가 있으므로 모형에 계속 남아있는 것입니다.
3. Hosmer Lemeshow Goodness-of-fit test를 R에서도 할수 있지만 앞서 말씀드린대로 웹R에서는 모형의 선택은 AIC를 기준으로 한 후에 anova로 모형을 비교해줍니다. 굳이 GOF 검정을 할 필요는 없는 것 같습니다.