Nothing
boot.sd1.sub <-
function(x.sub, x, weight.sub, weight, kp, nsim,
boot.index=c("r.hc", "r.is"), gamma){
if (boot.index == "r.is"){
r.is1 = r.is.sub(x.sub, x, weight.sub, weight, kp)
n <- length(x.sub)
Rbb <- NULL
for(i in 1:nsim){
s <- sample(1:n, n, replace = T)
ss <- x.sub[s]
wss <- weight.sub[s]
rownames(ss) <- NULL
R <- r.is.sub(ss, x, wss, weight, kp)
Rb <- R
Rbb <- rbind(Rb,Rbb)
}
r.se <- sd(Rbb)
ci1_up <- r.is1 + qnorm(gamma)*r.se
ci1_low <- r.is1 - qnorm(gamma)*r.se
ci2_up <- quantile(Rbb,1-(1-gamma)/2)
ci2_low <- quantile(Rbb,(1-gamma)/2)
tab <- rbind(c(ci1_low, r.is1, ci1_up), c(ci2_low, r.is1, ci2_up))
colnames(tab) <- c("ci.low", "r.is", "ci.up")
rownames(tab) <- c("norm", "quantile")
outlist <- list(se.r.is = r.se, summary = tab, boot.ind = as.vector(Rbb))
}
if (boot.index == "r.hc"){
r.hc1 = r.hc.sub(x.sub, x, weight.sub, weight, kp)$r.hc
n <- length(x.sub)
Rbb <- NULL
for(i in 1:nsim){
s <- sample(1:n, n, replace=T)
ss <- x.sub[s]
wss <- weight.sub[s]
rownames(ss) <- NULL
R <- r.hc.sub(ss, x, wss, weight, kp)$r.hc
Rb <- R
Rbb <- rbind(Rb,Rbb)
}
r.se <- sd(Rbb)
ci1_up <- r.hc1 + qnorm(gamma)*r.se
ci1_low <- r.hc1 - qnorm(gamma)*r.se
ci2_up <- quantile(Rbb,1-(1-gamma)/2)
ci2_low <- quantile(Rbb,(1-gamma)/2)
tab <- rbind(c(ci1_low, r.hc1, ci1_up), c(ci2_low, r.hc1, ci2_up))
colnames(tab) <-c("ci.low", "r.hc", "ci.up")
rownames(tab) <-c("norm", "quantile")
outlist <- list(se.r.hc = r.se, summary = tab, boot.ind = as.vector(Rbb))
}
return(outlist)
}
Any scripts or data that you put into this service are public.
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.