<div dir="ltr">Não seria o caso de você obter esses valores analiticamente, ao invés de usar simulação? Veja a documentação da função pnorm().<div><br></div><div>att,</div><div>FH</div></div><div class="gmail_extra"><br><br>
<div class="gmail_quote">2014-07-28 21:14 GMT-05:00 Marcos Bissoli <span dir="ltr"><<a href="mailto:mbissoli@gmail.com" target="_blank">mbissoli@gmail.com</a>></span>:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
<div dir="ltr">Boa noite a todos,<div><br></div><div>Preciso gerar uma sequência de números aleatórios, de uma distribuição normal. No entanto, a função <rnorm> não permite definir limites mínimo e máximo para o vetor de dados. Assim, criei uma função para tal fim:</div>
<div><br></div><div><div>>Gerador.Aleatorio <- function(n,menor,maior,x,s)</div><div>+{</div><div>+ dados <- rnorm(n,x,s)</div><div>+ while ((min(dados)<menor)|(max(dados)>maior))</div><div>+ {</div><div>
+ dados <- rnorm(n,x,s)</div><div>+ }</div><div>+ dados <- dados</div><div>+ print(dados)</div><div>+}</div><div><br></div><div>Onde:</div><div>- n é o número de valores que quero gerar</div><div>- menor é o limite mínimo do vetor</div>
<div>- maior é o limite máximo do vetor</div><div>- x é a média do vetor</div><div>- s é o desvio padrão do vetor</div><div><br></div><div>O problema é que, como usei um laço ("loop"), o custo de tempo computacional dessa função tornou-a simplesmente impraticável. Até 100 valores eu consegui gerar. Mas agora tentei com 284 e já jantei, fumei, tratei de assuntos familiares e... Nada! E minha demanda, na verdade, é para gerar mais de 1500 valores, simulando estudo epidemiológico observacional, a partir de valores já descritos em um estudo piloto que realizei (ou seja, eu tenho já estimativas de n, mínimo, máximo, média e desvio padrão).</div>
<div><br></div><div>Há alguma função, pacote ou sintaxe que torne isso possível? E qual seriam esses "limites do custo computacional"?</div><div><br></div><div>Grato desde já por qualquer ajuda.</div><div><br></div>
<div>Abraços fraternos,</div><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></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></div>