[R-br] Ordernar eixo y (indivíduos) em gráfico xyplot Lattice

Fernando Souza nandodesouza em gmail.com
Terça Janeiro 23 13:40:09 -02 2018


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 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/20180123/d7bee79a/attachment.html>


Mais detalhes sobre a lista de discussão R-br