<html>
<head>
<meta http-equiv="Content-Type" content="text/html; charset=iso-8859-1">
<style type="text/css" style="display:none;"> P {margin-top:0;margin-bottom:0;} </style>
</head>
<body dir="ltr">
<div style="color: rgb(0, 0, 0); font-family: Calibri,Helvetica,sans-serif; font-size: 12pt;">
Boa tarde,</div>
<div style="color: rgb(0, 0, 0); font-family: Calibri,Helvetica,sans-serif; font-size: 12pt;">
Estou tentando modelar as funções gompertz e logístico porem quero somar as duas equações e gerar um único gráfico. Porem a equação terá 6 parâmetros teria alguma rotina para esse caso. Irei anexar a rotina utilizada quando uma equação é modelada de cada vez
e os dados MFF que eu utilizei para modelar essas funções. Meu novo arquivo de dados GG (y são dois caracteres juntos) por isso que na rotina teria que somar a equação de gompertz (y~b0*exp(-exp(b1-b2*x) porem quando somo essa equação com a mesma equação de
gompertz da erro. Eu terei 6 parâmetros nesse caso , já tentei indicar de varias formas o b0, b1 e b2 pois além desses tenho mais um b0, b1 e b2. Espero que eu tenha me feito entender.</div>
<div style="color: rgb(0, 0, 0); font-family: Calibri,Helvetica,sans-serif; font-size: 12pt;">
Att Cláudia de Bem </div>
<div>
<div id="appendonsend"></div>
<div style="color:rgb(0,0,0); font-family:Calibri,Helvetica,sans-serif; font-size:12pt">
<br>
</div>
<hr tabindex="-1" style="display:inline-block; width:98%">
<div id="divRplyFwdMsg" dir="ltr"><font color="#000000" face="Calibri, sans-serif" style="font-size:11pt"><b>De:</b> R-br <r-br-bounces@listas.c3sl.ufpr.br> em nome de ASANTOS por (R-br) <r-br@listas.c3sl.ufpr.br><br>
<b>Enviado:</b> terça-feira, 29 de janeiro de 2019 18:45<br>
<b>Para:</b> r-br@listas.c3sl.ufpr.br; r-br-request@listas.c3sl.ufpr.br<br>
<b>Cc:</b> ASANTOS<br>
<b>Assunto:</b> [R-br] Ajuda para plotar IC para modelo glm com junção de níveis</font>
<div> </div>
</div>
<div class="BodyFragment"><font size="2"><span style="font-size:11pt">
<div class="PlainText">Prezados Membros do R-br,<br>
<br>
Estou tentando plotar um modelo glmcom stat_smooth() no ggplot2 sem <br>
sucesso. Minha ideia é representar o intervalo de confiança para um <br>
modelo onde realizei a junção de dois níveis para o fator Feature. Ou <br>
seja, modelos iguais são representados por uma única curva, porém a <br>
função stat_smooth() não "enxerga" os níveis que foram juntados. Segundo <br>
exemplo fictício do meu problema:<br>
<br>
#Simulação de um banco de dados, onde a variável resposta Production tem <br>
distribuição Gamma<br>
set.seed(123)<br>
d<-NULL<br>
N<-50<br>
d$Production <- rgamma(N,10)<br>
d$Feature <- ifelse(d$Production >7 & d$Production<10, <br>
c("green"),ifelse(d$Production>11,<br>
c("red"), c("blue")))<br>
d$Temp<-rnorm(N,20,5)<br>
d<-as.data.frame(d)<br>
#<br>
<br>
# Ajusto um modelo de Gamma completo<br>
mG<- glm(Production~ Feature + Temp, family= Gamma, data = d)<br>
summary(mG)<br>
anova(mG,test="Chi")<br>
<br>
# Comparo os níveis do fator Feature<br>
library(multcomp)<br>
PW<-summary(glht(mG, linfct = mcp(Feature= "Tukey")))<br>
PW<br>
cld(PW)<br>
<br>
# Faço de conta que os níveis green = blue e realizo um novo ajuste<br>
Feature2<-d$Feature<br>
levels(Feature2)<br>
levels(Feature2)[1]<-"blue&green"<br>
levels(Feature2)[2]<-"blue&green"<br>
levels(Feature2)<br>
d$Feature2<-Feature2<br>
mG2<- glm(Production~ Feature2 + Temp, family= Gamma, data = d)<br>
<br>
# Crio os valores de predição<br>
pred.data = data.frame(<br>
Feature2<-d$Feature2,<br>
Temp<-d$Temp<br>
)<br>
pred.data$Production = predict(mG2, newdata=pred.data, type="response")<br>
<br>
#Realizo o plot<br>
library("ggplot2")<br>
ggplot(d, aes(Temp, Production, colour = Feature2)) +<br>
geom_point() +<br>
geom_line(data=pred.data) +<br>
stat_smooth(method = "glm", formula = y ~ x, family = Gamma)<br>
#<br>
<br>
E tenho uma representação errada do intervalo de confiança do meu <br>
modelo mG2 porque não consigo fazer a função stat_smooth() considerar o <br>
ajuste mG2. Como posso resolver isso?<br>
<br>
Obrigado<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>
e-mails:alexandresantosbr@yahoo.com.br<br>
alexandre.santos@cas.ifmt.edu.br<br>
Lattes: <a href="http://lattes.cnpq.br/1360403201088680">http://lattes.cnpq.br/1360403201088680</a><br>
OrcID: orcid.org/0000-0001-8232-6722 - ResearcherID: A-5790-2016<br>
Researchgate: <a href="http://www.researchgate.net/profile/Alexandre_Santos10">www.researchgate.net/profile/Alexandre_Santos10</a><br>
LinkedIn: br.linkedin.com/in/alexandre-dos-santos-87961635<br>
Mendeley:www.mendeley.com/profiles/alexandre-dos-santos6/<br>
======================================================================<br>
<br>
_______________________________________________<br>
R-br mailing list<br>
R-br@listas.c3sl.ufpr.br<br>
<a href="https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br">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">http://www.leg.ufpr.br/r-br-guia</a>) e fornea cdigo mnimo reproduzvel.</div>
</span></font></div>
</div>
</body>
</html>