## File Name: lsem_bootstrap_inference.R
## File Version: 0.112
lsem_bootstrap_inference <- function(parameters_boot, est, repl_factor=NULL,
bc_square=NULL)
{
R <- ncol(parameters_boot)
est_boot <- rowMeans(parameters_boot, na.rm=TRUE)
if (is.null(repl_factor)){
repl_factor <- 1/(R-1)
}
se_boot <- sqrt( rowSums( ( parameters_boot - est_boot )^2 ) * repl_factor )
bias_boot <- (est_boot - est)*repl_factor*(R-1)
est_bc <- est - bias_boot
if (!is.null(bc_square)){
pb2 <- parameters_boot[bc_square,,drop=FALSE]^2
est_boot <- rowMeans( pb2, na.rm=TRUE )
bias_boot <- (est_boot - est[bc_square]^2)*repl_factor*(R-1)
v1 <- est[bc_square]^2 - bias_boot
est_bc[bc_square] <- ifelse( v1>0, sqrt(v1), 0 )
}
#-- output
res <- list(mean_boot=est_boot, se_boot=se_boot, est_bc=est_bc,
bias_boot=bias_boot, est=est)
return(res)
}
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.