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

[Cesar Rabak]

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