Bom, existe alguma forma mais elegante de fazer esta ordenação? Pode ser por fator e população ou fator, região e população. Tentei sem sucesso com. dados[order(dados$fator,dados$pop2007, decreasing=T),] dput(dados) structure(list(fator = structure(c(2L, 10L, 12L, 5L, 6L, 8L, 11L, 9L, 13L, 1L, 7L, 3L, 4L), .Label = c("Campinas", "Feira de Santana", "Guarulhos", "Jaboatão" , "Niterói" , "Nova Iguaçu" , "Osasco", "Recife", "Rio de Janeiro", "Salvador", "Santo André" , "São Paulo" , "Sorocaba"), class = "factor"), região = structure(c(1L, 1L, 2L, 2L, 2L, 1L, 2L, 2L, 2L, 2L, 2L, 2L, 1L), .Label = c(" NE", " SE"), class = "factor"), estado = structure(c(1L, 1L, 5L, 3L, 3L, 2L, 4L, 3L, 4L, 4L, 4L, 4L, 2L), .Label = c(" BA", " PE", " RJ", " SP", " SP "), class = "factor"), pop2007 = c(544113L, 2714119L, 11104712L, 476669L, 858150L, 1528970L, 676846L, 6178762L, 590846L, 1073020L, 724368L, 1289047L, 661901L)), .Names = c("fator", "região", "estado", "pop2007"), class = "data.frame", row.names = c(NA, -13L)) dados[order(dados$estado,dados$pop2007,decreasing=T),] fator região estado pop2007 9 Sorocaba SE SP 590846 3 São Paulo SE SP 11104712 7 Santo André SE SP 676846 2 Salvador NE BA 2714119 8 Rio de Janeiro SE RJ 6178762 6 Recife NE PE 1528970 11 Osasco SE SP 724368 5 Nova Iguaçu SE RJ 858150 4 Niterói SE RJ 476669 13 Jaboatão NE PE 661901 12 Guarulhos SE SP 1289047 1 Feira de Santana NE BA 544113 10 Campinas SE SP 1073020 André Oliveira Souza. Graduação em Matemática, mestrado em estatística aplicada.Instituto Federal de Educação, Ciência e Tecnologia do Espirito Santo. IFES