Boa tarde! Senhores, peço sua ajuda para criar um código que gere um número 'n' de amostras de tamanho 'n'. Atualmente, consigo fazer apenas uma amostra de tamanho 'n', mas gostaria de fazer de um modo prático a geração de outras 'n-1' amostras. No final gostaria de armazenar tudo numa matriz 'nxn'. Segue o código que estou usando para a geração de uma amostra: ## Tamanho da Amostra n=100 ## Parâmetros 1 a <- c(4.5,1.5) b <- c(3.0,9.0) p <- c(0.5,0.5) ## Função Acumulada Inversa fmenos1_kuma = function(x,a,b){ (1-(1-x)^(1/b))^(1/a) } ## Mistura mkuma = function(x){ (a[1]*b[1]*(x^(a[1]-1))*(1-x^a[1])^(b[1]-1))*p[1] + (a[2]*b[2]*(x^(a[2]-1))*(1-x^a[2])^(b[2]-1))*p[2] } ## Simulação da amostra u1 <- runif(n) u2 <- runif(n) x <- numeric(n) for (i in 1:n){ if(u1[i] < p[1]) x[i]= fmenos1_kuma(u2[i],a[1],b[1]) if(u1[i] >= p[1]) x[i]= fmenos1_kuma(u2[i],a[2],b[2]) } O que eu pensei em fazer foi adicionar um outro 'For' no passo da simulação da amostra, mas não consegui alimentar um objeto 'y' com os valores gerados para cada amostra. Obrigado. -- Att, Thiago Morais de Carvalho (61 - 8169-5094)