
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:nandodesouza@gmail.com <e-mail%3Anandodesouza@gmail.com> Lattes: http://lattes.cnpq.br/6519538815038307 Blog: https://producaoanimalcomr.wordpress.com/ ==========================================