Você não recebeu nenhum erro na chamada a nlsList() não? Você tentou ver o resultado, com summary ou print?
2017-09-08 19:19 GMT-03:00 Fernando Souza via R-br <r-br@listas.c3sl.ufpr.br >: > Caros colegas ao rodar a função augPred {nlme} recebo a seguinte mensagem: > > "Error in tapply(as.character(object[[nm]]), groups, FUN[[dClass]]) : > arguments must have same length" > > Meus dados possuem diferentes números de individuos por tratamento e há > alguns NA's. Acredito que é por isso que nao funciona. > > Há uma forma de contornar essa questão > > #install.packages("tidyr") > #install.packages("nlme") > > library(tidyr) #pacote para organizar seu banco de dados > library(nlme) > DADOS <- read.csv("https://www.dropbox.com/s/b4cckybgrcg86me/ > galinhas.csv?raw=1",head = TRUE) > colnames(DADOS) <- c("0","1","2","3","4","5","6","TRAT") > DADOS <- DADOS %>% gather('0','1','2','3','4','5','6',key ="SEMANA",value > = "PESO", -TRAT) > ID <-rep(seq(1,370), times=7) > DADOS<-cbind(ID,DADOS) > GALINHA <- groupedData(PESO~SEMANA |ID,data = DADOS) > #----------------------------------------------------------- > -------------------------------------------- > install.packages("nlme") > library(nlme) > 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) > > modelo.nlme.0<-nlme(modelo1.list) > > plot(augPred(modelo.nlme.0)) > > o plot acima retorna a seguinte mensagem de erro, que é devido a função > augPred > > " Error in tapply(as.character(object[[nm]]), groups, FUN[[dClass]]) : > arguments must have same length" > > -- > ========================================= > Fernando Souza > Zootecnista, DSc. Produção e Alimentação Animal > Celular: (31)99796-8781 (Vivo) > E-mail:nandodeso...@gmail.com <e-mail%3anandodeso...@gmail.com> > Lattes: http://lattes.cnpq.br/6519538815038307 > Blog: https://producaoanimalcomr.wordpress.com/ > ========================================== > > _______________________________________________ > R-br mailing list > R-br@listas.c3sl.ufpr.br > https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br > Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça > código mínimo reproduzível. >
_______________________________________________ R-br mailing list R-br@listas.c3sl.ufpr.br https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forne�a c�digo m�nimo reproduz�vel.