data(affluence)
affluence$weight <- rep(1, nrow(affluence))
aff.sub <- subset(affluence, education == 2)
boot.sd1.sub(aff.sub$income, affluence$income, aff.sub$weight, affluence$weight,
0.9, 10, "r.is", 0.95)
boot.sd1.sub(aff.sub$income, affluence$income, aff.sub$weight, affluence$weight,
0.9, 10, "r.hc", 0.95)
Run the code above in your browser using DataLab