<div dir="ltr"><div>Boa tarde!</div><div><br></div><div>Senhores, peço sua ajuda para criar um código que gere um número 'n' de amostras de tamanho 'n'.</div><div><br></div><div>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'.</div><div><br></div><div>Segue o código que estou usando para a geração de uma amostra:</div><div><br></div><div>## Tamanho da Amostra<br>n=100</div><div><br></div><div>## Parâmetros 1<br>a <- c(4.5,1.5)<br>b <- c(3.0,9.0)<br>p <- c(0.5,0.5)</div><div><br></div><div>## Função Acumulada Inversa<br>fmenos1_kuma = function(x,a,b){<br> (1-(1-x)^(1/b))^(1/a)<br>}</div><div><br></div><div>## Mistura <br>mkuma = function(x){<br> (a[1]*b[1]*(x^(a[1]-1))*(1-x^a[1])^(b[1]-1))*p[1] + <br> (a[2]*b[2]*(x^(a[2]-1))*(1-x^a[2])^(b[2]-1))*p[2]<br>}</div><div><br></div><div>## Simulação da amostra <br></div><div> u1 <- runif(n)<br> u2 <- runif(n)<br> x  <- numeric(n)<br> for (i in 1:n){<br>  if(u1[i] <  p[1]) x[i]= fmenos1_kuma(u2[i],a[1],b[1])<br>  if(u1[i] >= p[1]) x[i]= fmenos1_kuma(u2[i],a[2],b[2])<br> }<br><br clear="all"></div><div><br></div><div>O que eu pensei em fazer foi adicionar um outro 'For' no passo da simulação da </div><div>amostra, mas  não consegui alimentar um objeto 'y' com os valores gerados para cada amostra.</div><div><br></div><div><br></div><div>Obrigado.</div><div><br>-- <br></div><div class="gmail_signature">Att,<br>Thiago Morais de Carvalho (61 - 8169-5094)<br></div>
</div>