Pessoal plotei o mapa de MG, gostaria de selecionar e colocar em cima só as cidades da zona da mata de MG! São em torno de 75 cidades das 853. Estou usando o comando subset. Tem algo mais simples de selecionar estas 75 cidades?

require(rgdal)
mg=readOGR("31MUE250GC_SIR.shp",layer="31MUE250GC_SIR",encoding = "latin1")
plot(mg,axes=T)
zmata=subset(mg, NM_MUNICIP%in%c("ABRE CAMPO","ACAIACA","ALVINOPOLIS","AMPARO DO SERRA"))
plot(zmata,col="blue",add=TRUE)

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