w

Author: JCfly3000

.DS_Store

@@ -393,7 +393,7 @@ p+ theme_fivethirtyeight()
 # convert ggplot to plotly
 
 ```{r}
-p=ggplot(tips, aes(tip, total_bill,color=sex)) + geom_point()+ scale_x_continuous(name="new x name")+ scale_y_continuous(name="new y name")
+p=ggplot(data=tips, aes(tip, total_bill,color=sex)) + geom_point()+ scale_x_continuous(name="new x name")+ scale_y_continuous(name="new y name")
 
 pp=ggplotly(p)
 
@@ -461,6 +461,92 @@ p
 ```
 
 
+# ggplot2 extension packages
+
+https://www.youtube.com/watch?v=efj6-aawubs&list=PLBnFxG6owe1FcJnrNbaIJyALHMyoKW27T
+
+## patchwork
+
+The goal of patchwork is to make it ridiculously simple to combine separate ggplots into the same graphic. 
+
+```{r}
+#pak::pak('patchwork')
+```
+
+```{r}
+library(ggplot2)
+library(patchwork)
+
+p1 <- ggplot(mtcars) + geom_point(aes(mpg, disp))
+p2 <- ggplot(mtcars) + geom_boxplot(aes(gear, disp, group = gear))
+
+p3 <- ggplot(mtcars) + geom_smooth(aes(disp, qsec))
+p4 <- ggplot(mtcars) + geom_bar(aes(carb))
+
+(p1 | p2 | p3) /
+      p4
+```
+
+## legendry
+
+
+The goal of legendry is to provide additional guide functionality to the ggplot2 ecosystem.
+
+
+```{r}
+#pak::pak('legendry')
+```
+
+```{r}
+library(ggplot2)
+library(legendry)
+```
+
+
+```{r}
+
+base <- ggplot(mpg, aes(displ, hwy, colour = cty)) +
+  geom_point() +
+  labs(
+    x = "Engine displacement",
+    y = "Highway miles per gallon",
+    col = "City miles\nper gallon"
+  ) +
+  theme(axis.line = element_line())
+```
+
+
+```{r}
+# A partial guide to display a bracket
+efficient_bracket <- primitive_bracket(
+  # Keys determine what is displayed
+  key = key_range_manual(start = 25, end = Inf, name = "Efficient"),
+  bracket = "square",
+  # We want vertical text
+  theme = theme(
+    legend.text = element_text(angle = 90, hjust = 0.5),
+    axis.text.y.left = element_text(angle = 90, hjust = 0.5)
+  )
+)
+
+base + guides(y = guide_axis_stack("axis", efficient_bracket))
+```
+
+
+## marquee
+
+## ggiraph
+
+## geomtextpath
+
+## ggpattern
+
+## ggbump
+
+## gghighlight
+
+## ggdist
+
 
 
 # reference:

Plot/.DS_Store

@@ -18,8 +18,7 @@ format:
     code-block-border-left: "#31BAE9"
     code-copy: true
 ---
-
-# data
+# US gdp by State
 
 
 ```{r}
@@ -37,77 +36,43 @@ library(maps)
 mapStates = map("state", fill = TRUE, plot = FALSE)
 ```
 
+US location data
 ```{r}
 json_data=read_sf("us-states.json")
 ```
 
-
 ```{r}
 map_df <- as.data.frame(json_data)
 ```
 
 
-
-# map
-
-```{r}
-data001=read_csv('./data/us shooting dagta001.csv')
-```
-
-
-
 ```{r}
-data001=sample_n(data001,10) %>% clean_names() %>% mutate(full_adress=paste(address,city_or_county,state))
+glimpse(mapStates)
 ```
 
-```{r}
-glimpse(data001)
-```
+US 2-22 GDP data
 
-# search all address to latitude, longitude
 ```{r}
-# geocode the addresses
-data002 <- data001 %>%
-  geocode(full_adress, method = 'arcgis', lat = latitude , long = longitude)
-
-glimpse(data002)
-```
+library(openxlsx)
+library(readxl)
 
-```{r}
-data003=data002 %>% filter(is.na(latitude)==FALSE,(is.na(longitude)==FALSE))
+per_gdp_usd=read_excel('US state gpd 2022.xlsx') %>% mutate(names=str_replace(names,'\\*','') %>% str_trim() )
 ```
 
-# map
 
 ```{r}
-m <- leaflet(data=mapStates) %>%
-   # Add default OpenStreetMap map tiles
-  addTiles() %>%  
-  # add markers
-  addMarkers( lng=data003$longitude, lat=data003$latitude, popup=data003$state)%>% 
-  addPolygons(fillColor = topo.colors(10, alpha = NULL), stroke = FALSE)
-
-m  
+glimpse(per_gdp_usd)
 ```
 
 
-# US gdp by State
 
 ```{r}
-glimpse(mapStates)
-```
-
-```{r}
-library(openxlsx)
-library(readxl)
-
-per_gdp_usd=read_excel('US state gpd 2022.xlsx') %>% mutate(names=str_replace(names,'\\*','') %>% str_trim() )
+map_df002 <- merge(map_df, per_gdp_usd, by.x = "name", by.y = "names")
 ```
 
 
-
 ```{r}
-map_df002 <- merge(map_df, per_gdp_usd, by.x = "name", by.y = "names")
+glimpse(map_df002)
 ```
 
 
@@ -117,14 +82,21 @@ map_sf002 <- sf::st_as_sf(map_df002, sf_column_name = "geometry")
 
 ```{r}
 
-pal <- colorNumeric("magma", NULL)
+#pal <- colorNumeric("magma", NULL)
+pal <- colorNumeric(
+  palette = "Blues",
+  domain = map_df002$per_gdp_usd
+  )
 
+binpal <- colorBin("Blues",map_df002$per_gdp_usd, 20, pretty = FALSE)
 
 m <- leaflet(map_sf002) %>%
    # Add default OpenStreetMap map tiles
   addTiles() %>%  
    addPolygons(smoothFactor = 0.3, fillOpacity = 0.5,weight = 1,
-    fillColor = ~ pal(per_2022_usd)
+    #fillColor = ~ pal(per_2022_usd)
+    #color = ~pal(per_2022_usd)
+   color = ~binpal(per_2022_usd)
     ,popup = ~ paste0(
       "地区：", name, "<br/>",
       "<hr/>",

Plot/1 ggplot2.qmd

@@ -0,0 +1,110 @@
+---
+title: "World map with GPD per capita"
+subtitle: "with leaflet"
+author: "Tony Duan"
+
+execute:
+  warning: false
+  error: false
+
+format:
+  html:
+    toc: true
+    toc-location: right
+    code-fold: show
+    code-tools: true
+    number-sections: true
+    code-block-bg: true
+    code-block-border-left: "#31BAE9"
+    code-copy: true
+---
+
+# Download data
+
+```{r}
+library(tidyverse)
+library(rvest)
+library(chromote)
+library(janitor)
+library(tidygeocoder)
+library(leaflet)
+library(sf)
+```
+
+
+world geojson
+
+```{r}
+# From http://data.okfn.org/data/datasets/geo-boundaries-world-110m
+#countries <- sf::read_sf("https://rstudio.github.io/leaflet/json/countries.geojson")
+countries <- sf::read_sf("https://raw.githubusercontent.com/johan/world.geo.json/refs/heads/master/countries.geo.json")
+
+
+```
+
+world gpd per capita data
+
+```{r}
+world_gdp_per_capita=read.csv("https://gist.githubusercontent.com/vr-23/d6a4a0aadcf3a2640091ca43c25e1955/raw/57b9cb74b24dc87b489c74cc4e11209cd5720642/world-data-2023.csv")
+```
+
+
+```{r}
+# data source: https://data.worldbank.org/
+download.file("https://api.worldbank.org/v2/en/indicator/NY.GDP.PCAP.CD?downloadformat=csv",destfile='world_data.zip')
+```
+```{r}
+unzip('world_data.zip')
+```
+
+```{r}
+world_gdp_per_capita=read.csv("API_NY.GDP.PCAP.CD_DS2_en_csv_v2_76317.csv",skip = 3) |> clean_names()
+```
+
+
+
+```{r}
+data001=countries |> left_join(world_gdp_per_capita,by = join_by(id ==country_code))
+```
+
+
+
+```{r}
+pal <- colorNumeric(palette = c('red','orange','yellow','pink','blue', 'green'),
+                    domain = data001$x2023
+  )
+
+
+
+m <- leaflet(data001) %>%
+   # Add default OpenStreetMap map tiles
+  addTiles() %>%  
+   addPolygons(smoothFactor = 0.3, fillOpacity = 0.5,weight = 1,
+    fillColor = ~ pal(x2023)
+    ,color = ~pal(x2023)
+   #color = ~binpal(x2023)
+    ,popup = ~ paste0(
+      "地区：", name, "<br/>",
+      "<hr/>",
+      "人均gpd：", x2023, "（美元）", "<br/>"
+      
+    )
+     ,label = lapply(paste0(
+      "地区：", "<b>", data001$country_name, "</b>", "<br/>",
+      "人均gpd 美元：", round(data001$x2023), " USD<br/>"
+     # , "总gpd 美元：", round(data001$total_2022_usd), "M <br/>"
+
+    ), htmltools::HTML)
+    )%>% addLegend(
+    position = "bottomright", title = "人均gpd(美元)",
+    pal = pal, values = ~x2023, opacity = 1.0
+    )
+```
+
+
+```{r}
+m  |>  setView(lng = -10, lat = 45, zoom = 2)
+```
+
+# resouce:
+