Segue CMR. require(lattice) df <- data.frame(id=c(1,1,2,2,1,1,3,3,4,4), grp=c(1,1,1,1,2,2,2,2,2,2), x=1:10, y=1:10, z=jitter(1:10)) df <- within(df, { id <- factor(id) grp <- factor(grp) x <- factor(x) y <- factor(y) }) with(df, nlevels(id)*nlevels(grp)*nlevels(x)) with(df, nlevels(id)*nlevels(grp)) with(df, nlevels(x)) ## Opção 1: colocar um layout que seja multiplo. xyplot(z~y|id*grp*x, data=df, layout=c(8, 10), type="b", main="Test", ylab="y axis", xlab="x axis") ## Opção 2: criar novos fatores pela combinação dos envolvidos. df <- transform(df, id.grp=interaction(id, grp), id.grp.x=interaction(id, grp, x, drop=FALSE)) xyplot(z~y|id.grp*x, data=df, type="b", main="Test", ylab="y axis", xlab="x axis") xyplot(z~y|id.grp.x, data=df, type="b", main="Test", ylab="y axis", xlab="x axis") ## Visto que nem todas as combinações existem, o último gráfico é mais ## apropriado porque salva espaço. À disposição. Walmes.