-
Notifications
You must be signed in to change notification settings - Fork 1
/
Copy path.Rhistory
101 lines (101 loc) · 2.96 KB
/
.Rhistory
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
knitr::opts_chunk$set(echo = TRUE)
library(sf)
library(tidyverse)
library(terra)
library(here)
danang_shp <- st_read(here::here("data", "gadm41_VNM_shp", "gadm41_VNM_1.shp"))
View(danang_shp)
# load shapefile
# source: https://gadm.org/download_country.html
vietnam_shp <- st_read(here::here("data", "gadm41_VNM_shp", "gadm41_VNM_1.shp"))
# extract Da Nang
danang_shp <- vietnam_shp %>%
dplyr::filter(VARNAME_1 == "Da Nang")
View(danang_shp)
plot(danang_shp)
pop_2020 <- terra::raster(here::here("data", "pop_data", "vnm_ppp_2020_UNadj.tif"))
pop_2020 <- terra::rast(here::here("data", "pop_data", "vnm_ppp_2020_UNadj.tif"))
pop_2000 <- terra::rast(here::here("data", "pop_data", "vnm_ppp_2000_UNadj.tif"))
pop_2005 <- terra::rast(here::here("data", "pop_data", "vnm_ppp_2005_UNadj.tif"))
pop_2010 <- terra::rast(here::here("data", "pop_data", "vnm_ppp_2010_UNadj.tif"))
pop_2015 <- terra::rast(here::here("data", "pop_data", "vnm_ppp_2015_UNadj.tif"))
pop_2020 <- terra::rast(here::here("data", "pop_data", "vnm_ppp_2020_UNadj.tif"))
pop_2000
pop_2020
pop <- c(pop_2000, pop_2005, pop_2010, pop_2015, pop_2020)
pop
names(pop)
danang_pop <- pop %>%
terra::crop(.) %>%
terra::mask(.)
danang_pop <- pop %>%
terra::crop(., danang_shp) %>%
terra::mask(., danang_shp)
danang_pop
plot(danang_pop)
View(pop_2000)
sum(danang_pop$vnm_ppp_2000_UNadj)
sum(danang_pop$vnm_ppp_2000_UNadj)
terra::sum(danang_pop$vnm_ppp_2000_UNadj)
population <- danang_pop$vnm_ppp_2020_UNadj %>%
as_data_frame()
View(population)
population <- danang_pop$vnm_ppp_2020_UNadj %>%
as_data_frame() %>%
sum()
population <- danang_pop %>%
as_data_frame() %>%
sum()
population <- danang_pop$vnm_ppp_2020_UNadj %>%
as_data_frame() %>%
sum()
danang_pop$vnm_ppp_2020_UNadj %>%
as_data_frame() %>%
sum()
danang_pop$vnm_ppp_2000_UNadj %>%
as_data_frame() %>%
sum()
population <- pop %>%
sum(., par = TRUE)
population <- pop %>%
as_tibble()
gc()
population <- pop %>%
as_data_frame()
xaringan:::inf_mr()
xaringan:::inf_mr()
xaringan:::inf_mr()
xaringan:::inf_mr()
options(htmltools.dir.version = FALSE)
library(RefManageR)
BibOptions(check.entries = FALSE,
bib.style = "authoryear",
cite.style = "authoryear",
style = "markdown",
hyperlink = TRUE,
dashed = FALSE,
no.print.fields=c("doi", "url", "urldate", "issn"))
myBib <- ReadBib("./Bib.bib", check = FALSE)
library(knitcitations)
knitr::include_graphics('img/timeline.png')
knitr::include_graphics('img/timeline.png')
```
knitr::include_graphics('img/timeline.png')
```
knitr::include_graphics('img/timeline.png')
knitr::include_graphics('img/timeline.png')
knitr::include_graphics('img/timeline.png')
cost <- read.csv(file = costs)
xaringan:::inf_mr()
options(htmltools.dir.version = FALSE)
library(RefManageR)
BibOptions(check.entries = FALSE,
bib.style = "authoryear",
cite.style = "authoryear",
style = "markdown",
hyperlink = TRUE,
dashed = FALSE,
no.print.fields=c("doi", "url", "urldate", "issn"))
myBib <- ReadBib("./Bib.bib", check = FALSE)
library(knitcitations)
knitr::include_graphics('img/venn_diagram2.png')