Caros amigos, no exemplo abaixo, como faço para colocar as identificações do indivíduo em ordem crescente? att library(nlme) library(lattice) install.packages("tidyr") library(tidyr) DADOS.CURTO <- read.csv("https://www.dropbox.com/s/b4cckybgrcg86me/ galinhas.csv?raw=1 <https://n1.nylas.com/link/MESSAGE_ID/0?redirect=https%3A%2F%2Fwww.dropbox.com%2Fs%2Fb4cckybgrcg86me%2Fgalinhas.csv%3Fraw%3D1>",head = TRUE) head(DADOS.CURTO) colnames(DADOS.CURTO) <- c("0","1","2","3","4","5","6","TRAT") DADOS<- DADOS.CURTO %>% gather('0','1','2','3','4','5','6',key ="SEMANA",value = "PESO", -TRAT) head(DADOS) str(DADOS) ID <-rep(seq(1,370), times=7) DADOS<-cbind(ID,DADOS) DADOS <- transform(DADOS, SEMANA = as.numeric(SEMANA),ID=factor(ID)) 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 =DADOS) modelo.nlme.0<-nlme(modelo1.list) modelo.nlme.2<-update(modelo.nlme.0,fixed=list(a~TRAT,b + c~1),start=c(4.10,0,0,4.77,0.34),control=controle2) (corr2 <- plot(ACF(modelo.nlme.2,resType="pearson"),alpha = 0.05,ylab="autocorrelação",col = "gray",main = "B")) MM3.3<- update(modelo.nlme.2,weights = varExp(form = ~fitted(.)|SEMANA),correlation = corARMA(form = ~SEMANA |ID,q = 4)) #-------------------Gráfico---------------------------------------------------------------- s <- seq(1, nlevels(GALINHA$ID), by = 10) (RANDOMEFFECTS <- plot(ranef(MM3.3),id=0.05,idResType="normalized",layout=c(3,1),scales = list(y = list(tck = 1, at = s, labels = levels(GALINHA$ID)[s]),x=list(relation="free")),abline=c(v=0),ylab="Indivíduo",strip=TRUE)) -- ========================================= 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/ ==========================================