R block resampling by unique identifier for bootstrap

Question

I am attempting to block bootstrap a dataset using R. I have a data frame of firms in counties. I want to sample counties with replacement, then build a dataset with all firms in that sample of counties (with replacement). I run a regression on the new dataset. Then I sample again.

I have a for loop that works like so:

for(j in 1:10000){
y=NULL
for(i in 1:length(unique(data$firm_id))){
    y=rbind(y, data[which(data$county_id==sample(unique(data$county_id), replace=T)[i]),])
}
    a=rbind(a, lm(profit~employees, data=y)$coefficients)
}

Unfortunately, this sort of for loop in R is extremely slow and computationally expensive. Is it possible to implement this using a more efficient apply function?

There is a `library(boot)` package that may answer your question. However, you can probably do the trick with `sample` and simplifying your code. Without a small example, not sure how to help out. — Drew75, Feb 02 '14 at 08:55
Also look at http://www.burns-stat.com/pages/Tutor/R_inferno.pdf with a few tips on how to optimize your code. — Roman Luštrik, Feb 02 '14 at 09:08
http://stackoverflow.com/questions/11919808/block-bootstrap-from-subject-list might help — user20650, Feb 02 '14 at 12:48
Thanks! Not sure how I missed that question. That plus the code below helped me figure it out. — DannyMatt, Feb 06 '14 at 05:29

score 2 · Accepted Answer · answered Feb 02 '14 at 09:43

2

something like this could help:

positions<-replicate(1000, sample(1:nrow(df), nrow(df), T))

apply(positions, 2, function(i) lm(yvar[i]~xvar[i], df)$coef)

answered Feb 02 '14 at 09:43

Davide Passaretti

2,741
1
21
32

R block resampling by unique identifier for bootstrap

1 Answers1