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 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 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 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 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 href="mailto:R-br@listas.c3sl.ufpr.br">R-br@listas.c3sl.ufpr.br</a><br>
<a 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 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>