<html>
  <head>
    <meta http-equiv="Content-Type" content="text/html; charset=UTF-8">
  </head>
  <body text="#000000" bgcolor="#FFFFFF">
    <p>Bom dia Marcus e obrigado pela ajuda, <br>
    </p>
    <p>Mas ainda tem um problema, pois agora com a inserção de 
      ggplot(mapping=aes(x=Income, y=value, color = Gender2)), as linhas
      para o modelo predito ziguezagueiam ao invés de se comportarem
      como curvas. A ideia de utilizar a função gather() é justamente
      funcionar com a função order() quando eu fazia gráficos com a
      função plot() padrão e aparentemente isto esta bem difícil de ser
      resolvido.<br>
    </p>
    <pre class="moz-signature" cols="72">-- 
======================================================================
Alexandre dos Santos
Proteção Florestal 
IFMT - Instituto Federal de Educação, Ciência e Tecnologia de Mato Grosso
Campus Cáceres
Caixa Postal 244
Avenida dos Ramires, s/n
Bairro: Distrito Industrial 
Cáceres - MT                      CEP: 78.200-000
Fone: (+55) 65 99686-6970 (VIVO) (+55) 65 3221-2674 (FIXO)
<a class="moz-txt-link-abbreviated" href="mailto:e-mails:alexandresantosbr@yahoo.com.br">e-mails:alexandresantosbr@yahoo.com.br</a> 
        <a class="moz-txt-link-abbreviated" href="mailto:alexandre.santos@cas.ifmt.edu.br">alexandre.santos@cas.ifmt.edu.br</a> 
Lattes: <a class="moz-txt-link-freetext" href="http://lattes.cnpq.br/1360403201088680">http://lattes.cnpq.br/1360403201088680</a> 
OrcID: orcid.org/0000-0001-8232-6722   -   ResearcherID: A-5790-2016
Researchgate: <a class="moz-txt-link-abbreviated" href="http://www.researchgate.net/profile/Alexandre_Santos10">www.researchgate.net/profile/Alexandre_Santos10</a>                       
LinkedIn: br.linkedin.com/in/alexandre-dos-santos-87961635
Mendeley:www.mendeley.com/profiles/alexandre-dos-santos6/
======================================================================</pre>
    <div class="moz-cite-prefix">Em 19/06/2019 11:08, Marcus Nunes
      escreveu:<br>
    </div>
    <blockquote type="cite"
cite="mid:CA+QGQvuyJFiVnbv51B7QtJWt+Zz7=fTqDjYRD-nHU9_Mv=pscQ@mail.gmail.com">
      <meta http-equiv="content-type" content="text/html; charset=UTF-8">
      <div dir="ltr">O problema está na linha
        <div><br>
        </div>
        <div>  ggplot(mapping=aes(x=type, y=value, color = Gender2)) +<br>
        </div>
        <div><br>
        </div>
        <div>Não faz muito sentido criar um gráfico de dispersão com uma
          variável categórica de apenas um nível como variável
          preditora. Rode o código abaixo que ele deve ficar mais
          próximo do desejado.</div>
        <div><br>
        </div>
        <div>df3  %>%<br>
            tidyr::gather(type, value, Consumption) %>%<br>
            ggplot(mapping=aes(x=Income, y=value, color = Gender2)) +<br>
            geom_smooth(mapping=aes(ymin = lcl, ymax = ucl), stat =
          "identity", colour = "black") +<br>
            geom_point(df,mapping=aes(x=Income, y=Consumption, color =
          Gender2), size = 2) +<br>
            geom_line(mapping=aes(x=Income, y=pred), colour = "black")<br>
        </div>
        <div><br>
        </div>
        <div>Eu achei esquisito que o nível Fem_Her está muito distante
          do nível Male no gráfico, mas como a income dos Male é em
          torno de 20 vezes a dos Fem_Her, então o plot faz sentido.</div>
        <div><br>
        </div>
        <div>Ah, troquei a cor da linha dos valores preditos porque ela
          estava sendo confundida com alguns pontos. Talvez preto não
          seja a melhor opção, então escolha uma cor que se harmonize
          melhor à tua paleta de cores.</div>
        <div><br>
        </div>
        <div>Atenciosamente,<br clear="all">
          <div>
            <div dir="ltr" class="gmail_signature"
              data-smartmail="gmail_signature">
              <div dir="ltr">
                <div>
                  <div dir="ltr">
                    <div>
                      <div dir="ltr">
                        <div>
                          <div dir="ltr">
                            <div dir="ltr">
                              <div dir="ltr">
                                <div dir="ltr">
                                  <div dir="ltr">--</div>
                                  <div dir="ltr">
                                    <div dir="ltr">Marcus Nunes</div>
                                    <div dir="ltr">Professor Adjunto</div>
                                    <div dir="ltr">Universidade Federal
                                      do Rio Grande do Norte</div>
                                    <div dir="ltr">Centro de Ciências
                                      Exatas e da Terra</div>
                                    <div dir="ltr">Departamento de
                                      Estatística</div>
                                    <div dir="ltr">Laboratório de
                                      Estatística Aplicada</div>
                                    <div dir="ltr"><a
                                        href="mailto:marcus@marcusnunes.me"
                                        target="_blank"
                                        moz-do-not-send="true">marcus@marcusnunes.me</a></div>
                                    <div dir="ltr"><a
                                        href="https://marcusnunes.me/"
                                        target="_blank"
                                        moz-do-not-send="true">https://marcusnunes.me/</a></div>
                                    <div dir="ltr"><a
                                        href="http://shiny.estatistica.ccet.ufrn.br"
                                        target="_blank"
                                        moz-do-not-send="true">http://shiny.estatistica.ccet.ufrn.br</a></div>
                                    <div dir="ltr"><br>
                                    </div>
                                  </div>
                                </div>
                              </div>
                            </div>
                          </div>
                        </div>
                      </div>
                    </div>
                  </div>
                </div>
              </div>
            </div>
          </div>
          <br>
        </div>
      </div>
      <br>
      <div class="gmail_quote">
        <div dir="ltr" class="gmail_attr">On Wed, Jun 19, 2019 at 11:45
          AM ASANTOS por (R-br) <<a
            href="mailto:r-br@listas.c3sl.ufpr.br"
            moz-do-not-send="true">r-br@listas.c3sl.ufpr.br</a>>
          wrote:<br>
        </div>
        <blockquote class="gmail_quote" style="margin:0px 0px 0px
          0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">Prezados
          Membros,<br>
          <br>
          Estou tentando plotar sem sucesso o intervalo em um modelo GLM
          de <br>
          Poisson onde fiz a junção de níveis para a variável categórica
          Gender. <br>
          Alguém baseado no CRM abaixo poderia dar uma luz?<br>
          <br>
          #Pacotes<br>
          library(ggplot2)<br>
          library(dplyr)<br>
          library(tidyverse)<br>
          <br>
          #Dados inventados com uma variável resposta (Consumption) e
          duas <br>
          explicativas, sendo uma categórica e outra quali<br>
          Consumption <- c(501, 502, 503, 504, 26, 27, 50, 56, 63,
          60, 72, 93, 78, <br>
          43, 59, 70, 53, 80)<br>
          Gender <- gl(n = 3, k = 6, length = 3*6, labels = c("Male",
          <br>
          "Female","hermaphrodite"), ordered = FALSE)<br>
          Income <- c(5010, 5020, 5030, 5040, 260, 270, 550, 560,
          680, 690, 720, <br>
          550, 560, 680, 690, 720,500,512)<br>
          df3 <- data.frame(Consumption, Gender, Income)<br>
          df3<br>
          <br>
          # GLM de Poisson<br>
          fm1 <- glm(Consumption~Gender+Income, data=df3,
          family=poisson)<br>
          summary(fm1)<br>
          <br>
          # ANOVA do modelo ajustado<br>
          anova(fm1,test="Chi")<br>
          <br>
          #Comparo a variável Gender<br>
          sort(tapply(df3$Consumption,df3$Gender,mean))<br>
          Gender2<-df3$Gender<br>
          levels(Gender2)<br>
          levels(Gender2)[2]<-"Fem_Her"<br>
          levels(Gender2)[3]<-"Fem_Her"<br>
          levels(Gender2)<br>
          fm2<-glm(Consumption~Gender2+Income, data=df3,
          family=poisson)<br>
          anova(fm1,fm2,test="Chi")<br>
          # 0.7824 Female/Hermaphrodite são iguais então eu junto<br>
          <br>
          #Faço a predição sobre o modelo final e dos intervalos de
          confiança<br>
          <br>
          pred <- predict(fm2, type="response", se.fit = TRUE)<br>
          df3 = cbind(df3, pred = pred$fit)<br>
          df3 = cbind(df3, se = pred$se.fit)<br>
          df3 = cbind(df3, ucl=df3$pred + 1.96*df3$se)<br>
          df3 = cbind(df3, lcl=df3$pred - 1.96*df3$se)<br>
          df3 = cbind(df3, Gender2)<br>
          <br>
          df<-df3 %>%<br>
             dplyr::group_by(Income, Gender2) %>%<br>
             dplyr::summarize(Consumption = mean(Consumption, na.rm =
          TRUE))<br>
          df<-as.data.frame(df)<br>
          <br>
          #Faço o plot usando o ggplot2<br>
          df3  %>%<br>
             tidyr::gather(type, value, Consumption) %>%<br>
             ggplot(mapping=aes(x=type, y=value, color = Gender2)) +<br>
                 geom_smooth(mapping=aes(ymin = lcl, ymax = ucl), stat =
          <br>
          "identity") +<br>
             geom_point(df,mapping=aes(x=Income, y=Consumption, color =
          Gender2)) +<br>
             geom_line(mapping=aes(x=Income, y=pred))<br>
          <br>
          #<br>
          <br>
          Obrigado,<br>
          <br>
          Alexandre<br>
          <br>
          -- <br>
======================================================================<br>
          Alexandre dos Santos<br>
          Proteção Florestal<br>
          IFMT - Instituto Federal de Educação, Ciência e Tecnologia de
          Mato Grosso<br>
          Campus Cáceres<br>
          Caixa Postal 244<br>
          Avenida dos Ramires, s/n<br>
          Bairro: Distrito Industrial<br>
          Cáceres - MT                      CEP: 78.200-000<br>
          Fone: (+55) 65 99686-6970 (VIVO) (+55) 65 3221-2674 (FIXO)<br>
          <a href="mailto:e-mails%3Aalexandresantosbr@yahoo.com.br"
            target="_blank" moz-do-not-send="true">e-mails:alexandresantosbr@yahoo.com.br</a><br>
                   <a href="mailto:alexandre.santos@cas.ifmt.edu.br"
            target="_blank" moz-do-not-send="true">alexandre.santos@cas.ifmt.edu.br</a><br>
          Lattes: <a href="http://lattes.cnpq.br/1360403201088680"
            rel="noreferrer" target="_blank" moz-do-not-send="true">http://lattes.cnpq.br/1360403201088680</a><br>
          OrcID: <a href="http://orcid.org/0000-0001-8232-6722"
            rel="noreferrer" target="_blank" moz-do-not-send="true">orcid.org/0000-0001-8232-6722</a> 
           -   ResearcherID: A-5790-2016<br>
          Researchgate: <a
            href="http://www.researchgate.net/profile/Alexandre_Santos10"
            rel="noreferrer" target="_blank" moz-do-not-send="true">www.researchgate.net/profile/Alexandre_Santos10</a><br>
          LinkedIn: <a
            href="http://br.linkedin.com/in/alexandre-dos-santos-87961635"
            rel="noreferrer" target="_blank" moz-do-not-send="true">br.linkedin.com/in/alexandre-dos-santos-87961635</a><br>
          Mendeley:<a
            href="http://www.mendeley.com/profiles/alexandre-dos-santos6/"
            rel="noreferrer" target="_blank" moz-do-not-send="true">www.mendeley.com/profiles/alexandre-dos-santos6/</a><br>
======================================================================<br>
          <br>
          _______________________________________________<br>
          R-br mailing list<br>
          <a href="mailto:R-br@listas.c3sl.ufpr.br" target="_blank"
            moz-do-not-send="true">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" moz-do-not-send="true">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" moz-do-not-send="true">http://www.leg.ufpr.br/r-br-guia</a>)
          e fornea cdigo mnimo reproduzvel.</blockquote>
      </div>
    </blockquote>
  </body>
</html>