<html><head></head><body><div class="ydp5556b043yahoo-style-wrap" style="font-family: Helvetica Neue, Helvetica, Arial, sans-serif; font-size: 13px;"><div><div><br></div><div>Não funcionou a função prediction. O engraçado é quem nem mensagem de erro mostra.</div><div><br></div><div><br></div><div><span><div>lixo$previsao <-predict(simples, newdata=lixo,interval = "prediction")</div><div><span><span style="color: rgb(0, 0, 0); font-family: Helvetica Neue, Helvetica, Arial, sans-serif;">lixo$previsao <-predict(simples, newdata=lixo,interval = "prediction", level=0.95)</span></span><br></div><div><br></div><div><br></div><div>segue meu codigo,</div><div><br></div><div><br></div><div><span><div>attach(waste)</div><div>names(waste)</div><div>simples <-lm(waste$wastee ~ restrnts)</div><div>simples</div><div>ajuste<-fitted(simples)</div><div>waste$ajuste <-ajuste</div><div>erro <- (waste$ajuste - waste$wastee)</div><div>waste$erro<-erro</div><div>prova <- waste$ajuste-waste$erro</div><div>waste$prova <-prova</div><div>head(waste)</div><div><br></div><div>y<- (0.14685 + 0.01014 * waste$restrnts)</div><div>waste$y<-y</div><div><br></div><div>plot(wastee ~ ajuste)</div><div><br></div><div>head(waste)</div><div><br></div><div><br></div><div>#Prevendo em novos dados.</div><div><br></div><div>library(readxl)</div><div>lixo <- read_excel("Dieretorio/lixo.xlsx")</div><div>View(lixo)</div><div><br></div><div>attach(lixo)</div><div><br></div><div>lixo$previsao <-predict(simples, newdata=lixo,interval="prediction",level=0.95)</div><div>lixo</div></span><br></div></span><br></div><div>eu consigo fazer a previsão, mas preciso do intervalo de confiança para previsão.</div><div><br></div><div class="ydp5556b043signature"><div><font color="#0000ff">Edmar</font></div></div></div>
<div>nsegue </div><div><br></div>
</div><div id="yahoo_quoted_4860631892" class="yahoo_quoted">
<div style="font-family:'Helvetica Neue', Helvetica, Arial, sans-serif;font-size:13px;color:#26282a;">
<div>
Em segunda-feira, 10 de dezembro de 2018 12:00:18 BRST, <r-br-request@listas.c3sl.ufpr.br> escreveu:
</div>
<div><br></div>
<div><br></div>
<div><div dir="ltr">Enviar submissões para a lista de discussão R-br para <br></div><div dir="ltr"> <a ymailto="mailto:r-br@listas.c3sl.ufpr.br" href="mailto:r-br@listas.c3sl.ufpr.br">r-br@listas.c3sl.ufpr.br</a><br></div><div dir="ltr"><br></div><div dir="ltr">Para se cadastrar ou descadastrar via WWW, visite o endereço<br></div><div dir="ltr"> <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><br></div><div dir="ltr">ou, via email, envie uma mensagem com a palavra 'help' no assunto ou<br></div><div dir="ltr">corpo da mensagem para <br></div><div dir="ltr"> <a ymailto="mailto:r-br-request@listas.c3sl.ufpr.br" href="mailto:r-br-request@listas.c3sl.ufpr.br">r-br-request@listas.c3sl.ufpr.br</a><br></div><div dir="ltr"><br></div><div dir="ltr">Você poderá entrar em contato com a pessoa que gerencia a lista pelo<br></div><div dir="ltr">endereço<br></div><div dir="ltr"> <a ymailto="mailto:r-br-owner@listas.c3sl.ufpr.br" href="mailto:r-br-owner@listas.c3sl.ufpr.br">r-br-owner@listas.c3sl.ufpr.br</a><br></div><div dir="ltr"><br></div><div dir="ltr">Quando responder, por favor edite sua linha Assunto assim ela será<br></div><div dir="ltr">mais específica que "Re: Contents of R-br digest..."<br></div><div dir="ltr"><br></div><div dir="ltr"><br></div><div dir="ltr">Tópicos de Hoje:<br></div><div dir="ltr"><br></div><div dir="ltr"> 1. Re: IC para previsao (Elias T. Krainski)<br></div><div dir="ltr"> 2. Problema para plotar IC para modelo glm Gamma usando<br></div><div dir="ltr"> stat_smooth() no ggplot2 (ASANTOS)<br></div><div dir="ltr"><br></div><div dir="ltr"><br></div><div dir="ltr">----------------------------------------------------------------------<br></div><div dir="ltr"><br></div><div dir="ltr">Message: 1<br></div><div dir="ltr">Date: Sun, 9 Dec 2018 16:54:45 +0000 (UTC)<br></div><div dir="ltr">From: "Elias T. Krainski" <<a ymailto="mailto:eliaskrainski@yahoo.com.br" href="mailto:eliaskrainski@yahoo.com.br">eliaskrainski@yahoo.com.br</a>><br></div><div dir="ltr">To: a lista Brasileira oficial de discussão do programa R.<br></div><div dir="ltr"> <<a ymailto="mailto:r-br@listas.c3sl.ufpr.br" href="mailto:r-br@listas.c3sl.ufpr.br">r-br@listas.c3sl.ufpr.br</a>><br></div><div dir="ltr">Subject: Re: [R-br] IC para previsao<br></div><div dir="ltr">Message-ID: <<a ymailto="mailto:815218723.1276327.1544374485235@mail.yahoo.com" href="mailto:815218723.1276327.1544374485235@mail.yahoo.com">815218723.1276327.1544374485235@mail.yahoo.com</a>><br></div><div dir="ltr">Content-Type: text/plain; charset="utf-8"<br></div><div dir="ltr"><br></div><div dir="ltr"><br></div><div dir="ltr">Coloca "prediction" <br></div><div dir="ltr"> <br></div><div dir="ltr"> On Fri, Dec 7, 2018 at 22:39, Edmar Caldas por (R-br)<<a ymailto="mailto:r-br@listas.c3sl.ufpr.br" href="mailto:r-br@listas.c3sl.ufpr.br">r-br@listas.c3sl.ufpr.br</a>> wrote: Preciso de uma ajuda para criar o IC para ajuste do modelo de regressão linear simples.<br></div><div dir="ltr">não dá nenhuma mensagem de erro, cria somente a previsao do modelo em novos dados.<br></div><div dir="ltr"><br></div><div dir="ltr">lixo$previsao <-predict(simples, newdata=lixo,interval="confidence",level=0.95)lixo<br></div><div dir="ltr"><br></div><div dir="ltr"><br></div><div dir="ltr">Edmar_______________________________________________<br></div><div dir="ltr">R-br mailing list<br></div><div dir="ltr"><a ymailto="mailto:R-br@listas.c3sl.ufpr.br" href="mailto:R-br@listas.c3sl.ufpr.br">R-br@listas.c3sl.ufpr.br</a><br></div><div dir="ltr"><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><br></div><div dir="ltr">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. <br></div><div dir="ltr">-------------- Próxima Parte ----------<br></div><div dir="ltr">Um anexo em HTML foi limpo...<br></div><div dir="ltr">URL: <<a href="http://listas.inf.ufpr.br/pipermail/r-br/attachments/20181209/9931f2c9/attachment-0001.html" target="_blank">http://listas.inf.ufpr.br/pipermail/r-br/attachments/20181209/9931f2c9/attachment-0001.html</a>><br></div><div dir="ltr">-------------- Próxima Parte ----------<br></div><div dir="ltr">Um texto embutido e sem conjunto de caracteres especificado foi limpo...<br></div><div dir="ltr">Nome: Untitled<br></div><div dir="ltr">Url: <<a href="http://listas.inf.ufpr.br/pipermail/r-br/attachments/20181209/9931f2c9/attachment-0001.ksh" target="_blank">http://listas.inf.ufpr.br/pipermail/r-br/attachments/20181209/9931f2c9/attachment-0001.ksh</a>><br></div><div dir="ltr"><br></div><div dir="ltr">------------------------------<br></div><div dir="ltr"><br></div><div dir="ltr">Message: 2<br></div><div dir="ltr">Date: Mon, 10 Dec 2018 01:47:12 -0300<br></div><div dir="ltr">From: ASANTOS <<a ymailto="mailto:alexandre.santos@cas.ifmt.edu.br" href="mailto:alexandre.santos@cas.ifmt.edu.br">alexandre.santos@cas.ifmt.edu.br</a>><br></div><div dir="ltr">To: "<a ymailto="mailto:r-br@listas.c3sl.ufpr.br" href="mailto:r-br@listas.c3sl.ufpr.br">r-br@listas.c3sl.ufpr.br</a>" <<a ymailto="mailto:r-br@listas.c3sl.ufpr.br" href="mailto:r-br@listas.c3sl.ufpr.br">r-br@listas.c3sl.ufpr.br</a>><br></div><div dir="ltr">Subject: [R-br] Problema para plotar IC para modelo glm Gamma usando<br></div><div dir="ltr"> stat_smooth() no ggplot2<br></div><div dir="ltr">Message-ID: <<a ymailto="mailto:8d326872-10b9-06f7-cc1b-564caf54f3e7@cas.ifmt.edu.br" href="mailto:8d326872-10b9-06f7-cc1b-564caf54f3e7@cas.ifmt.edu.br">8d326872-10b9-06f7-cc1b-564caf54f3e7@cas.ifmt.edu.br</a>><br></div><div dir="ltr">Content-Type: text/plain; charset=utf-8; format=flowed<br></div><div dir="ltr"><br></div><div dir="ltr">Prezados Membros do R-br,<br></div><div dir="ltr"><br></div><div dir="ltr"> Estou tentando plotar um modelo glm Gamma com stat_smooth() no <br></div><div dir="ltr">ggplot2 sem sucesso. Minha ideia é representar o intervalo de confiança <br></div><div dir="ltr">para um modelo onde realizei a junção de dois níveis para o fator <br></div><div dir="ltr">Feature. Segundo exemplo fictício do meu problema:<br></div><div dir="ltr"><br></div><div dir="ltr">#Simulação de um banco de dados, onde a variável resposta Production tem <br></div><div dir="ltr">distribuição Gamma<br></div><div dir="ltr">set.seed(123)<br></div><div dir="ltr">d<-NULL<br></div><div dir="ltr">N<-50<br></div><div dir="ltr">d$Production <- rgamma(N,10)<br></div><div dir="ltr">d$Feature <- ifelse(d$Production >7 & d$Production<10, <br></div><div dir="ltr">c("green"),ifelse(d$Production>11,<br></div><div dir="ltr">c("red"), c("blue")))<br></div><div dir="ltr">d$Temp<-rnorm(N,20,5)<br></div><div dir="ltr">d<-as.data.frame(d)<br></div><div dir="ltr">#<br></div><div dir="ltr"><br></div><div dir="ltr"># Ajusto um modelo de Gamma completo<br></div><div dir="ltr">mG<- glm(Production~ Feature + Temp, family= Gamma, data = d)<br></div><div dir="ltr">summary(mG)<br></div><div dir="ltr">anova(mG,test="Chi")<br></div><div dir="ltr"><br></div><div dir="ltr"># Comparo os níveis do fator Feature<br></div><div dir="ltr">library(multcomp)<br></div><div dir="ltr">PW<-summary(glht(mG, linfct = mcp(Feature= "Tukey")))<br></div><div dir="ltr">PW<br></div><div dir="ltr">cld(PW)<br></div><div dir="ltr"><br></div><div dir="ltr"># Faço de conta que os níveis green = blue e realizo um novo ajuste<br></div><div dir="ltr">Feature2<-d$Feature<br></div><div dir="ltr">levels(Feature2)<br></div><div dir="ltr">levels(Feature2)[1]<-"blue&green"<br></div><div dir="ltr">levels(Feature2)[2]<-"blue&green"<br></div><div dir="ltr">levels(Feature2)<br></div><div dir="ltr">d$Feature2<-Feature2<br></div><div dir="ltr">mG2<- glm(Production~ Feature2 + Temp, family= Gamma, data = d)<br></div><div dir="ltr"><br></div><div dir="ltr"># Crio os valores de predição<br></div><div dir="ltr">pred.data = data.frame(<br></div><div dir="ltr">Feature2<-d$Feature2,<br></div><div dir="ltr">Temp<-d$Temp<br></div><div dir="ltr">)<br></div><div dir="ltr">pred.data$Production = predict(mG2, newdata=pred.data, type="response")<br></div><div dir="ltr"><br></div><div dir="ltr">#Realizo o plot<br></div><div dir="ltr">library("ggplot2")<br></div><div dir="ltr">ggplot(d, aes(Temp, Production, colour = Feature2)) +<br></div><div dir="ltr"> geom_point() +<br></div><div dir="ltr"> geom_line(data=pred.data) +<br></div><div dir="ltr"> stat_smooth(method = "glm", formula = y ~ x, family = Gamma)<br></div><div dir="ltr">#<br></div><div dir="ltr"><br></div><div dir="ltr"> E tenho uma representação errada do intervalo de confiança do meu <br></div><div dir="ltr">modelo mG2 porque não consigo fazer a função stat_smooth() considerar o <br></div><div dir="ltr">ajuste mG2. Como posso resolver isso?<br></div><div dir="ltr"><br></div><div dir="ltr">Obrigado<br></div><div dir="ltr"><br></div><div dir="ltr">-- <br></div><div dir="ltr">======================================================================<br></div><div dir="ltr">Alexandre dos Santos<br></div><div dir="ltr">Proteção Florestal<br></div><div dir="ltr">IFMT - Instituto Federal de Educação, Ciência e Tecnologia de Mato Grosso<br></div><div dir="ltr">Campus Cáceres<br></div><div dir="ltr">Caixa Postal 244<br></div><div dir="ltr">Avenida dos Ramires, s/n<br></div><div dir="ltr">Bairro: Distrito Industrial<br></div><div dir="ltr">Cáceres - MT CEP: 78.200-000<br></div><div dir="ltr">Fone: (+55) 65 99686-6970 (VIVO) (+55) 65 3221-2674 (FIXO)<br></div><div dir="ltr">e-mails:<a ymailto="mailto:alexandresantosbr@yahoo.com.br" href="mailto:alexandresantosbr@yahoo.com.br">alexandresantosbr@yahoo.com.br</a><br></div><div dir="ltr"> <a ymailto="mailto:alexandre.santos@cas.ifmt.edu.br" href="mailto:alexandre.santos@cas.ifmt.edu.br">alexandre.santos@cas.ifmt.edu.br</a><br></div><div dir="ltr">Lattes: <a href="http://lattes.cnpq.br/1360403201088680" target="_blank">http://lattes.cnpq.br/1360403201088680</a><br></div><div dir="ltr">OrcID: orcid.org/0000-0001-8232-6722 - ResearcherID: A-5790-2016<br></div><div dir="ltr">Researchgate: www.researchgate.net/profile/Alexandre_Santos10<br></div><div dir="ltr">LinkedIn: br.linkedin.com/in/alexandre-dos-santos-87961635<br></div><div dir="ltr">Mendeley:www.mendeley.com/profiles/alexandre-dos-santos6/<br></div><div dir="ltr">======================================================================<br></div><div dir="ltr"><br></div><div dir="ltr"><br></div><div dir="ltr"><br></div><div dir="ltr">------------------------------<br></div><div dir="ltr"><br></div><div dir="ltr">Subject: Legenda do Digest<br></div><div dir="ltr"><br></div><div dir="ltr">_______________________________________________<br></div><div dir="ltr">R-br mailing list<br></div><div dir="ltr"><a ymailto="mailto:R-br@listas.c3sl.ufpr.br" href="mailto:R-br@listas.c3sl.ufpr.br">R-br@listas.c3sl.ufpr.br</a><br></div><div dir="ltr"><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><br></div><div dir="ltr"><br></div><div dir="ltr"><br></div><div dir="ltr">------------------------------<br></div><div dir="ltr"><br></div><div dir="ltr">Fim da Digest R-br, volume 96, assunto 8<br></div><div dir="ltr">****************************************<br></div></div>
</div>
</div></body></html>