Comencé generando una muestra de 500 números aleatorios distribuidos uniformemente entre 0 y 1 usando el siguiente código:

set.seed(1234)
X<-runif(500, min=0, max=1)

Ahora, necesito escribir un psuedocode que genere 10000 muestras de N = 500 para una simulación de MC, calcular la media de mi X recién creada y almacenar el número de iteración y el valor de la media en un objeto de resultado. Nunca he intentado esto, y hasta ahora tengo esto:

n.iter <-(10000*500)
results <- matrix (0, n.iter, 4)

Finalmente, una vez que esto se logre, lo ejecutaré, luego obtendré la mediana, la media y el mínimo / máximo de las medias de muestra acumuladas y las guardaré en un marco de datos llamado MC.table. (También tenga en cuenta, arriba, no tengo idea de por qué hay un "4" en el código de matriz, estoy trabajando en los ejemplos anteriores). Cualquier consejo o ayuda sería muy apreciado.

EDIT: Tengo un ejemplo que puede funcionar, pero realmente no entiendo lo que está pasando con él, así que por favor explique su validez para esto:

Ni <- 10000
n <- 500
c <- 0

for (i in n){
for (j in 1:Ni){
c <- c+ 1
d <- data.frame (x= , y= )
results [c,1] <- c
results [c,2] <- j
results [c,3] <- i
results [c,4] <- something( d$x, d$y)
rm (d) } }

Si pudieras tomarte el tiempo de explicar lo que eso significa, ¡eso ayudaría mucho a ayudarme! ¡Gracias!