<div dir="ltr"><div><div><div><div><div><div><div><div>Prezados...<br><br></div>Executei o comando sugerido pelo Mauricio<br><br><b>mod1<-glm(PoliFar~ClassEco+Educ,family="poisson"(link="log"),data=dmha2)</b><br><br> e estou encontrando a seguinte mensagem de erro<br><b><br>Error in if (any(y < 0)) stop("negative values not allowed for the 'Poisson' family") : valor ausente onde TRUE/FALSE necessário<br>Além disso: Warning message:<br>In Ops.factor(y, 0) : ‘<’ not meaningful for factors</b><br><br></div></div>Meus dados estão no seguinte formato: (<b>data frame dmha2</b>):<br><br></div><b>ClassEco Educ PoliFar<br></b></div><b>A 0-2 S<br></b></div><b>D-E 3-5 N<br>B 12+ S<br></b></div><div><b>A 6-8 S<br></b></div><b>C 9-11 N<br></b></div><b>D-E 0-2 S<br></b><div><div><div><div><div><b>... ... ...</b><br><br>Vocês sabem o que está errado?<br><br></div><div>Att.,<br></div></div></div></div></div></div><div class="gmail_extra"><br clear="all"><div><div class="gmail_signature"><div dir="ltr"><b><br></b><div style="text-align:left"><i><b>Emerson</b></i><br></div></div></div></div>
<br><div class="gmail_quote">Em 19 de abril de 2016 19:03, Emerson Cotta Bodevan <span dir="ltr"><<a href="mailto:bodevan.ec@gmail.com" target="_blank">bodevan.ec@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>Obrigado Mauricio.<br><br></div><div>Vou executar aqui.<br><br></div><div>Att.,<br></div><div><br></div><div class="gmail_extra"><br clear="all"><div><div><div dir="ltr"><b><br></b><div style="text-align:left"><i><b>Emerson</b></i><br></div></div></div></div><div><div class="h5">
<br><div class="gmail_quote">Em 19 de abril de 2016 19:01, Mauricio Cardeal <span dir="ltr"><<a href="mailto:mcardeal2010@gmail.com" target="_blank">mcardeal2010@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 text="#000000" bgcolor="#FFFFFF">
Emerson, você pode tentar assim:<span><br>
<br>
modelo=glm(variavel dependente ~ variaveis
independentes,family="poisson"(link="log"))<br>
exp(modelo$coefficients)<br>
summary(modelo)<br>
<br></span>
Para os intervalos de confiança use a função confint:<br>
<br>
exp(confint(modelo)) # 95% CI for exponentiated coefficients<span><font color="#888888"><br>
<br>
Mauricio Cardeal<br>
UFBA</font></span><div><div><br>
<br>
<div>Em 19/04/2016 18:37, Emerson Cotta
Bodevan escreveu:<br>
</div>
<blockquote type="cite">
<div dir="ltr">
<div>
<div>
<div>
<div>Olá Marco!<br>
<br>
</div>
Obrigado pela dica. Depois de ler as referências e
sugestões, vi que (apesar de minha resposta ser
dicotômica) , como meu estudo é transversal, o melhor é
usar a Regressão de Poisson com variância Robusta.<br>
<br>
</div>
A dica do Maurício Cardeal. Agora... como construo o IC da
Razão de Prevalência?<br>
</div>
<br>
</div>
Agradeço a todos.<br>
</div>
<div class="gmail_extra"><br clear="all">
<div>
<div>
<div dir="ltr"><b><br>
</b>
<div style="text-align:left"><i><b>Emerson</b></i><br>
</div>
</div>
</div>
</div>
<br>
<div class="gmail_quote">Em 19 de abril de 2016 17:22, Marco
Nunes <span dir="ltr"><<a href="mailto:nunes.ma@outlook.com" target="_blank">nunes.ma@outlook.com</a>></span>
escreveu:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
<div>
<div dir="ltr">Olá Emerson<br>
<br>
O uso da razão de prevalência na regressão logística já
foi reapondido.<br>
Para o cálculo de razão de prevalência (assim como risco
relatico e razão de odds também) utilizo o pacote
"epiR".<br>
Encaminho um tutorial abaixo.<br>
<br>
<br>
# Exemplo de utilização do epi.2by2 para calculo da
razão de prevalencia<br>
<br>
library(epiR)<br>
<br>
dat <- matrix(c(13,2163,5,3349), nrow = 2, byrow =
TRUE)<br>
rownames(dat) <- c("DF+", "DF-")<br>
colnames(dat) <- c("FUS+", "FUS-")<br>
dat<br>
epi.2by2(dat = as.table(dat), method =
"cross.sectional", conf.level = 0.95, units = 100,
homogeneity = "breslow.day", outcome = "as.columns")<br>
<br>
<br>
Prof. Dr. Marco Antonio Prado Nunes<br>
Departamento de Medicina/UFS<br>
<br>
<div>
<hr>Date: Sat, 9 Apr 2016 15:49:02 -0300<br>
From: <a href="mailto:bodevan.ec@gmail.com" target="_blank">bodevan.ec@gmail.com</a><br>
To: <a href="mailto:r-br@listas.c3sl.ufpr.br" target="_blank">r-br@listas.c3sl.ufpr.br</a><br>
Subject: Re: [R-br] Razão de prevalência
<div>
<div><br>
<br>
<div dir="ltr">
<div>
<div>Obrigado César.<br>
<br>
</div>
Vou verificar.<br>
<br>
</div>
Abraço,<br>
</div>
<div><br clear="all">
<div>
<div>
<div dir="ltr"><b><br>
</b>
<div style="text-align:left"><i><b>Emerson</b></i><br>
</div>
</div>
</div>
</div>
<br>
<div>Em 8 de abril de 2016 14:30, Cesar Rabak <span dir="ltr"><<a href="mailto:cesar.rabak@gmail.com" target="_blank">cesar.rabak@gmail.com</a>></span>
escreveu:<br>
<blockquote style="border-left:1px #ccc solid;padding-left:1ex">
<div dir="ltr">Olá Emerson,
<div><br>
</div>
<div>Embora a referência apresentada pelo
Leonardo seja recente e interessante,
pode acontecer que devido a exigências
do veículo que você publicará,
orientação ou chefia, o cálculo tenha
que seguir com a RL regular para os IC
da RP.</div>
<div><br>
</div>
<div>Nesse caso, recomendo a leitura de um
artigo mais "clássico" e que tem boas
referências e código e R: <a href="http://www.scielo.br/pdf/csp/v31n3/0102-311X-csp-31-03-00487.pdf" target="_blank"></a><a href="http://www.scielo.br/pdf/csp/v31n3/0102-311X-csp-31-03-00487.pdf" target="_blank">http://www.scielo.br/pdf/csp/v31n3/0102-311X-csp-31-03-00487.pdf</a></div>
<div><br>
</div>
<div>HTH</div>
<div>--</div>
<div>Cesar Rabak</div>
</div>
<div><br>
<div>
<div>
<div>2016-04-07 15:21 GMT-03:00
Emerson Cotta Bodevan <span dir="ltr"><<a href="mailto:bodevan.ec@gmail.com" target="_blank"></a><a href="mailto:bodevan.ec@gmail.com" target="_blank">bodevan.ec@gmail.com</a>></span>:<br>
</div>
</div>
<blockquote style="border-left:1px #ccc solid;padding-left:1ex">
<div>
<div>
<div dir="ltr">
<div>
<div>
<div>
<div>Prezados, boa tarde.<br>
<br>
</div>
Estou usando regressão
logística, tendo como
resposta a ocorrência ou
não de uma doença comum
(alta prevalência).<br>
<br>
</div>
Como calcular a <b>razão de
prevalência</b> e seus
intervalos de confiança?
Estudo transversal.<br>
<br>
</div>
Algum pacote adequado?<br>
<br>
</div>
Agradeço qualquer ajuda.<span><font color="#888888"><br clear="all">
<div>
<div>
<div>
<div>
<div>
<div>
<div>
<div>
<div>
<div dir="ltr"><b><br>
</b>
<div style="text-align:left"><i><b>Emerson</b></i><br>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</font></span></div>
<br>
</div>
</div>
<span>_______________________________________________<br>
R-br mailing list<br>
<a href="mailto:R-br@listas.c3sl.ufpr.br" target="_blank">R-br@listas.c3sl.ufpr.br</a><br>
<a href="https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br" rel="noreferrer" 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" rel="noreferrer" target="_blank"></a><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>
</span></blockquote>
</div>
<br>
</div>
<br>
_______________________________________________<br>
R-br mailing list<br>
<a href="mailto:R-br@listas.c3sl.ufpr.br" target="_blank">R-br@listas.c3sl.ufpr.br</a><br>
<a href="https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br" rel="noreferrer" 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" rel="noreferrer" target="_blank"></a><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>
<br>
</div>
</div>
_______________________________________________
R-br mailing list
<a href="mailto:R-br@listas.c3sl.ufpr.br" target="_blank">R-br@listas.c3sl.ufpr.br</a>
<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>
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.</div>
</div>
</div>
<br>
_______________________________________________<br>
R-br mailing list<br>
<a href="mailto:R-br@listas.c3sl.ufpr.br" target="_blank">R-br@listas.c3sl.ufpr.br</a><br>
<a href="https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br" rel="noreferrer" 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" rel="noreferrer" 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>
<br>
<fieldset></fieldset>
<br>
<pre>_______________________________________________
R-br mailing list
<a href="mailto:R-br@listas.c3sl.ufpr.br" target="_blank">R-br@listas.c3sl.ufpr.br</a>
<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>
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.</pre>
</blockquote>
<br>
</div></div></div>
<br>_______________________________________________<br>
R-br mailing list<br>
<a href="mailto:R-br@listas.c3sl.ufpr.br" target="_blank">R-br@listas.c3sl.ufpr.br</a><br>
<a href="https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br" rel="noreferrer" 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" rel="noreferrer" 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></div></div></div>
</blockquote></div><br></div>