고수님들께 다시한번 도움을 청합니다....ㅎ
mallow cp를 구하려하는데, 페키지에 툴이 없는지 궁금합니다.
물론 구글에서 찾아봤는데,
leaps라는 함수가 있다고도 하고 Cp라고 쓰면 된다고도 합니다.
하지만 막상 해 보니까 library에 leaps가 없다고 하고 Cp도 안됩니다.
물론 x,y,beta_hat 등을 각각 불러들여 함수를 만들수도 있지만 너무 복잡하고
정해진 툴이 있는데 굳이 mallow cp를 구할때 마다 복잡하게 쓸 필요가 없을 것 같아서요....
부탁드립니다.... 이 사이트에서 정말 많은 도움을 받고 있어 너무 감사하다는 말씀 다시 올립니다.~~
첫댓글 install.packages("leaps")로 설치하셔서 library(leaps)로 불러오시면 됩니다. 다음 사이트 한번 참고하세요.
http://statweb.stanford.edu/~jtaylo/courses/stats203/R/selection/selection.R.html
개인적으로 강선찬님께 밥이라도 사야할 것 같네요. 정말 감사합니다...
(애가 아파서 응급실 갔다가 지금돌아왔어요...흑흑....T_T,,,,,그런데 이렇게 답글을 남겨 주셨네요... 고맙습니다.)
leaps를 이용해 회귀분석 모델링의 변수선택을 한 적이 있습니다. Mallow Cp, AIC, BIC 등을 옵션으로 주면 자동으로 변수 선택을 해주었던 함수였던 것으로 기억하는데..즉 Mallow Cp를 구하는 것이 목적이 아니라 변수 선택을 위한 함수인 것으로 기억합니다. 저는 이 자료 http://www.stat.columbia.edu/~martin/W2024/R10.pdf 를 참고 하였습니다.
다음 문서를 참고하셔도 도움이 되실겁니다. http://www.jstatsoft.org/v34/i12/paper