HW7

president<-read.csv('./1976-2020-president.csv')
#Clean Data
dt<-as.data.table(president)
data<-dt[ ,list( mostvote=max(candidatevotes), candidate=candidate, vote=candidatevotes, party=party_simplified, state_po=state_po, state_fips=state_fips), by=.(year, state)]
data$winner<-ifelse(data$vote==data$mostvote,1,NA)
data<-na.omit(data)
data$color=0
for(i in 1:nrow(data)){
  if(data$party[i]=="DEMOCRAT"){
    data$color[i]=1
  }
    if(data$party[i]=="REPUBLICAN"){
    data$color[i]=2
    }
    if(data$party[i]=="OTHER"){
    data$color[i]=3
  }
}
data <- mutate(data,color = factor(color))

data2<-data[which(year==2012)]

plot_usmap(data=data2, values= "color", color='grey')+ 
scale_fill_manual(values = c(`1` = "blue", `2` = "red", `3`='yellow'), name = "Presidential Election Results by Year", labels=c(`1`="Democrat", `2`="Republican", `3`="Other")) + theme(legend.position = "right" , plot.title =element_text(color="black", size=18, face="bold.italic", hjust=.5))+ggtitle("2012")

plot_usmap(data=data, values= "color", color='grey')+ 
scale_fill_manual(values = c(`1` = "blue", `2` = "red", `3`='yellow'), name = "Presidential Election Results by State and Year", labels=c(`1`="Democrat", `2`="Republican", `3`="Other")) + theme(legend.position = "right" , plot.title =element_text(color="black", size=18, face="bold.italic", hjust=.5)) + facet_wrap(facets = vars(year))

Part 2: Interactive Maps with Leaflet

myMap <- leaflet() %>%
addProviderTiles(providers$OpenStreetMap)
myMap

myMap <- leaflet() %>%
addProviderTiles(providers$Esri.NatGeoWorldMap) %>%
  setView(lat=33.947474, lng=-83.373671, zoom = 12)

myMap

sp <- st_read("G:/My Drive/Seans Drive/PhD/Classes/AQME/HW/220130_rrc_outline_block_al2/T220130_RRC_Outline_Block_AL2.shp")

Reading layer T220130_RRC_Outline_Block_AL2' from data sourceG:DriveDrive\220130_rrc_outline_block_al2220130_RRC_Outline_Block_AL2.shp’ using driver `ESRI Shapefile’ Simple feature collection with 178 features and 15 fields Geometry type: MULTIPOLYGON Dimension: XY Bounding box: xmin: 92.13366 ymin: 20.91856 xmax: 92.26863 ymax: 21.22292 Geodetic CRS: WGS 84

campShapeFile <- st_transform(sp, CRS("+proj=longlat +datum=WGS84 +no_defs"))

myMap <- leaflet() %>%
  addProviderTiles(providers$Esri.NatGeoWorldMap) %>%
  addPolygons(data=campShapeFile, fill=TRUE, stroke=T, weight=1, highlight = highlightOptions(fillOpacity = 0.7), label = campShapeFile$OBJECTID_1) %>%
  addWMSTiles(
"http://mesonet.agron.iastate.edu/cgi-bin/wms/nexrad/n0r.cgi",
layers = "nexrad-n0r-900913",
options = WMSTileOptions(format = "image/png", transparent = TRUE),
attribution = "Weather data © 2012 IEM Nexrad"
)

  
myMap