<html>
<head>
<meta content="text/html; charset=ISO-8859-1"
http-equiv="Content-Type">
</head>
<body bgcolor="#FFFFFF" text="#000000">
<div class="moz-cite-prefix">Eric<br>
<br>
Em 2008 eu precisei disto. Procurei e não achei.<br>
<br>
Como era um caso específico, calculei os quantis que precisava por
meio de bootstrap.<br>
<br>
Geração de numeros aleatórios é bem simples, bastaobservar no
código que o Walmes postou anteriormente.<br>
<pre class="moz-signature" cols="72">[]s
Leonard de Assis
assis <dot> leonard <at> gmail <dot> com
</pre>
Em 09/10/2012 14:39, Eric Ferreira escreveu:<br>
</div>
<blockquote
cite="mid:CA+R2kSjcBSCO=x_Dg+7MgYB9Lhr7EGN9CmELe6YneN-1ZFi1oQ@mail.gmail.com"
type="cite">Walmes,<br>
<br>
Eu realmente não fui muito claro. Procuro (se é que existe) um
pacote ou função que, dados duas normais, forneça os famosos 'r'
(random), 'd' (density), 'q' (quantile), 'p' (probability) da sua
razão, já levando em conta os diversos formatos possíveis da
densidade da razão como as bimodalidades, Cauchy e convergencias à
propria normal. Mas talvez não exista uma coisa muito pronta desse
jeito que preciso...<br>
<br>
Obrigado.<br>
<br>
Eric.<br>
<br>
<br>
<br>
<div class="gmail_quote">Em 9 de outubro de 2012 14:29, Walmes
Zeviani <span dir="ltr"><<a moz-do-not-send="true"
href="mailto:walmeszeviani@gmail.com" target="_blank">walmeszeviani@gmail.com</a>></span>
escreveu:<br>
<blockquote class="gmail_quote" style="margin:0 0 0
.8ex;border-left:1px #ccc solid;padding-left:1ex"><font
face="trebuchet ms,sans-serif">Eric,<br>
<br>
Não entendi bem sua mensagem. Sabemos da estatística
matemática que a razão de duas normais padrão é Cauchy. Ao
envolver outras normais não se tem um modelo único, mas
existe uma convergência para normal sob certas condições. De
qualquer forma, o que se pode fazer tranquilamente com o R é
estudar a distribuição por simulação. O que não ficou claro
é se você quer um procedimento para encontrar analiticamente
a função densidade ou por simulação.<br>
<br>
<span style="font-family:courier new,monospace"># caso de
razão de duas normais padrões<br>
n <- 1000<br>
x1 <- rnorm(n)<br>
x2 <- rnorm(n)<br>
r <- x1/x2<br>
plot(ecdf(r), xlim=c(-100,100))<br>
curve(pcauchy(x, 0, 1), add=TRUE, col=2)<br>
<br>
# caso de duas normais com média grande, r corvenge para
normal<br>
x1 <- rnorm(n, 30, 2)<br>
x2 <- rnorm(n, 30, 2)<br>
r <- x1/x2<br>
plot(ecdf(r))<br>
curve(pnorm(x, m=mean(r), sd=sd(r)), add=TRUE, col=2)<br>
qqnorm(r)<br>
den <- density(r)<br>
str(den)<br>
plot(den)<br>
fx <- with(den, approxfun(x=x, y=y))<br>
curve(fx, 0.5, 1.5)<br>
abline(v=c(0.75, 1.25))<br>
integrate(fx, 0.75, 1.25)<br>
</span><br>
À disposição.<br>
Walmes.<br>
<br clear="all">
</font><span style="font-family:trebuchet ms,sans-serif">==========================================================================</span><br
style="font-family:trebuchet ms,sans-serif">
<span style="font-family:trebuchet ms,sans-serif">Walmes
Marques Zeviani</span><br style="font-family:trebuchet
ms,sans-serif">
<span style="font-family:trebuchet ms,sans-serif">LEG
(Laboratório de Estatística e Geoinformação, 25.450418 S,
49.231759 W)</span><br style="font-family:trebuchet
ms,sans-serif">
<span style="font-family:trebuchet ms,sans-serif">Departamento
de Estatística - Universidade Federal do Paraná</span><br
style="font-family:trebuchet ms,sans-serif">
<span style="font-family:trebuchet ms,sans-serif">fone: <a
moz-do-not-send="true"
href="tel:%28%2B55%29%2041%203361%203573"
value="+554133613573" target="_blank">(+55) 41 3361 3573</a></span><br
style="font-family:trebuchet ms,sans-serif">
<span style="font-family:trebuchet ms,sans-serif">VoIP: (3361
3600) 1053 1173</span><br style="font-family:trebuchet
ms,sans-serif">
<span style="font-family:trebuchet ms,sans-serif">e-mail: <a
moz-do-not-send="true" href="mailto:walmes@ufpr.br"
target="_blank">walmes@ufpr.br</a><br>
skype: walmeszeviani<br style="font-family:trebuchet
ms,sans-serif">
</span><span style="font-family:trebuchet ms,sans-serif">twitter:
@walmeszeviani</span><br style="font-family:trebuchet
ms,sans-serif">
<span style="font-family:trebuchet ms,sans-serif">homepage: <a
moz-do-not-send="true"
href="http://www.leg.ufpr.br/%7Ewalmes" target="_blank">http://www.leg.ufpr.br/~walmes</a><span
onmouseout="cancel = false;
window.setTimeout(WRCHideContent, 1000);
clearTimeout(showTimer);" onmouseover=" var self = this;
showTimer =
window.setTimeout(function(){WRCShowContent({'rating':{'value':-1,'weight':0},'flags':{},'single':true,'ttl':7200,'expireTime':'20121009163448'},
self.className)},600);" class="wrc0"
style="padding-right:16px;width:16px;height:16px"></span></span><br
style="font-family:trebuchet ms,sans-serif">
<span style="font-family:trebuchet ms,sans-serif">linux user
number: 531218</span><br style="font-family:trebuchet
ms,sans-serif">
<span style="font-family:trebuchet ms,sans-serif">==========================================================================</span><br>
<br>
_______________________________________________<br>
R-br mailing list<br>
<a moz-do-not-send="true"
href="mailto:R-br@listas.c3sl.ufpr.br">R-br@listas.c3sl.ufpr.br</a><br>
<a moz-do-not-send="true"
href="https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br"
target="_blank">https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br</a><span
onmouseout="cancel = false;
window.setTimeout(WRCHideContent, 1000);
clearTimeout(showTimer);" onmouseover=" var self = this;
showTimer =
window.setTimeout(function(){WRCShowContent({'rating':{'value':-1,'weight':0},'flags':{},'single':true,'ttl':7200,'expireTime':'20121009163448'},
self.className)},600);" class="wrc0"
style="padding-right:16px;width:16px;height:16px"></span><br>
Leia o guia de postagem (<a moz-do-not-send="true"
href="http://www.leg.ufpr.br/r-br-guia" target="_blank">http://www.leg.ufpr.br/r-br-guia</a><span
onmouseout="cancel = false;
window.setTimeout(WRCHideContent, 1000);
clearTimeout(showTimer);" onmouseover=" var self = this;
showTimer =
window.setTimeout(function(){WRCShowContent({'rating':{'value':-1,'weight':0},'flags':{},'single':true,'ttl':7200,'expireTime':'20121009163448'},
self.className)},600);" class="wrc0"
style="padding-right:16px;width:16px;height:16px"></span>) e
forneça código mínimo reproduzível.<br>
</blockquote>
</div>
<br>
<br clear="all">
<br>
-- <br>
Dr Eric B Ferreira<br>
Exact Sciences Department<br>
Federal University of Alfenas<br>
Brazil<br>
<br>
<fieldset class="mimeAttachmentHeader"></fieldset>
<br>
<pre wrap="">_______________________________________________
R-br mailing list
<a class="moz-txt-link-abbreviated" href="mailto:R-br@listas.c3sl.ufpr.br">R-br@listas.c3sl.ufpr.br</a>
<a class="moz-txt-link-freetext" href="https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br">https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br</a>
Leia o guia de postagem (<a class="moz-txt-link-freetext" href="http://www.leg.ufpr.br/r-br-guia">http://www.leg.ufpr.br/r-br-guia</a>) e forneça código mínimo reproduzível.</pre>
</blockquote>
<br>
</body>
</html>