2

Estou tentando fazer um mapa no R pela primeira vez. Fiz todo o processo necessário com os dados e estou seguindo as instruções de um script disponível aqui mas ele não reconhece o objeto "longest".

Vocês podem ajudar?

O script está assim:

#install.packages("pacman")
rm(list=ls())

pacman::p_load(tidyverse, janitor, haven, forcats, 
       formattable,
       knitr,stringr, hrbrthemes, 
       rgdal, mapproj, plotly)

CP <- read.table("/Users/marciarangelcandido/CP_UF.csv", sep = ";", header = TRUE)

CP <- CP %>% 
  mutate(Estado = case_when(Estado == "A3re" ~ "AC",
                            TRUE ~ as.character(Estado))) %>% 
  rename("sigla" = "Estado")


mapa <- read.csv("https://raw.githubusercontent.com/thiago-ms-cp/programacao_iesp/master/mapa_est_geral.csv")

library(descr)

mapa$sigla <- toUTF8(mapa$sigla, "IBM850")

#palette(c("#779999", "#99bbbb", "#bbdddd", "#ddffff"))
#plot(mapa, col = mapa$sigla)
#title("Mapa de Publicações da Ciência Política")


CP <- CP %>% 
  select(sigla, Estrato)

mapa <- mapa %>%
  left_join(CP, by = "sigla")

mapa<-mapa %>%
  group_by(sigla) %>% 
  mutate(a1= sum(Estrato==10)) %>%  
  mutate(a2=sum(Estrato==9)) %>%
  mutate(b1=sum(Estrato==8)) %>%
  mutate(b2=sum(Estrato==7)) %>%
  mutate(b3=sum(Estrato==6)) %>%
  mutate(b4=sum(Estrato==5)) %>%
  mutate(b5=sum(Estrato==4)) %>%
  mutate(c=sum(Estrato==3))

mapa$texto <- paste("A1 =", CP$a1, "revistas",
                  "; A2 =", CP$a2, "revistas",
                  "; B1 =", CP$b1, "revistas",
                  "; B2 =", CP$b2, "revistas",
                  "; B3 =", CP$b3, "revistas",
                  "; B4 =", CP$b4, "revistas",
                  "; B5 =", CP$b5, "revistas",
                  "; C =", CP$c, "revistas") 


mapa %>%
  plot_ly(~longest, ~latest, group = ~group)


plot_ly(mapa, lon = **longest**, lat = latest, text = texto,
        color = Estrato, type = 'scattergeo', locationmode = 'country names') %>%
  layout(title = 'Populations<br>(Click legend to toggle)')


mapa$Estado <- cut(mapa$Estado, 4)
levels(mapa$Estado)
1
  • 3
    A dificuldade é que seu exemplo não é reproduzível. A base CP_UF não foi disponibilizada. Além disso, sugiro você dar uma lida na documentação do pacote sf. Vai facilitar muito a sua vida na manipulação de mapas. O plotly aceita simple features (sf) e a maioria dos htmlwidgets também.
    – José
    13/12/2017 às 18:29

1 Resposta 1

0

Se a ideia é meramente criar um mapa do Brasil, o seguinte código o produz de maneira bem simples:

library(geobr)
library(ggplot2)
library(cowplot)

uf <- read_state(code_state = "all")

no_axis <- theme(axis.title=element_blank(),
                 axis.text=element_blank(),
                 axis.ticks=element_blank())

plot_grafico <- ggplot() + 
  geom_sf(data=uf, fill="white", color="black", size=.15, show.legend = FALSE) +
  theme_nothing() +
  no_axis +
  labs(subtitle="Terras indígenas", size=8)

plot_grafico

A função no_axis é meramente para remover os eixos e deixar a imagem simples. O pacote geo_br permite produzir diversos mapas. Recomendo explorar as funcionalidades.

Esta não é a resposta que você está procurando? Pesquise outras perguntas com a tag .