<html xmlns:v="urn:schemas-microsoft-com:vml" xmlns:o="urn:schemas-microsoft-com:office:office" xmlns:w="urn:schemas-microsoft-com:office:word" xmlns:x="urn:schemas-microsoft-com:office:excel" xmlns:m="http://schemas.microsoft.com/office/2004/12/omml" xmlns="http://www.w3.org/TR/REC-html40"><head><meta http-equiv=Content-Type content="text/html; charset=iso-8859-1"><meta name=Generator content="Microsoft Word 14 (filtered medium)"><!--[if !mso]><style>v\:* {behavior:url(#default#VML);}
o\:* {behavior:url(#default#VML);}
w\:* {behavior:url(#default#VML);}
.shape {behavior:url(#default#VML);}
</style><![endif]--><style><!--
/* Font Definitions */
@font-face
{font-family:"Cambria Math";
panose-1:2 4 5 3 5 4 6 3 2 4;}
@font-face
{font-family:Cambria;
panose-1:2 4 5 3 5 4 6 3 2 4;}
@font-face
{font-family:Calibri;
panose-1:2 15 5 2 2 2 4 3 2 4;}
@font-face
{font-family:Tahoma;
panose-1:2 11 6 4 3 5 4 4 2 4;}
@font-face
{font-family:"Trebuchet MS";
panose-1:2 11 6 3 2 2 2 2 2 4;}
/* Style Definitions */
p.MsoNormal, li.MsoNormal, div.MsoNormal
{margin:0cm;
margin-bottom:.0001pt;
font-size:12.0pt;
font-family:"Times New Roman","serif";}
a:link, span.MsoHyperlink
{mso-style-priority:99;
color:blue;
text-decoration:underline;}
a:visited, span.MsoHyperlinkFollowed
{mso-style-priority:99;
color:purple;
text-decoration:underline;}
p.MsoAcetate, li.MsoAcetate, div.MsoAcetate
{mso-style-priority:99;
mso-style-link:"Texto de balão Char";
margin:0cm;
margin-bottom:.0001pt;
font-size:8.0pt;
font-family:"Tahoma","sans-serif";}
span.EstiloDeEmail17
{mso-style-type:personal-reply;
font-family:"Cambria","serif";
color:#C00000;}
span.TextodebaloChar
{mso-style-name:"Texto de balão Char";
mso-style-priority:99;
mso-style-link:"Texto de balão";
font-family:"Tahoma","sans-serif";
mso-fareast-language:PT-BR;}
.MsoChpDefault
{mso-style-type:export-only;
font-family:"Calibri","sans-serif";
mso-fareast-language:EN-US;}
@page WordSection1
{size:612.0pt 792.0pt;
margin:70.85pt 3.0cm 70.85pt 3.0cm;}
div.WordSection1
{page:WordSection1;}
--></style><!--[if gte mso 9]><xml>
<o:shapedefaults v:ext="edit" spidmax="1026" />
</xml><![endif]--><!--[if gte mso 9]><xml>
<o:shapelayout v:ext="edit">
<o:idmap v:ext="edit" data="1" />
</o:shapelayout></xml><![endif]--></head><body lang=PT-BR link=blue vlink=purple><div class=WordSection1><p class=MsoNormal><span style='font-family:"Cambria","serif";color:#C00000'>Caro Walmes, boa tarde.<o:p></o:p></span></p><p class=MsoNormal><span style='font-family:"Cambria","serif";color:#C00000'><o:p> </o:p></span></p><p class=MsoNormal><span style='font-family:"Cambria","serif";color:#C00000'>Li o comentário sobre “curva exponencial” no e-mail abaixo endereçado ao grupo.<o:p></o:p></span></p><p class=MsoNormal><span style='font-family:"Cambria","serif";color:#C00000'>Pois bem, estou com um probleminha para resolver, a saber:<o:p></o:p></span></p><p class=MsoNormal><span style='font-family:"Cambria","serif";color:#C00000'><o:p> </o:p></span></p><p class=MsoNormal><span style='font-family:"Cambria","serif";color:#C00000'>Tenho estas duas séries:<o:p></o:p></span></p><p class=MsoNormal><span style='font-family:"Cambria","serif";color:#C00000'><o:p> </o:p></span></p><p class=MsoNormal><span style='font-family:"Cambria","serif";color:#C00000'>x<-c(1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16,17,18,19,20,21,22,23)<o:p></o:p></span></p><p class=MsoNormal><span style='font-family:"Cambria","serif";color:#C00000'>y<-c(100,150,200,250,300,350,400,450,500,550,600,650,700,750,800,850,900,950,1000,1050,1100,1150,1200)<o:p></o:p></span></p><p class=MsoNormal><span style='font-family:"Cambria","serif";color:#C00000'><o:p> </o:p></span></p><p class=MsoNormal><span style='font-family:"Cambria","serif";color:#C00000'>Utilizei um script apresentado no Livro Conhecendo o R de Luiz A Pertenelli e Márcio P de Mello pág. 165 e não consegui resolver no R.<o:p></o:p></span></p><p class=MsoNormal><span style='font-family:"Cambria","serif";color:#C00000'><o:p> </o:p></span></p><p class=MsoNormal><span style='font-family:"Cambria","serif";color:#C00000'>Bom! Ajustei a curva y para exponencial Y no EXCEL conforme abaixo:<o:p></o:p></span></p><p class=MsoNormal><span style='font-family:"Cambria","serif";color:#C00000'><o:p> </o:p></span></p><p class=MsoNormal><span style='font-family:"Cambria","serif";color:#C00000'>Y<-c(187,206,227,250,275,303,334,368,405,446,491,541,596,657,723,797,877,966,1064,1172,1291,1422,1566)<o:p></o:p></span></p><p class=MsoNormal><span style='font-size:10.0pt;font-family:"Arial","sans-serif";color:black'><o:p> </o:p></span></p><p class=MsoNormal><span style='font-size:10.0pt;font-family:"Arial","sans-serif";color:black'><o:p> </o:p></span></p><p class=MsoNormal><span style='font-family:"Cambria","serif";color:#C00000'>Por favor, como ajusto à exponencial a série y no R? Como chegar ao resultado Y no R?<o:p></o:p></span></p><p class=MsoNormal><span style='font-family:"Cambria","serif";color:#C00000'><o:p> </o:p></span></p><p class=MsoNormal><span style='font-family:"Cambria","serif";color:#C00000'>#Script utilizado que não funcionou:<o:p></o:p></span></p><p class=MsoNormal><span style='font-family:"Cambria","serif";color:#C00000'>x<-c(1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16,17,18,19,20,21,22,23)<o:p></o:p></span></p><p class=MsoNormal><span style='font-family:"Cambria","serif";color:#C00000'>y<-c(100,150,200,250,300,350,400,450,500,550,600,650,700,750,800,850,900,950,1000,1050,1100,1150,1200)<o:p></o:p></span></p><p class=MsoNormal><span style='font-family:"Cambria","serif";color:#C00000'>dados<-data.frame(x,y)<o:p></o:p></span></p><p class=MsoNormal><span style='font-family:"Cambria","serif";color:#C00000'>plot(dados)<o:p></o:p></span></p><p class=MsoNormal><span style='font-family:"Cambria","serif";color:#C00000'>funcao<-y~a*exp(b*x)<o:p></o:p></span></p><p class=MsoNormal><span style='font-family:"Cambria","serif";color:#C00000'>exponencial<-nls(funcao, dados, start=c(a=1,b=1))<o:p></o:p></span></p><p class=MsoNormal><span style='font-family:"Cambria","serif";color:#C00000'>exponencial<o:p></o:p></span></p><p class=MsoNormal><span style='font-family:"Cambria","serif";color:#C00000'><o:p> </o:p></span></p><p class=MsoNormal><span style='font-family:"Cambria","serif";color:#C00000'>=====Erro no R ==================<o:p></o:p></span></p><p class=MsoNormal><span style='font-family:"Cambria","serif";color:#C00000'>10,11,12,13,14,15,16,17,18,19,20,21,22,23)<o:p></o:p></span></p><p class=MsoNormal><span style='font-family:"Cambria","serif";color:#C00000'>> y<-c(100,150,200,250,300,350,400,450,500,550,600,650,700,750,800,850,900,950,1000,1050,1100,1150,1200)<o:p></o:p></span></p><p class=MsoNormal><span style='font-family:"Cambria","serif";color:#C00000'>> dados<-data.frame(x,y)<o:p></o:p></span></p><p class=MsoNormal><span style='font-family:"Cambria","serif";color:#C00000'>> plot(dados)<o:p></o:p></span></p><p class=MsoNormal><span style='font-family:"Cambria","serif";color:#C00000'>> funcao<-y~a*exp(b*x)<o:p></o:p></span></p><p class=MsoNormal><span style='font-family:"Cambria","serif";color:#C00000'>> exponencial<-nls(funcao, dados, start=c(a=1,b=1))<o:p></o:p></span></p><p class=MsoNormal><b><span style='font-family:"Cambria","serif"'>Erro em nls(funcao, dados, start = c(a = 1, b = 1)) : gradiente singular<o:p></o:p></span></b></p><p class=MsoNormal><span style='font-family:"Cambria","serif";color:#C00000'>> exponencial<o:p></o:p></span></p><p class=MsoNormal><b><span style='font-family:"Cambria","serif"'>Erro: objeto 'exponencial' não encontrado<o:p></o:p></span></b></p><p class=MsoNormal><span style='font-family:"Cambria","serif";color:#C00000'>><o:p> </o:p></span></p><p class=MsoNormal><span style='font-family:"Cambria","serif";color:#C00000'><o:p> </o:p></span></p><p class=MsoNormal><span style='font-family:"Cambria","serif";color:#C00000'>Onde estou errando?<o:p></o:p></span></p><p class=MsoNormal><span style='font-family:"Cambria","serif";color:#C00000'><o:p> </o:p></span></p><p class=MsoNormal><span style='font-family:"Cambria","serif";color:#C00000'><o:p> </o:p></span></p><p class=MsoNormal><span style='font-size:11.0pt;font-family:"Cambria","serif";color:#C00000'>Abs.<o:p></o:p></span></p><div class=MsoNormal align=center style='text-align:center'><span style='font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1F497D'><hr size=2 width="100%" align=center></span></div><p class=MsoNormal><a href="http://www.posicional.com/"><span style='font-size:11.0pt;font-family:"Calibri","sans-serif";color:#C00000;text-decoration:none'><img border=0 width=409 height=94 id="Imagem_x0020_2" src="cid:image001.gif@01CEEC52.A5DC3A70" alt="assinatura_email_sites"></span></a><span style='font-size:11.0pt;font-family:"Calibri","sans-serif";color:#C00000'><o:p></o:p></span></p><p class=MsoNormal><span style='font-size:11.0pt;font-family:"Calibri","sans-serif";color:#C00000'><o:p> </o:p></span></p><p class=MsoNormal><b><span style='font-size:14.0pt;font-family:"Cambria","serif";color:black'>Cícero C Nunes<o:p></o:p></span></b></p><p class=MsoNormal><span style='font-size:8.0pt;font-family:"Cambria","serif";color:#C00000'>Diretor Técnico Operacional<o:p></o:p></span></p><p class=MsoNormal><span style='font-size:8.0pt;font-family:"Cambria","serif";color:#C00000'>Divisão Remuneração<o:p></o:p></span></p><p class=MsoNormal><span style='font-size:10.0pt;font-family:"Cambria","serif";color:#C00000'><a href="mailto:c2n@posicional.com"><span style='color:blue'>c2n@posicional.com</span></a><o:p></o:p></span></p><p class=MsoNormal><span style='font-size:11.0pt;font-family:"Cambria","serif";color:#1F497D'><o:p> </o:p></span></p><p class=MsoNormal><b><span style='font-size:9.0pt;font-family:"Cambria","serif";color:#C00000'>+55(16)3397.0226<o:p></o:p></span></b></p><p class=MsoNormal><b><span style='font-size:9.0pt;font-family:"Cambria","serif";color:#C00000'>+55(11)3280.0226<o:p></o:p></span></b></p><p class=MsoNormal><b><span style='font-size:9.0pt;font-family:"Cambria","serif";color:#C00000'>+55(11)9.9978.0026<o:p></o:p></span></b></p><p class=MsoNormal><span style='font-size:11.0pt;font-family:"Cambria","serif";color:#1F497D'><o:p> </o:p></span></p><p class=MsoNormal><b><span style='font-size:9.0pt;font-family:"Cambria","serif";color:#1F497D'>SKYPE</span></b><span style='font-size:9.0pt;font-family:"Cambria","serif";color:#1F497D'>: </span><span style='font-size:9.0pt;font-family:"Cambria","serif";color:black'>consultormp</span><span style='font-size:9.0pt;font-family:"Cambria","serif";color:#1F497D'><o:p></o:p></span></p><p class=MsoNormal><span style='font-size:9.0pt;font-family:"Cambria","serif";color:#595959'><o:p> </o:p></span></p><p class=MsoNormal><span style='font-size:9.0pt;font-family:"Cambria","serif";color:#595959'>Rua Padre Duarte, 1360 – salas 6/7 – Centro – Araraquara, SP – CEP.: 14.801-310<o:p></o:p></span></p><p class=MsoNormal><b><i><span style='font-size:9.0pt;font-family:"Cambria","serif";color:#C00000'>“Antes de imprimir pense em sua responsabilidade com o meio ambiente – Quem ama cuida!”<o:p></o:p></span></i></b></p><p class=MsoNormal><span style='font-family:"Cambria","serif";color:#C00000'><o:p> </o:p></span></p><p class=MsoNormal><b><span style='font-size:10.0pt;font-family:"Tahoma","sans-serif"'>De:</span></b><span style='font-size:10.0pt;font-family:"Tahoma","sans-serif"'> R-br [mailto:r-br-bounces@listas.c3sl.ufpr.br] <b>Em nome de </b>walmes .<br><b>Enviada em:</b> sexta-feira, 8 de novembro de 2013 18:33<br><b>Para:</b> r-br@listas.c3sl.ufpr.br<br><b>Assunto:</b> Re: [R-br] Contrastes em modelos de regressão de poisson inflacionados por zeros<o:p></o:p></span></p><p class=MsoNormal><o:p> </o:p></p><div><div><p class=MsoNormal><span style='font-family:"Trebuchet MS","sans-serif"'>Você não deve esperar ver uma curva exponencial para a média das contagens em função do tempo e uma sigmoide para a proporção de zeros não Poisson para os zeros se não simulou os dados com efeito de tempo e não declarou efeito de tempo na porção binomial. Para que isso aconteça a simulação tem que estar condizente, tem que ter os efeitos para que as curvas não sejam as "curvas nulas".<br><br></span><span style='font-family:"Courier New"'>#------------------------------------------------------------------<br># Definições da sessão.<br><br>rm(list=ls())<br>require(pscl)<br>require(VGAM)<br>require(multcomp)<br>require(lattice)<br>require(latticeExtra)<br><br>#------------------------------------------------------------------<br># Dados artificiais.<br><br>da <- expand.grid(tempo=rep(1:10), trat=gl(3,10))<br><br>## simula porção Poisson<br>X <- model.matrix(~trat+tempo, da)<br>colnames(X)<br>beta.pois <- c(-1,0.25,0.5,0.25)<br>eta.pois <- X%*%beta.pois<br>da$y.pois <- rpois(nrow(X), lambda=exp(eta.pois))<br>xyplot(y.pois~tempo|trat, da)<br><br>## simula porção binomial<br>X <- model.matrix(~tempo, da)<br>beta.bern <- c(0,0.3)<br>eta.bern <- X%*%beta.bern<br>da$y.bern <- rbinom(nrow(X), size=1, prob=1/(1+exp(-eta.bern)))<br><br>xyplot(y.bern~tempo|trat, da)<br><br>## monta a contagem com inflação de zeros<br>da$y <- with(da, ifelse(y.bern==0, 0, y.pois))<br>xyplot(y~tempo|trat, da)<br><br>#------------------------------------------------------------------<br># Modelo completo.<br><br>compl.mod <- zeroinfl(y~trat+tempo|trat+tempo, data=da)<br>coef(compl.mod)<br><br>#------------------------------------------------------------------<br># Predição do modelo considerando as duas porções.<br><br>pred <- expand.grid(tempo=rep(1:10), trat=gl(3,1))<br>X <- model.matrix(~trat+tempo, pred)<br>i <- grep("^count\\_", names(coef(compl.mod)))<br>eta <- X%*%coef(compl.mod)[i]<br>pred$m.pois <- exp(eta)<br>i <- grep("^zero\\_", names(coef(compl.mod)))<br>eta <- X%*%coef(compl.mod)[i]<br>pred$p.zero <- exp(eta)/(1+exp(eta))<br><br>xyplot(y~tempo|trat, data=da, jitter.x=TRUE)+<br> as.layer(xyplot(m.pois~tempo|trat, data=pred, type="l"))+<br> as.layer(xyplot(p.zero~tempo|trat, data=pred,<br> type="l", lty=2, lwd=2))+<br> layer(panel.abline(h=1, lty=2))<br><br>xyplot(p.zero~tempo|trat, data=pred,<br> type="l", lty=2, lwd=2)<br><br>#------------------------------------------------------------------<br></span><span style='font-family:"Trebuchet MS","sans-serif"'><br><br>À disposição.<br>Walmes.<o:p></o:p></span></p></div><div><p class=MsoNormal><br clear=all><o:p></o:p></p><div><div><p class=MsoNormal><span style='font-family:"Trebuchet MS","sans-serif"'>==========================================================================<br>Walmes Marques Zeviani<br>LEG (Laboratório de Estatística e Geoinformação, 25.450418 S, 49.231759 W)<br>Departamento de Estatística - Universidade Federal do Paraná<br>fone: (+55) 41 3361 3573<br>skype: walmeszeviani<br>homepage: <a href="http://www.leg.ufpr.br/%7Ewalmes" target="_blank">http://www.leg.ufpr.br/~walmes</a><br>linux user number: 531218<br>==========================================================================</span><o:p></o:p></p></div></div><p class=MsoNormal><o:p> </o:p></p></div></div></div></body></html>