Comecei gerando uma amostra de 500 números aleatórios uniformemente distribuídos entre 0 e 1 usando o código abaixo:

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

Agora, preciso escrever um psuedocode que gere 10000 amostras de N = 500 para uma simulação de MC, calcule a média do meu X recém-criado e armazene o número da iteração e o valor médio em um objeto de resultado. Eu nunca tentei isso, e até agora eu tenho isso:

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

Finalmente, uma vez que isso seja feito, eu devo executá-lo, depois obter mediana, média e min / max das médias de amostra acumuladas e salvá-las em um quadro de dados chamado MC.table. (Observe também, acima, eu não tenho idéia do porque há um "4" no código da matriz --- estou trabalhando fora dos exemplos anteriores). Qualquer conselho ou ajuda seria muito apreciado.

EDIT: Eu tenho um exemplo que pode funcionar, mas eu realmente não entendo o que está acontecendo com isso, por favor, elaborar sobre a sua validade para isso:

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) } }

Se você pudesse dedicar um tempo para explicar o que isso significa, isso ajudaria muito! Obrigado!