[R-br] Dificuldade na utilização de um loop duplo para a criação de várias amostras

Thiago Touya thiagotouya em gmail.com
Quarta Novembro 19 16:23:15 BRST 2014


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)
-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20141119/6193dd1e/attachment.html>


Mais detalhes sobre a lista de discussão R-br