[R-br] Pedindo ajuda para simulação em R

[Mauro Sznelwar]

1) Este algoritmo da distribuição geométrica(p)

Conhecendo a cdf: F(x)=1-(1-p)^x+1
                F(x-1)=1-(1-p)^x
Usando o método da transformada inversa
                                 F(x-1)<R<=F(x)
1-(1-p)^x<R<=1-(1-p)^x+1
(1-p)^x+1<R-1<=1-(1-p)^x
ln(1-R)/ln(1-p)-1<=X<ln(1-R)/ln(1-p)
X=[(ln(1-R)/ln(1-p))-1]

2)Gerar a distribuição de Poisson através do processo de Poisson que pode ser simulado gerando variáveis aleatórias exponenciais até sua soma exceder 1.
Lembrar o método da inversa para distribuição exponencial
X1= -1/A(ln(Ri))
Então A1+....+Ax<=1<A1+A2+....+Ax+Ax+1
torna-se 
-somatória de i=1 até x(1/A(ln(Ri)<=1< -somatória(i=1 até x+1(1/lambda(ln(Ri))
somatória de i=1 até x((ln(Ri)>=lambda>somatória(i=1 até x+1(ln(Ri))
ln(produtória(i=1 até x(Ri)>= -lambda>ln(produtória(i=1 até x(Ri)
(produtória(i=1 até x(Ri)>=e^ -lambda>(produtória(i=1 até x(Ri)
Isto leva ao algoritmo da rejeição

Usando 
(produtória(i=1 até x(Ri)>=e^ -lambda>(produtória(i=1 até x(Ri)
Passo 1: N=0 ; P=1
Passo 2: Gere número aleatório Rn-1 e substitua P por P*Rn-1
Passo 3: Se P<e^-lambda, então aceita N=n, senão rejeita o n, soma n por um e
retorne ao passo 2. 
3) Gerar a distribuição de Pareto que é definida por f(x|alfa)=alfa*x^-alfa-1 sobre (1,infinito). Mostre que pode ser gerada como potencia -1/alfa de uma variável uniforme. 
Eu tentei fazer a inversa mas não sei se é isto: a funçao de distribuição deve ser
F(x)=1-x^-alfa e ao fazer a tranformada inversa: log(1-x)=log(u)/-alfa
4)Gerar a distribuição Beta onde U e V são i.i.d e a distribuição é (U^1/alfa)/((U^1/alfa)+V^(1/Beta)), condicional em U^1/alfa+V^1/B<=1.Tem que ser gerado com a potencia -1/alfa de uma uniforme.
Desde já agradeço




[edit]


--------------------------------------------------------------------------------


_______________________________________________
R-br mailing list
R-br em listas.c3sl.ufpr.br
https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código mínimo reproduzível.
-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20110821/7fcf5632/attachment.html>