회귀분석을 할 때 잔차를 다시 sampling하는 방식으로 bootstraping을 해서 회귀계수를 추정하고 싶습니다.
프로그래밍이 처음이라 재미는 있는데 어떤 식으로 진행해야 하는지 감이 잘 안 잡히네요!
> Bootstrap<-function(){K<-sample(residual,size=935,replace=TRUE)
+ Y<-reg$fitted.values
+ S<-Y+K
+ RegS<-lm(S~V5)
+ return(summary(RegS))
+ }
이런 식으로 sampling을 한번 해서 보여주는 함수는 만들었는데, 이런 과정을 100번 1000번 해서 회귀계수의 평균을 추정하는 함수는 어떻게 만들어야 하는지 조언 부탁드립니다 (_ _)
본인삭제금지
글쓰기