교수님 안녕하세요

교수님이 쓰신 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을 사용하는 법도 언급하셨는데요.

어느 방법의 사용을 추천하시는지도 궁금합니다-!

감사합니다-!!