Plotar mapa por fator categórico

Pessoal estou tentando plotar um mapa por um fator (X) e por uma densidade (Y). O gráfico pela variável contínua (Y) está bom, mas pela variável categórica (X) dá os avisos, o gráfico não sai colorido e plota preto e branco com aviso. Já busquei na lista e no são google, mas não consegui resolver. O arquivo é um .shp do Mato Grosso que busquei no IBGE.require(rgdal) mt=readOGR("51MUE250GC_SIR.shp",layer="51MUE250GC_SIR",encoding = "latin1") n <- nrow(mt@data) mt@data$Y <- rnorm(n)head(mt@data) ID CD_GEOCODM NM_MUNICIP Y X 0 218 5100102 ACORIZAL -2.1099611 D 1 219 5100201 ÁGUA BOA 0.9998267 D 2 220 5100250 ALTA FLORESTA -0.7521304 D 3 221 5100300 ALTO ARAGUAIA -1.5353345 B 4 222 5100359 ALTO BOA VISTA 1.1464625 C 5 223 5100409 ALTO GARÇAS -1.8388984 D spplot(mt,"Y")mt@data$X <- sample(LETTERS[1:4], size = n, replace = T)spplot(mt, "X", col.regions = heat.colors(4))Houve 32 avisos (use warnings() para vê-los)warnings() Mensagens de aviso: 1: In Ops.factor(z, at[i]) : ‘>=’ não faz sentido para fatores 2: In Ops.factor(z, at[i + 1]) : ‘<’ não faz sentido para fatores 3: In Ops.factor(z, at[i]) : ‘>=’ não faz sentido para fatores 4: In Ops.factor(z, at[i + 1]) : ‘<’ não faz sentido para fatores 5: In Ops.factor(z, at[i]) : ‘>=’ não faz sentido para fatores 6: In Ops.factor(z, at[i + 1]) : ‘<’ não faz sentido para fatores 7: In Ops.factor(z, at[i]) : ‘>=’ não faz sentido para fatores 8: In Ops.factor(z, at[i + 1]) : ‘<’ não faz sentido para fatores 9: In Ops.factor(z, at[i]) : ‘>=’ não faz sentido para fatores 10: In Ops.factor(z, at[i + 1]) : ‘<’ não faz sentido para fatores 11: In Ops.factor(z, at[i]) : ‘>=’ não faz sentido para fatores 12: In Ops.factor(z, at[i + 1]) : ‘<’ não faz sentido para fatores 13: In Ops.factor(z, at[i]) : ‘>=’ não faz sentido para fatores 14: In Ops.factor(z, at[i + 1]) : ‘<’ não faz sentido para fatores 15: In Ops.factor(z, at[i]) : ‘>=’ não faz sentido para fatores 16: In Ops.factor(z, at[i + 1]) : ‘<’ não faz sentido para fatores 17: In Ops.factor(z, at[i]) : ‘>=’ não faz sentido para fatores 18: In Ops.factor(z, at[i + 1]) : ‘<’ não faz sentido para fatores 19: In Ops.factor(z, at[i]) : ‘>=’ não faz sentido para fatores 20: In Ops.factor(z, at[i + 1]) : ‘<’ não faz sentido para fatores 21: In Ops.factor(z, at[i]) : ‘>=’ não faz sentido para fatores 22: In Ops.factor(z, at[i + 1]) : ‘<’ não faz sentido para fatores 23: In Ops.factor(z, at[i]) : ‘>=’ não faz sentido para fatores 24: In Ops.factor(z, at[i + 1]) : ‘<’ não faz sentido para fatores 25: In Ops.factor(z, at[i]) : ‘>=’ não faz sentido para fatores 26: In Ops.factor(z, at[i + 1]) : ‘<’ não faz sentido para fatores 27: In Ops.factor(z, at[i]) : ‘>=’ não faz sentido para fatores 28: In Ops.factor(z, at[i + 1]) : ‘<’ não faz sentido para fatores 29: In Ops.factor(z, at[i]) : ‘>=’ não faz sentido para fatores 30: In Ops.factor(z, at[i + 1]) : ‘<’ não faz sentido para fatores 31: In Ops.factor(z, at[i]) : ‘>=’ não faz sentido para fatores 32: In Ops.factor(z, at[i + 1]) : ‘<’ não faz sentido para fatores Alguém poderia me ajudar no erro? obrigado 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

Andre, bom dia! Deve funcionar se transformar "X" em factor(). mt@data$X <- as.factor(sample(LETTERS[1:4], size = n, replace=T)) Éder Comunello <c <comunello.eder@gmail.com>omunello.eder@gmail.com> Dourados, MS - [22 16.5'S, 54 49'W]

Não sei usar a função spplot, mas pela mensagem de erro, será que não seria o caso de transformar a variável em um fator ordenado? mt@data$X <- as.ordered(sample(LETTERS[1:4], size = n, replace=T)) Leonardo Ferreira Fontenelle[1] Em Dom 19 abr. 2015, às 09:11, Éder Comunello escreveu:
Andre, bom dia!
Deve funcionar se transformar "X" em factor().
mt@data$X <- as.factor(sample(LETTERS[1:4], size = n, replace=T))
Éder Comunello <c[2]omunello.eder@gmail.com> Dourados, MS - [22 16.5'S, 54 49'W]
_________________________________________________ R-br mailing list R-br@listas.c3sl.ufpr.br https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código mínimo reproduzível.
Links: 1. http://lattes.cnpq.br/9234772336296638 2. mailto:comunello.eder@gmail.com

Leonardo, boa tarde! Fiquei na dúvida quanto ao uso de as.ordered() e rodei um pequeno teste. Parece funcionar do mesmo modo que com as.factor(). ### <code r> getwd() url0 <- " ftp://geoftp.ibge.gov.br/malhas_digitais/municipio_2007/escala_2500mil/proj_geografica_sirgas2000/uf/mt/51mu2500gsr.zip " bName <- sub(".zip$", "", basename(url0)) setInternet2(FALSE) ### download a partir de ftp no RStudio Windows if (!file.exists(basename(url0))) download.file(url0, dest=basename(url0), mode="wb") unzip(basename(url0), list=T) ### conteúdo do zip unzip(basename(url0)) require(rgdal) mt <- readOGR(dsn=paste0(bName,".shp"), layer=bName, encoding="latin1") str(mt, max=2) plot(mt) head(mt@data) n <- nrow(mt@data); n set.seed(756); x <- sample(LETTERS[1:3], size = n, replace=T) mt@data$X1 <- as.factor(x) mt@data$X2 <- as.ordered(x) str(mt@data) head(mt@data) spplot(mt, "X1", col.regions=heat.colors(3)) spplot(mt, "X2", col.regions=heat.colors(3)) ### </code> Éder Comunello <c <comunello.eder@gmail.com>omunello.eder@gmail.com> Dourados, MS - [22 16.5'S, 54 49'W]

Que bom, então! Leonardo Ferreira Fontenelle[1] Em Ter 21 abr. 2015, às 17:44, Éder Comunello escreveu:
Leonardo, boa tarde!
Fiquei na dúvida quanto ao uso de as.ordered() e rodei um pequeno teste. Parece funcionar do mesmo modo que com as.factor().
### <code r> getwd() url0 <- "ftp://geoftp.ibge.gov.br/malhas_digitais/municipio_2007/escala_2500mil/proj_geografica_sirgas2000/uf/mt/51mu2500gsr.zip" bName <- sub(".zip$", "", basename(url0))
setInternet2(FALSE) ### download a partir de ftp no RStudio Windows if (!file.exists(basename(url0))) download.file(url0, dest=basename(url0), mode="wb") unzip(basename(url0), list=T) ### conteúdo do zip unzip(basename(url0))
require(rgdal) mt <- readOGR(dsn=paste0(bName,".shp"), layer=bName, encoding="latin1") str(mt, max=2)
plot(mt) head(mt@data) n <- nrow(mt@data); n
set.seed(756); x <- sample(LETTERS[1:3], size = n, replace=T) mt@data$X1 <- as.factor(x) mt@data$X2 <- as.ordered(x)
str(mt@data) head(mt@data) spplot(mt, "X1", col.regions=heat.colors(3)) spplot(mt, "X2", col.regions=heat.colors(3)) ### </code>
Éder Comunello <c[2]omunello.eder@gmail.com> Dourados, MS - [22 16.5'S, 54 49'W]
_________________________________________________ R-br mailing list R-br@listas.c3sl.ufpr.br https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código mínimo reproduzível.
Links: 1. http://lattes.cnpq.br/9234772336296638 2. mailto:comunello.eder@gmail.com
participantes (3)
-
Andre Oliveira
-
Leonardo Ferreira Fontenelle
-
Éder Comunello