tu pode colocar isso dentro de uma função e repetir a função 3 vezes. Substitui o comando que tu nos passou por esse, a principio ele deveria funcionar. Depois disso é só dividir as matrizes outX, outY e outZ nas 3 variaveis que tu quer, ou trabalhar nelas mesmo.
func <- function(X) {
ic_inf_X [j] = ic_inf_X [j] = ic_inf_X [j] = rep(0,n)
for (j in 1:n) {
X_e = rep(X[j], length(m))+ (runif(m, -U,U))
funX = function(X_e){mean(X_e)}
X_boot = one.boot(X_e,funX,R)
X_boot_ic=
boot.ci(X_boot,conf = c(0.95),type = c("norm"))
ic_inf_X [j]= X_boot_ic$normal[2]
ic_sup_X [j] = X_boot_ic$normal[3]
X_b[j]= X_boot$t0
}
return(cbind(ic_inf=ic_inf,ic_sup=ic_sup,b=b))
}
outX = func(X)
outY = func(Y)
outZ = func(Z)
[]'s