<html>
<head>
<meta content="text/html; charset=utf-8" http-equiv="Content-Type">
</head>
<body text="#000000" bgcolor="#FFFFFF">
<tt>Emerson, segue um exemplo sem erro:<br>
<br>
</tt><tt>>
a1=glm(ce~fx+pp+di+sx+tb+hf,family="poisson"(link="log"))</tt><tt><br>
</tt><tt>> summary(a1,cor=F)</tt><tt><br>
</tt><tt><br>
</tt><tt>Call:</tt><tt><br>
</tt><tt>glm(formula = ce ~ fx + pp + di + sx + tb + hf, family =
poisson(link = "log"))</tt><tt><br>
</tt><tt><br>
</tt><tt>Deviance Residuals: </tt><tt><br>
</tt><tt>  Min     1Q Median     3Q    Max </tt><tt><br>
</tt><tt>-1,406Â -0,864Â -0,703Â Â 0,593Â Â 1,576Â </tt><tt><br>
</tt><tt><br>
</tt><tt>Coefficients:</tt><tt><br>
</tt><tt>           Estimate Std. Error z value       Â
Pr(>|z|)Â Â Â </tt><tt><br>
</tt><tt>(Intercept)Â Â -1,399Â Â Â Â Â 0,196Â Â -7,14 0,00000000000094
***</tt><tt><br>
</tt><tt>fx            0,455     0,211   2,16           0,031 *Â
</tt><tt><br>
</tt><tt>pp            0,518     0,263   1,97           0,049 *Â
</tt><tt><br>
</tt><tt>di            0,414     0,207   2,00           0,046 *Â
</tt><tt><br>
</tt><tt>sx            0,287     0,202   1,42           0,156  Â
</tt><tt><br>
</tt><tt>tb           -1,177     1,009  -1,17           0,243  Â
</tt><tt><br>
</tt><tt>hf           -0,338     0,421  -0,80           0,423  Â
</tt><tt><br>
</tt><tt>---</tt><tt><br>
</tt><tt>Signif. codes: 0 ‘***’ 0,001 ‘**’ 0,01 ‘*’ 0,05 ‘.’ 0,1 ‘
’ 1</tt><tt><br>
</tt><tt><br>
</tt><tt>(Dispersion parameter for poisson family taken to be 1)</tt><tt><br>
</tt><tt><br>
</tt><tt>Â Â Â Null deviance: 183,26Â on 250Â degrees of freedom</tt><tt><br>
</tt><tt>Residual deviance: 160,98 on 244 degrees of freedom ####
Aqui Residual deviance não é maior do que os graus de liberdade
(sugere não haver overdispersion)</tt><tt><br>
</tt><tt>Â (16 observations deleted due to missingness)</tt><tt><br>
</tt><tt>AIC: 383</tt><tt><br>
</tt><tt><br>
</tt><tt>Number of Fisher Scoring iterations: 5</tt><tt><br>
</tt><tt><br>
#Razões de prevalências<br>
<br>
</tt><tt>> exp(a1$coefficients)</tt><tt><br>
</tt><tt>(Intercept)         fx         pp         di        Â
sx         tb         hf </tt><tt><br>
</tt><tt>Â Â Â 0,24673Â Â Â Â 1,57594Â Â Â Â 1,67850Â Â Â Â 1,51360Â Â Â Â
1,33183Â Â Â Â 0,30818Â Â Â Â 0,71343 </tt><tt><br>
</tt><tt>> exp(confint(a1)) # 95% CI for exponentiated
coefficients</tt><tt><br>
</tt><tt>Waiting for profiling to be done...</tt><tt><br>
</tt><tt>Â Â Â Â Â Â Â Â Â Â Â Â Â Â 2,5 %Â 97,5 %</tt><tt><br>
</tt><tt>(Intercept) 0,165063 0,35644</tt><tt><br>
</tt><tt>fx         1,035317 2,36868</tt><tt><br>
</tt><tt>pp         0,975754 2,74838</tt><tt><br>
</tt><tt>di         1,015014 2,29550</tt><tt><br>
</tt><tt>sx         0,890559 1,97234</tt><tt><br>
</tt><tt>tb         0,017425 1,39886</tt><tt><br>
</tt><tt>hf         0,278106 1,49427</tt><tt><br>
</tt><tt><br>
<br>
# Verificando fator de inflação das variancia<br>
<br>
> vif(a1)</tt><tt><br>
</tt><tt>   fx    pp    di    sx    tb    hf </tt><tt><br>
</tt><tt>1,0911 1,0722 1,0131 1,0238 1,0084 1,0040 </tt><tt><br>
</tt><tt>> sqrt(vif(a1))</tt><tt><br>
</tt><tt>   fx    pp    di    sx    tb    hf </tt><tt><br>
</tt><tt>1,0446 1,0355 1,0065 1,0118 1,0042 1,0020 </tt><tt><br>
</tt><tt>> sqrt(vif(a1)) > 2</tt><tt><br>
</tt><tt>  fx   pp   di   sx   tb   hf </tt><tt><br>
</tt><tt>FALSE FALSE FALSE FALSE FALSE FALSE <br>
<br>
# Análise de residuos studentizados<br>
</tt><tt></tt><tt>> #residuos</tt><tt><br>
</tt><tt>> library(MASS)</tt><tt><br>
</tt><tt>> sresid <- studres(a1)</tt><tt><br>
</tt><tt>> mean(sresid)</tt><tt><br>
</tt><tt>[1] 0,00022105</tt><tt><br>
</tt><tt>> var(sresid)</tt><tt><br>
</tt><tt>[1] 1,0093</tt><tt><br>
</tt><tt>> <br>
<br>
#variável dependente<br>
> tab(ce)<br>
 freq.abs freq.rel<br>
0Â Â Â Â Â 149Â Â Â Â 58,9<br>
1Â Â Â Â Â 104Â Â Â Â 41,1<br>
Total:Â 253 <br>
<br>
Mauricio Cardeal<br>
UFBA<br>
<br>
</tt>
<div class="moz-cite-prefix">Em 19/04/2016 22:00, Emerson Cotta
Bodevan escreveu:<br>
</div>
<blockquote
cite="mid:CAP5i-syDp_PAekeqiKUtpoB3_xbmFKtjz_Sf28RVu9dshLd72w@mail.gmail.com"
type="cite">
<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 moz-do-not-send="true"
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
moz-do-not-send="true"
href="mailto:mcardeal2010@gmail.com"
target="_blank"><a class="moz-txt-link-abbreviated" href="mailto:mcardeal2010@gmail.com">mcardeal2010@gmail.com</a></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
moz-do-not-send="true"
href="mailto:nunes.ma@outlook.com"
target="_blank"><a class="moz-txt-link-abbreviated" href="mailto:nunes.ma@outlook.com">nunes.ma@outlook.com</a></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
moz-do-not-send="true"
href="mailto:bodevan.ec@gmail.com"
target="_blank"><a class="moz-txt-link-abbreviated" href="mailto:bodevan.ec@gmail.com">bodevan.ec@gmail.com</a></a><br>
To: <a
moz-do-not-send="true"
href="mailto:r-br@listas.c3sl.ufpr.br"
target="_blank"><a class="moz-txt-link-abbreviated" href="mailto:r-br@listas.c3sl.ufpr.br">r-br@listas.c3sl.ufpr.br</a></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
moz-do-not-send="true"
href="mailto:cesar.rabak@gmail.com" target="_blank"><a class="moz-txt-link-abbreviated" href="mailto:cesar.rabak@gmail.com">cesar.rabak@gmail.com</a></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
moz-do-not-send="true"
href="http://www.scielo.br/pdf/csp/v31n3/0102-311X-csp-31-03-00487.pdf"
target="_blank"><a class="moz-txt-link-freetext" href="http://www.scielo.br/pdf/csp/v31n3/0102-311X-csp-31-03-00487.pdf">http://www.scielo.br/pdf/csp/v31n3/0102-311X-csp-31-03-00487.pdf</a></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
moz-do-not-send="true" href="mailto:bodevan.ec@gmail.com"
target="_blank"><a class="moz-txt-link-abbreviated" href="mailto:bodevan.ec@gmail.com">bodevan.ec@gmail.com</a></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
moz-do-not-send="true"
href="mailto:R-br@listas.c3sl.ufpr.br" target="_blank"><a class="moz-txt-link-abbreviated" href="mailto:R-br@listas.c3sl.ufpr.br">R-br@listas.c3sl.ufpr.br</a></a><br>
<a
moz-do-not-send="true"
href="https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br"
rel="noreferrer"
target="_blank"><a class="moz-txt-link-freetext" href="https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br">https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br</a></a><br>
Leia o guia de
postagem (<a
moz-do-not-send="true"
href="http://www.leg.ufpr.br/r-br-guia" target="_blank"><a class="moz-txt-link-freetext" href="http://www.leg.ufpr.br/r-br-guia">http://www.leg.ufpr.br/r-br-guia</a></a>)
e forneça
código mÃnimo
reproduzÃvel.<br>
</span></blockquote>
</div>
<br>
</div>
<br>
_______________________________________________<br>
R-br mailing list<br>
<a
moz-do-not-send="true"
href="mailto:R-br@listas.c3sl.ufpr.br" target="_blank"><a class="moz-txt-link-abbreviated" href="mailto:R-br@listas.c3sl.ufpr.br">R-br@listas.c3sl.ufpr.br</a></a><br>
<a
moz-do-not-send="true"
href="https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br"
rel="noreferrer"
target="_blank"><a class="moz-txt-link-freetext" href="https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br">https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br</a></a><br>
Leia o guia de
postagem (<a
moz-do-not-send="true"
href="http://www.leg.ufpr.br/r-br-guia" target="_blank"><a class="moz-txt-link-freetext" href="http://www.leg.ufpr.br/r-br-guia">http://www.leg.ufpr.br/r-br-guia</a></a>)
e forneça código
mÃnimo
reproduzÃvel.<br>
</blockquote>
</div>
<br>
</div>
<br>
</div>
</div>
_______________________________________________
R-br mailing list <a
moz-do-not-send="true"
href="mailto:R-br@listas.c3sl.ufpr.br"
target="_blank"><a class="moz-txt-link-abbreviated" href="mailto:R-br@listas.c3sl.ufpr.br">R-br@listas.c3sl.ufpr.br</a></a>
<a moz-do-not-send="true"
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
moz-do-not-send="true"
href="http://www.leg.ufpr.br/r-br-guia"
target="_blank"><a class="moz-txt-link-freetext" href="http://www.leg.ufpr.br/r-br-guia">http://www.leg.ufpr.br/r-br-guia</a></a>)
e forne�a c�digo m�nimo
reproduz�vel.</div>
</div>
</div>
<br>
_______________________________________________<br>
R-br mailing list<br>
<a moz-do-not-send="true"
href="mailto:R-br@listas.c3sl.ufpr.br"
target="_blank">R-br@listas.c3sl.ufpr.br</a><br>
<a moz-do-not-send="true"
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
moz-do-not-send="true"
href="http://www.leg.ufpr.br/r-br-guia"
rel="noreferrer" target="_blank"><a class="moz-txt-link-freetext" href="http://www.leg.ufpr.br/r-br-guia">http://www.leg.ufpr.br/r-br-guia</a></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 moz-do-not-send="true" href="mailto:R-br@listas.c3sl.ufpr.br" target="_blank">R-br@listas.c3sl.ufpr.br</a>
<a moz-do-not-send="true" 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 moz-do-not-send="true" 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 moz-do-not-send="true"
href="mailto:R-br@listas.c3sl.ufpr.br"
target="_blank">R-br@listas.c3sl.ufpr.br</a><br>
<a moz-do-not-send="true"
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
moz-do-not-send="true"
href="http://www.leg.ufpr.br/r-br-guia"
rel="noreferrer" target="_blank"><a class="moz-txt-link-freetext" href="http://www.leg.ufpr.br/r-br-guia">http://www.leg.ufpr.br/r-br-guia</a></a>)
e forneça código mÃnimo reproduzÃvel.<br>
</blockquote>
</div>
<br>
</div>
</div>
</div>
</div>
</blockquote>
</div>
<br>
</div>
<br>
<fieldset class="mimeAttachmentHeader"></fieldset>
<br>
<pre wrap="">_______________________________________________
R-br mailing list
<a class="moz-txt-link-abbreviated" href="mailto:R-br@listas.c3sl.ufpr.br">R-br@listas.c3sl.ufpr.br</a>
<a class="moz-txt-link-freetext" href="https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br">https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br</a>
Leia o guia de postagem (<a class="moz-txt-link-freetext" href="http://www.leg.ufpr.br/r-br-guia">http://www.leg.ufpr.br/r-br-guia</a>) e forne�a c�digo m�nimo reproduz�vel.</pre>
</blockquote>
<br>
</body>
</html>