[R-br] Modelo Misto - plot (augPred()) mensagem de erro

[Fernando Souza]

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 em gmail.com <e-mail%3Anandodesouza em gmail.com>
Lattes: http://lattes.cnpq.br/6519538815038307
Blog: https://producaoanimalcomr.wordpress.com/
==========================================
-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20170908/f8776aef/attachment.html>