<!DOCTYPE HTML PUBLIC "-//W3C//DTD HTML 4.0 Transitional//EN">
<HTML><HEAD>
<META http-equiv=Content-Type content="text/html; charset=iso-8859-1">
<META content="MSHTML 6.00.2900.3429" name=GENERATOR>
<STYLE></STYLE>
</HEAD>
<BODY bgColor=#ffffff>
<DIV>1) Este algoritmo da distribuição geométrica(p)</DIV>
<DIV><FONT face=Arial size=2></FONT> </DIV>
<DIV><FONT face=Arial size=2>Conhecendo a cdf:
F(x)=1-(1-p)^x+1<BR>
F(x-1)=1-(1-p)^x<BR>Usando o método da transformada
inversa<BR>
F(x-1)<R<=F(x)<BR>1-(1-p)^x<R<=1-(1-p)^x+1<BR>(1-p)^x+1<R-1<=1-(1-p)^x<BR>ln(1-R)/ln(1-p)-1<=X<ln(1-R)/ln(1-p)<BR>X=[(ln(1-R)/ln(1-p))-1]</FONT></DIV>
<DIV> </DIV>
<DIV>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.</DIV>
<DIV><FONT face=Arial size=2>Lembrar o método da inversa para distribuição
exponencial<BR>X1= -1/A(ln(Ri))<BR>Então
A1+....+Ax<=1<A1+A2+....+Ax+Ax+1<BR>torna-se <BR>-somatória de i=1 até
x(1/A(ln(Ri)<=1< -somatória(i=1 até x+1(1/lambda(ln(Ri))<BR>somatória de
i=1 até x((ln(Ri)>=lambda>somatória(i=1 até
x+1(ln(Ri))<BR>ln(produtória(i=1 até x(Ri)>= -lambda>ln(produtória(i=1 até
x(Ri)<BR>(produtória(i=1 até x(Ri)>=e^ -lambda>(produtória(i=1 até
x(Ri)<BR>Isto leva ao algoritmo da rejeição</FONT></DIV>
<DIV><FONT face=Arial size=2></FONT> </DIV>
<DIV><FONT face=Arial size=2>Usando <BR>(produtória(i=1 até x(Ri)>=e^
-lambda>(produtória(i=1 até x(Ri)<BR>Passo 1: N=0 ; P=1<BR>Passo 2: Gere
número aleatório Rn-1 e substitua P por P*Rn-1<BR>Passo 3: Se P<e^-lambda,
então aceita N=n, senão rejeita o n, soma n por um e<BR>retorne ao passo 2.
</FONT></DIV>
<DIV><FONT face=Arial size=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. </FONT></DIV>
<DIV><FONT face=Arial size=2>Eu tentei fazer a inversa mas não sei se é isto: a
funçao de distribuição deve ser</FONT></DIV>
<DIV><FONT face=Arial size=2>F(x)=1-x^-alfa e ao fazer a tranformada inversa:
log(1-x)=log(u)/-alfa</FONT></DIV>
<DIV><FONT face=Arial size=2>4)<FONT face="Times New Roman" size=3>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.</FONT></FONT></DIV>
<DIV><FONT face=Arial size=2><FONT face="Times New Roman" size=3>Desde já
agradeço</FONT></DIV>
<DIV> </DIV></FONT>
<DIV><FONT face=Arial size=2><SPAN class=Apple-style-span
style="WORD-SPACING: 0px; FONT: 13px/19px sans-serif; TEXT-TRANSFORM: none; COLOR: rgb(0,0,0); TEXT-INDENT: 0px; WHITE-SPACE: normal; LETTER-SPACING: normal; BACKGROUND-COLOR: rgb(255,255,255); orphans: 2; widows: 2; webkit-text-decorations-in-effect: none; webkit-text-size-adjust: auto; webkit-text-stroke-width: 0px">
<DL style="MARGIN-TOP: 0.2em; MARGIN-BOTTOM: 0.5em">
<DD style="MARGIN-BOTTOM: 0.1em; MARGIN-LEFT: 1.6em; LINE-HEIGHT: 1.5em"><FONT
size=2></FONT> </DD></DL>
<P style="MARGIN: 0.4em 0px 0.5em; LINE-HEIGHT: 1.5em"><FONT
size=2></FONT> </P>
<H2
style="OVERFLOW-Y: hidden; FONT-WEIGHT: normal; FONT-SIZE: 19px; BACKGROUND-IMAGE: none; OVERFLOW-X: hidden; PADDING-BOTTOM: 0.17em; MARGIN: 0px 0px 0.6em; WIDTH: auto; COLOR: black; PADDING-TOP: 0.5em; BORDER-BOTTOM: rgb(170,170,170) 1px solid; background-origin: initial; background-clip: initial"><SPAN
class=editsection style="FONT-SIZE: 13px; FLOAT: right; MARGIN-LEFT: 5px">[<A
title="Edit section: Generating bounded Pareto random variables"
style="BACKGROUND-IMAGE: none; COLOR: rgb(6,69,173); TEXT-DECORATION: none; background-origin: initial; background-clip: initial"
href="http://en.wikipedia.org/w/index.php?title=Pareto_distribution&action=edit§ion=21">edit</A>]</SPAN></H2></SPAN></FONT></DIV>
<DIV> </DIV>
<P>
<HR>
<P></P>_______________________________________________<BR>R-br mailing
list<BR>R-br@listas.c3sl.ufpr.br<BR>https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br<BR>Leia
o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código mínimo
reproduzível.</BODY></HTML>