<div dir="ltr"><div>> packages <- c("tidyr", "nlme""latticeExtra","lattice")</div><div>Erro: unexpected string constant in "packages <- c("tidyr", "nlme""latticeExtra""</div><div>> </div><div>> #Instala pacotes necessários, se ja tiver instalado ativa-o</div><div>> for(p in packages){</div><div>+ if(!require(p,character.only = TRUE)) install.packages(p)</div><div>+ library(p,character.only = TRUE)</div><div>+ }</div><div>Error in packages : objeto 'packages' não encontrado</div><div>> </div><div><br></div></div><div class="gmail_extra"><br><div class="gmail_quote">2018-07-13 17:48 GMT-03:00 Fernando Souza via R-br <span dir="ltr"><<a href="mailto:r-br@listas.c3sl.ufpr.br" target="_blank">r-br@listas.c3sl.ufpr.br</a>></span>:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div>Caros colegas</div><br><div>Desejo fazer um gráfico, semelhante ao retornado pela função augPred{nlme}, contendo os valores observados, os valores preditos pelo modelo fixo, e os valores individuais preditos. como esse exemplo (Só que para modelo não linear)</div><div><a href="https://rdrr.io/cran/nlme/man/plot.augPred.html" title="https://rdrr.io/cran/nlme/man/plot.augPred.html" target="_blank">https://rdrr.io/cran/nlme/man/<wbr>plot.augPred.html</a><br></div><div>por algum motivo a função augPred() não está funcionando com o meu modelo (Acredito por ser desbalanceado )</div><br><div>A função <a href="http://predict.nlme" title="predict.nlme" target="_blank">predict.nlme</a> {nlme}, me retorna um dataframe com as informação que quero, porém como transformá-los em gráfico?</div><br><div>Sei que é possível faze-lo através do xyplot{lattice} porém não estou conseguindo, devido a minha limitação no uso desse pacote.</div><br><div>Alguém poderia me ajudar?</div><br><div>segue CMR</div><br><br><br><br><div>packages <- c("tidyr", "nlme""latticeExtra","lattice"<wbr>)</div><br><div>#Instala pacotes necessários, se ja tiver instalado ativa-o</div><div>for(p in packages){</div><div>  if(!require(p,character.only = TRUE)) install.packages(p)</div><div>  library(p,character.only = TRUE)</div><div>}</div><br><div>##============================<wbr>==============================<wbr>==============================<wbr>============</div><div>##                              **Preparação do banco de dados**</div><div>##============================<wbr>==============================<wbr>==============================<wbr>============</div><div>DADOS.CURTO <- read.csv("<a href="https://www.dropbox.com/s/b4cckybgrcg86me/galinhas.csv?raw=1" target="_blank">https://www.dropbox.<wbr>com/s/b4cckybgrcg86me/<wbr>galinhas.csv?raw=1</a>",head = TRUE)</div><div>colnames(DADOS.CURTO) <- c("0","1","2","3","4","5","6",<wbr>"TRAT")</div><div>DADOS.CURTO <-DADOS.CURTO[with(DADOS.<wbr>CURTO,order(TRAT)),]</div><div>DADOS.CURTO$ID <- c(1:93,94:280,281:370)</div><div>DADOS<- DADOS.CURTO %>% gather('0','1','2','3','4','5'<wbr>,'6',key ="SEMANA",value = "PESO", -TRAT,-ID)</div><div>DADOS <- transform(DADOS, SEMANA = as.numeric(SEMANA),ID=factor(<wbr>ID))</div><div>DADOS</div><div>##----------------------------<wbr>------------------------------<wbr>----------</div><div>GALINHA <- groupedData(PESO~SEMANA |ID,order.groups=FALSE,data = DADOS)</div><div>##----------------------------<wbr>------------------------------<wbr>------------------------------<wbr>---</div><div>modelo1.list <- nlsList(PESO ~ a * exp(-b * exp( -c * SEMANA)),start = c(a = 4.10,b = 4.53,c = 0.37 ),na.action=na.omit,data =GALINHA)</div><div>##--------------Definição do modelo de efeito Fixo--------------------------<wbr>------</div><div>controle2<-nlmeControl(<wbr>maxIter=300,msMaxIter = 300,niterEM= 100)</div><div>modelo.nlme.0<-nlme(modelo1.<wbr>list)</div><div>modelo.nlme.2<-update(modelo.<wbr>nlme.0,fixed=list(a~TRAT,b + c~1),start=c(4.10,0,0,4.77,0.<wbr>34),control=controle2)</div><div>##----------------------------<wbr>------------------------------<wbr>------------------------------<wbr>--------------</div><div>controle3<-nlmeControl(<wbr>maxIter=300,msMaxIter=1000)</div><div>MM3.3<- update(modelo.nlme.2,weights = varExp(form = ~fitted(.)|SEMANA),control=<wbr>controle3,correlation = corARMA(form = ~SEMANA |ID,q = 4))</div><div>summary(MM3.3)</div><div>##----------------------------<wbr>------------------------------<wbr>-------</div><div>DADOS.ORIGINAL <- getData(MM3.3)</div><div>predito<-predict(MM3.3)</div><div>##----------------------------<wbr>------------------------------<wbr>--------</div><br><div><u></u><div><div>Sent from <a href="https://getmailspring.com/" target="_blank">Mailspring</a>, the best free email app for work</div></div><u></u></div><br>______________________________<wbr>_________________<br>
R-br mailing list<br>
<a href="mailto:R-br@listas.c3sl.ufpr.br">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/<wbr>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-<wbr>guia</a>) e forneça código mínimo reproduzível.<br></blockquote></div><br></div>