# Generate Data
set.seed(10)
n <- 500
beta <- c(1,1,rep(0,8))
X <- mvtnorm::rmvnorm(n,sigma=diag(1,10))
eps <- rnorm(n)
Y <- (X%*%beta)**3+eps
# Apply bootstrap SIR
res = SIR_bootstrap(Y, X, H = 10, B = 10)
# Eigen values
plot(res,choice="eigvals")
# Estimated index versus Y
plot(res,choice="estim_ind")
Run the code above in your browser using DataLab