<div dir="ltr"><div dir="ltr">Edmar,<div><br></div><div>O CMR precisa atender a cláusula <b>reproduzível </b>se você espera que a gente possa entender seu problema:<br><div><br></div><div>R version 3.5.1 (2018-07-02) -- "Feather Spray"</div><div>Copyright (C) 2018 The R Foundation for Statistical Computing</div><div>Platform: i386-w64-mingw32/i386 (32-bit)</div><div><br></div><div>R é um software livre e vem sem GARANTIA ALGUMA.</div><div>Você pode redistribuí-lo sob certas circunstâncias.</div><div>Digite 'license()' ou 'licence()' para detalhes de distribuição.</div><div><br></div><div>R é um projeto colaborativo com muitos contribuidores.</div><div>Digite 'contributors()' para obter mais informações e</div><div>'citation()' para saber como citar o R ou pacotes do R em publicações.</div><div><br></div><div>Digite 'demo()' para demonstrações, 'help()' para o sistema on-line de ajuda,</div><div>ou 'help.start()' para abrir o sistema de ajuda em HTML no seu navegador.</div><div>Digite 'q()' para sair do R.</div><div><br></div><div>> attach(waste)<br></div><div>Error in attach(waste) : objeto 'waste' não encontrado</div><div>></div><br></div><div>Concorda que não dá para seguir daqui?</div><div><br></div><div><br></div></div></div><br><div class="gmail_quote"><div dir="ltr">On Mon, Dec 10, 2018 at 12:08 PM Edmar Caldas por (R-br) <<a href="mailto:r-br@listas.c3sl.ufpr.br">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"><div><div class="gmail-m_-6698138850097310284ydp5556b043yahoo-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="gmail-m_-6698138850097310284ydp5556b043signature"><div><font color="#0000ff">Edmar</font></div></div></div>
        <div>nsegue </div><div><br></div>
        
        </div><div id="gmail-m_-6698138850097310284yahoo_quoted_4860631892" class="gmail-m_-6698138850097310284yahoo_quoted">
            <div style="font-family:"Helvetica Neue",Helvetica,Arial,sans-serif;font-size:13px;color:rgb(38,40,42)">
                
                <div>
                    Em segunda-feira, 10 de dezembro de 2018 12:00:18 BRST,  <<a href="mailto:r-br-request@listas.c3sl.ufpr.br" target="_blank">r-br-request@listas.c3sl.ufpr.br</a>> 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 href="mailto:r-br@listas.c3sl.ufpr.br" target="_blank">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 href="mailto:r-br-request@listas.c3sl.ufpr.br" target="_blank">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 href="mailto:r-br-owner@listas.c3sl.ufpr.br" target="_blank">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 href="mailto:eliaskrainski@yahoo.com.br" target="_blank">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 href="mailto:r-br@listas.c3sl.ufpr.br" target="_blank">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 href="mailto:815218723.1276327.1544374485235@mail.yahoo.com" target="_blank">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 href="mailto:r-br@listas.c3sl.ufpr.br" target="_blank">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 href="mailto:R-br@listas.c3sl.ufpr.br" target="_blank">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 href="mailto:alexandre.santos@cas.ifmt.edu.br" target="_blank">alexandre.santos@cas.ifmt.edu.br</a>><br></div><div dir="ltr">To: "<a href="mailto:r-br@listas.c3sl.ufpr.br" target="_blank">r-br@listas.c3sl.ufpr.br</a>" <<a href="mailto:r-br@listas.c3sl.ufpr.br" target="_blank">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 href="mailto:8d326872-10b9-06f7-cc1b-564caf54f3e7@cas.ifmt.edu.br" target="_blank">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 href="mailto:alexandresantosbr@yahoo.com.br" target="_blank">alexandresantosbr@yahoo.com.br</a><br></div><div dir="ltr">         <a href="mailto:alexandre.santos@cas.ifmt.edu.br" target="_blank">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: <a href="http://orcid.org/0000-0001-8232-6722" target="_blank">orcid.org/0000-0001-8232-6722</a>   -   ResearcherID: A-5790-2016<br></div><div dir="ltr">Researchgate: <a href="http://www.researchgate.net/profile/Alexandre_Santos10" target="_blank">www.researchgate.net/profile/Alexandre_Santos10</a><br></div><div dir="ltr">LinkedIn: <a href="http://br.linkedin.com/in/alexandre-dos-santos-87961635" target="_blank">br.linkedin.com/in/alexandre-dos-santos-87961635</a><br></div><div dir="ltr">Mendeley:<a href="http://www.mendeley.com/profiles/alexandre-dos-santos6/" target="_blank">www.mendeley.com/profiles/alexandre-dos-santos6/</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">------------------------------<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 href="mailto:R-br@listas.c3sl.ufpr.br" target="_blank">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></div>_______________________________________________<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.</blockquote></div>