교수님 안녕하세요
교수님이 쓰신 r통계와 그래프 및 웹에서 하는 r 통계를 참고로 많이 배우고 있습니다.
지금 초보적이지만 r을 쓰려고 노력하고 있는데요
여기 web-r에서 하는 것처럼 multivariable logistic regression에서 p value가 0.2이상인 것을 자동으로 골라서
multivariable analysis를 r 에서도 구현하고 싶습니다.
현재는 제가 하고 있는 방법은 각각 univariate로 (glm) 돌려서 그중에서 매뉴얼로 제가 p가 0.2이상인것을 골라서
full model을 glm로 만든후 step을 사용해서 시행하고 있습니다
(spss에서 돌리는 것과 비슷하게 하고 있습니다)
full.model=glm(Murder-a+b+c+d ,data=states)
reduced .model=step (full.model ,direction= "backward" )
하지만 변수가 많아지니까 30개 넘는 것에서 고르는 곳도 생각보다 시간도 오래걸리고 실수할 확률도 높아지는거 같아서요
혹시 full.model에서 전부 넣고 step을 사용하면 자동으로 p value 0.2이상인 것만 multivariable에 들어가게 되는건가요??
web-r처럼 p value 를 조정하려면 어떻게 해야하나요?;
그리고 자동이라면 multivable로 들어가는 것은 무엇인지 어떻게 알게되나요??
그리고 r통계와 그래프에서 보면 AIC로 하는 방법도 있다고 하시고 또한 shrinkageO 패키지와 all subset regression을 사용하는 법도 언급하셨는데요.
어느 방법의 사용을 추천하시는지도 궁금합니다-!
감사합니다-!!
Comment 2
-
cardiomoon
2017.07.14 15:18
-
권혜미
2017.07.14 17:26
그렇다면 혹시......
univariate에서 p<0.2이하인 것을 자동으로 선택하는 패키지나 코드가 잇을까요?ㅠㅠ
현재는 variable들 glm(A~variable1개, data=abc)을 lapply로 해서 나온 결과를 결과 나온걸 제가 매뉴얼로 고르고 잇습니다.ㅠ
No. | Subject | Author | Date |
---|---|---|---|
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 |
1174 | group-based trajectory modeling | swpapa | 2023.03.19 |
1) step함수는 p값과 상관없이 AIC 값을 기준으로 모형을 골라줍니다.
2) full model에 어떤 설명변수를 넣을지는 연구자가 정하는 것입니다. 웹R에서는 univariate analysis를 해서 사용자가 정해준 p 값의 threshold 이하인 것을 골라서 full model에 넣어줍니다.
(하나하나 비교하는 것입니다.)
3) step() 함수의 결과와 all subset regression의 결과, Cp plot의 결과가 일치하면 좋겠지만 일치하지 않는다면 연구자가 골라야 하겠죠.^^