Re: [R-br] Ajuda para elaboração de boxplot em painel

Prezad@s Estou com dificuldade de elaborar um boxplot por grupo (ou painel), no formato do que está no arquivo pdf em anexo a esse E-mail. Segue o CMR do banco de dados (Documentação do SAS) O arquivo pdf pode ser acessado em https://www.dropbox.com/s/5q65geuq7e9m5ey/Dados%26grafico-DocSAS.pdf?dl=0 Desde já agradeço Parts=structure(list(Machine = structure(c(1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L), .Label = c("A386", "A455", "C334"), class = "factor"), Sample = c(1L, 1L, 1L, 1L, 2L, 2L, 2L, 2L, 3L, 3L, 3L, 3L, 4L, 4L, 4L, 4L, 5L, 5L, 5L, 5L, 6L, 6L, 6L, 6L, 7L, 7L, 7L, 7L, 8L, 8L, 8L, 8L, 9L, 9L, 9L, 9L, 10L, 10L, 10L, 10L, 11L, 11L, 11L, 11L, 12L, 12L, 12L, 12L, 13L, 13L, 13L, 13L, 14L, 14L, 14L, 14L, 15L, 15L, 15L, 15L, 16L, 16L, 16L, 16L, 17L, 17L, 17L, 17L, 18L, 18L, 18L, 18L, 19L, 19L, 19L, 19L, 20L, 20L, 20L, 20L, 21L, 21L, 21L, 21L, 22L, 22L, 22L, 22L, 23L, 23L, 23L, 23L, 24L, 24L, 24L, 24L), Diam = c(4.32, 4.55, 4.16, 4.44, 4.49, 4.3, 4.52, 4.61, 4.44, 4.32, 4.25, 4.5, 4.55, 4.15, 4.42, 4.49, 4.21, 4.3, 4.29, 4.63, 4.56, 4.61, 4.29, 4.56, 4.63, 4.3, 4.41, 4.58, 4.38, 4.65, 4.43, 4.44, 4.12, 4.49, 4.3, 4.36, 4.45, 4.56, 4.38, 4.51, 4.62, 4.67, 4.7, 4.58, 4.33, 4.23, 4.34, 4.58, 4.29, 4.38, 4.28, 4.41, 4.15, 4.35, 4.28, 4.23, 4.21, 4.3, 4.32, 4.38, 4.16, 4.28, 4.31, 4.59, 4.14, 4.18, 4.08, 4.21, 4.51, 4.2, 4.28, 4.19, 4.1, 4.33, 4.37, 4.47, 3.99, 4.09, 4.47, 4.25, 4.24, 4.54, 4.43, 4.38, 4.23, 4.48, 4.31, 4.57, 4.27, 4.4, 4.32, 4.56, 4.7, 4.65, 4.49, 4.38)), .Names = c("Machine", "Sample", "Diam"), class = "data.frame", row.names = c(NA, -96L)) Gilênio Borges Fernandes Professor Associado IV (Aposentado) Professor Adjunto A (Substituto) Universidade Federal da Bahia Instituto de Matemática e Estatística Departamento de Estatística Av. Adhemar de Barros, s/n – Ondina. 40.170-110 - Salvador - BA, Brasil Tel.: (071)3283-6340/6341/6337 Fax: (071)3283-6336 Skype: gilenio.fernandes Lattes: http://lattes.cnpq.br/6764860618464860 Em 3 de maio de 2017 15:58, Gilenio Borges Fernandes <gilenio@ufba.br> escreveu:
Prezad@s
Estou com dificuldade de elaborar um boxplot por grupo (ou painel), no formato do que está no arquivo pdf em anexo a esse E-mail. Segue o CMR do banco de dados (Documentação do SAS)
Desde já agradeço
Parts=structure(list(Machine = structure(c(1L, 1L, 1L, 1L, 1L, 1L,
1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L,
1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 2L, 2L,
2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L,
2L, 2L, 2L, 2L, 2L, 2L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L,
3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L,
3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L), .Label = c("A386", "A455",
"C334"), class = "factor"), Sample = c(1L, 1L, 1L, 1L, 2L, 2L,
2L, 2L, 3L, 3L, 3L, 3L, 4L, 4L, 4L, 4L, 5L, 5L, 5L, 5L, 6L, 6L,
6L, 6L, 7L, 7L, 7L, 7L, 8L, 8L, 8L, 8L, 9L, 9L, 9L, 9L, 10L,
10L, 10L, 10L, 11L, 11L, 11L, 11L, 12L, 12L, 12L, 12L, 13L, 13L,
13L, 13L, 14L, 14L, 14L, 14L, 15L, 15L, 15L, 15L, 16L, 16L, 16L,
16L, 17L, 17L, 17L, 17L, 18L, 18L, 18L, 18L, 19L, 19L, 19L, 19L,
20L, 20L, 20L, 20L, 21L, 21L, 21L, 21L, 22L, 22L, 22L, 22L, 23L,
23L, 23L, 23L, 24L, 24L, 24L, 24L), Diam = c(4.32, 4.55, 4.16,
4.44, 4.49, 4.3, 4.52, 4.61, 4.44, 4.32, 4.25, 4.5, 4.55, 4.15,
4.42, 4.49, 4.21, 4.3, 4.29, 4.63, 4.56, 4.61, 4.29, 4.56, 4.63,
4.3, 4.41, 4.58, 4.38, 4.65, 4.43, 4.44, 4.12, 4.49, 4.3, 4.36,
4.45, 4.56, 4.38, 4.51, 4.62, 4.67, 4.7, 4.58, 4.33, 4.23, 4.34,
4.58, 4.29, 4.38, 4.28, 4.41, 4.15, 4.35, 4.28, 4.23, 4.21, 4.3,
4.32, 4.38, 4.16, 4.28, 4.31, 4.59, 4.14, 4.18, 4.08, 4.21, 4.51,
4.2, 4.28, 4.19, 4.1, 4.33, 4.37, 4.47, 3.99, 4.09, 4.47, 4.25,
4.24, 4.54, 4.43, 4.38, 4.23, 4.48, 4.31, 4.57, 4.27, 4.4, 4.32,
4.56, 4.7, 4.65, 4.49, 4.38)), .Names = c("Machine", "Sample",
"Diam"), class = "data.frame", row.names = c(NA, -96L))
Gilênio Borges Fernandes Professor Associado IV (Aposentado) Professor Adjunto A (Substituto) Universidade Federal da Bahia Instituto de Matemática e Estatística Departamento de Estatística Av. Adhemar de Barros, s/n – Ondina. 40.170-110 - Salvador - BA, Brasil Tel.: (071)3283-6340/6341/6337 Fax: (071)3283-6336 Skype: gilenio.fernandes Lattes: http://lattes.cnpq.br/6764860618464860

Segue CMR. library(latticeExtra) Parts$Sample <- factor(Parts$Sample, levels = unique(sort(Parts$Sample))) str(Parts) nlev <- function(x) nlevels(droplevels(x)) resizePanels( bwplot(Diam ~ Sample | Machine, data = Parts, pch = "|", layout = c(NA, 1), scales = list(x = list(relation = "free")), par.settings = list( box.rectangle = list(col = 1, fill = c("gray70")), box.umbrella = list(col = 1, lty = 1))), w = with(Parts, tapply(Sample, Machine, nlev))) À disposição. Walmes.
participantes (2)
-
Gilenio Borges Fernandes
-
Walmes Zeviani