Continuo "quebrando" a cabeça para tentar produzir o gráfico que preciso. Consegui chegar próximo mas ainda não é o ideal.
1) acrescentar a legenda para cada uma das linhas (estou denominando cada uma delas de m1,m2,m3 e m4);
2) substituir os números (1, 2 e 3) que aparecem nas divisões em amarelo por nomes (pode ser intercepto, tratamento B e tratamento C, respectivamente).
m1a<-c(94.1,91.9,92.1,92.8,92.3)
m2a<-c(96.3,95.7,95.8,95.1,95.6)
m3a<-c(91.6,91.0,91.1,90.0,88.6)
m4a<-c(95.9,92.9,100.0,90.9,91.1)
mma=c(m1a,m2a,m3a,m4a)
mat1=matrix(mma,nrow = 5, ncol = 4, byrow = FALSE,dimnames = list(c("L1","L2","L3","L4","L5"),c("m1","m2","m3","m4")))
m1b<-c(14.1,11.9,16.1,12.8,9.3)
m2b<-c(16.3,15.7,19.8,19.1,9.6)
m3b<-c(11.6,11.0,19.1,9.0,8.6)
m4b<-c(15.9,12.9,10.0,19.9,11.1)
mmb=c(m1b,m2b,m3b,m4b)
mat2=matrix(mmb,nrow = 5, ncol = 4, byrow = FALSE,dimnames = list(c("L1","L2","L3","L4","L5"),c("m1","m2","m3","m4")))
m1c<-c(34.1,51.9,42.1,32.8,62.3)
m2c<-c(36.3,45.7,35.8,55.1,75.6)
m3c<-c(51.6,71.0,71.1,70.0,68.6)
m4c<-c(45.9,42.9,40.0,30.9,51.1)
mmc=c(m1c,m2c,m3c,m4c)
mat3=matrix(mmc,nrow = 5, ncol = 4, byrow = FALSE,dimnames = list(c("L1","L2","L3","L4","L5"),c("m1","m2","m3","m4")))
matriz=matrix(c(mat1,mat2,mat3),5,12)
a=ts(matriz,frequency = 0.5, start = c(2.5, 2))
require(lattice)
xyplot(a, screens = c(1,1,1,1,2,2,2,2,3,3,3,3),xlab = "medidas de tempo",ylab = "%",col=c(1:4))