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",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))
--