<div dir="ltr">Prezados,<div><br></div><div>Agradeço imensamente a todos que colaboraram. As propostas de Benilton e Ivan funcionaram com muita velocidade, mesmo quando solicitei uma simulação de 1534 dados. No entanto, quando apliquei o teste de aderência (Shapiro-Wilk), o valor de p indicou que os dados NÃO obedeciam a uma distribuição normal.</div>
<div><br></div><div>Testei tanto as duas propostas quanto o meu código usando um n=100, pensando se seria um problema do tamanho da amostra. Com o meu "loop", o teste indicou distribuição Normal (p>0.05), mas com as duas propostas deu um p<0.05.</div>
<div><br></div><div>Creio que meus conhecimentos ainda iniciais no software não me permitiram compreender a proposta de Walmes para eu testar. Não consegui identificar onde eu deveria incluir meus parâmetros de n, valor máx, valor mín, média e desvio padrão.</div>
<div><br></div><div>De qualquer forma, agradeço muito a disposição de todos em ajudar.</div><div><br></div><div>Abraços fraternos,</div></div><div class="gmail_extra"><br><br><div class="gmail_quote">Em 29 de julho de 2014 14:19, walmes . <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"><div dir="ltr"><div class="gmail_default" style="font-family:trebuchet ms,sans-serif">Se eu fosse gerar da normal truncada, eu iria proceder usando o método da inversa da acumulada. Talvez seja mais rápido. O com while pode demorar demais se a taxa de rejeição for grande. Por exemplo, o tempo para gerar sem números entre 0 e 1 é maior do que para gerar entre -1 e 1 porque a taxa de rejeição no primeiro caso é maior.<br>

<br><span style="font-family:courier new,monospace">par(mfrow=c(1,2))<br>curve(dnorm(x, 0, 1), -3, 3)<br>abline(v=c(-0.5, 2), lty=2)<br>curve(pnorm(x, 0, 1), -3, 3)<br>abline(v=c(-0.5, 2), lty=2)<br>layout(1)<br><br>lim <- pnorm(c(-0.5, 2), 0, 1)<br>

u <- runif(100, lim[1], lim[2])<br>y <- qnorm(u, 0, 1)<br><br>par(mfrow=c(1,2))<br>curve(dnorm(x, 0, 1), -3, 3)<br>abline(v=c(-0.5, 2), lty=2)<br>rug(y)<br>curve(pnorm(x, 0, 1), -3, 3)<br>abline(v=c(-0.5, 2), lty=2)<br>

rug(u, side=2) <br>rug(y)<br>layout(1)<br></span><br>À disposição.<span class="HOEnZb"><font color="#888888"><br>Walmes.<br><br></font></span></div></div>
<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><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>) e forneça código mínimo reproduzível.<br></blockquote></div><br><br clear="all"><div><br></div>-- <br>
<div dir="ltr">MARCOS BISSOLI<br>Faculdade de Nutrição<br>Universidade Federal de Alfenas<br><br>E-mail: <a href="mailto:mbissoli@gmail.com" target="_blank">mbissoli@gmail.com</a><br>Twitter: #mbissoli<div>Facebook: <a href="https://www.facebook.com/MarcosBissoli" target="_blank">https://www.facebook.com/MarcosBissoli</a><br>
<br>Alfenas, Minas Gerais, Brasil<br><br><br>*****Pense na Natureza antes de Imprimir*****<br>Divulgue ON-LINE<br><br>Eu apoio a ENEN "na luta por um Brasil sem fome"<br><br>"por ĉiu popolo ties propran lingvon, por ĉiuj popoloj la esperantan"<br>
(para cada povo sua própria língua, para todos os povos o Esperanto)</div><div><br></div><div>E nunca votarei no PSDB/DEM!</div></div>
</div>