<div dir="ltr">Obrigado pela dica Walmes. O que você sugere para eu obter esses dados? Simulação MC?<div><br></div><div>Eu gostaria de colocar isso dentro de uma função para que eu possa obter esse resultado para um conjunto grande de valores.<div>
<br></div><div>Segue um exemplo:<div><br></div><div><div>idade=c(40,45,50,55,60,65,60,30,35,22)</div><div>sd=c(10.2,15.1,6,6.5,10,12.5,8,9,9.5,10)</div><div>npart=c(1000,100,1500,10000,500,200,300,400,500,150)</div><div>dados = cbind(as.data.frame(idade),as.data.frame(sd),as.data.frame(npart))</div>
</div></div><div><br></div><div>Ficaria assim a função?</div><div><br></div><div><div style="font-family:arial,sans-serif;font-size:13px">prop <- function(x){</div><div style="font-family:arial,sans-serif;font-size:13px">
for (i in 1:nrow(x))</div><div style="font-family:arial,sans-serif;font-size:13px">z[i]=<span style="font-family:'courier new',monospace">100*pnorm(45, </span>x$idade[i]<span style="font-family:'courier new',monospace">, </span>x$sd[i]<span style="font-family:'courier new',monospace">, lower.tail=FALSE)</span></div>
<div style="font-family:arial,sans-serif;font-size:13px">resultado <- z<br></div><div style="font-family:arial,sans-serif;font-size:13px">print(resultado)}</div></div><div><br></div><div>prop(dados)<br></div><div><br></div>
<div>Mas ainda não estou obtendo os valores o que preciso modificar nessa função.</div></div></div><div class="gmail_extra"><br><br><div class="gmail_quote">Em 30 de novembro de 2013 12:29, 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 você tem média e desvio-padrão e acredita que X, sua variável aleatória, tem esses valores como parâmetros de uma distribuição normal, você não precisa simular e calcular frequências, use direto as probabilidades teóricas e multiplique pelo número de pessoas. Ou seja,<br>

<br><span style="font-family:courier new,monospace">## Seja X ~ Normal(mean=2, dp=1).<br>## Qual P(X>2.5)?<br>## Se n = 100, qual esperado para n*(X>2.5)? E = n*p<br><br>## por simulação (sujeito à erro Monte Carlo),<br>

## vai convergir quando n for para infinito.<br>x <- rnorm(100, 2, 1)<br>table(x>2.5)<br><br>## por teoria<br>100*pnorm(2.5, 2, 1, lower.tail=FALSE)<br></span><br></div><div class="gmail_default" style="font-family:trebuchet ms,sans-serif">

Você tá assumindo uma distribuição teórica a partir da qual pode-se usar, portanto, os resultados teóricos mas tá resolvendo o problema na força bruta, simulando da variável e obtendo a distribuição de frequência. É como se você tivesse uma moto pronta para andar mas resolveu retirar o seu motor e adaptar em uma bicicleta. Vai te custar tempo, não vai funcionar igual, não vai ser tão rápido nem tão seguro, quem ver vai achar esquisito e te encher de perguntas, muito embora você vai de A para B nessa bicicleta. Mas você não precisa fazer isso. Vá de moto.<br>

</div><div class="gmail_default" style="font-family:trebuchet ms,sans-serif"><br>À disposição.<span class="HOEnZb"><font color="#888888"><br></font></span></div><span class="HOEnZb"><font color="#888888"><div class="gmail_default" style="font-family:trebuchet ms,sans-serif">
Walmes.</div><br></font></span></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">Sérgio Henrique Almeida da Silva Junior<br>Doutorando em Epidemiologia em Saúde Pública<br>Escola Nacional de Saúde Pública Sérgio Arouca - ENSP/FIOCRUZ<br><a href="http://lattes.cnpq.br/1611345552843383" target="_blank">http://lattes.cnpq.br/1611345552843383</a> <br>
Tel: (21) 68463637<br><div><a href="http://www.linkedin.com/profile/view?id=250437145&trk=tab_pro" target="_blank">http://www.linkedin.com/profile/view?id=250437145&trk=tab_pro</a><br></div></div>
</div>