<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>