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