-
Notifications
You must be signed in to change notification settings - Fork 1
/
Copy path5-diversity_ubi_compl_essays.R
109 lines (95 loc) · 3.01 KB
/
5-diversity_ubi_compl_essays.R
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
102
103
104
105
106
107
108
109
#load: dplyr, EconGeo
#import data:
all_function<- read.csv("/Users/aferloni/switchdrive/Andrea-Celine/These Andrea/Paper 2 - Urban relatedness/DATA/lur_data_function.csv", header = T, stringsAsFactors = F)
nodes<- read.csv("/Users/aferloni/switchdrive/Andrea-Celine/These Andrea/Paper 2 - Urban relatedness/DATA/nodes.csv", header = T, stringsAsFactors = F)
#calculate diversity
#prepare data----
#distinguish periods
period1<- all_function%>%
filter(time==1)%>%
distinct(lur, group_id, count)%>%
filter(group_id %in% nodes$id)%>%
arrange(lur)
period2<- all_function%>%
filter(time==2)%>%
distinct(lur, group_id, count)%>%
filter(group_id %in% nodes$id)%>%
arrange(lur)
period3<- all_function%>%
filter(time==3)%>%
distinct(lur, group_id, count)%>%
filter(group_id %in% nodes$id)%>%
arrange(lur)
period4<- all_function%>%
filter(time==4)%>%
distinct(lur, group_id, count)%>%
filter(group_id %in% nodes$id)%>%
arrange(lur)
#create a city industry matrix
period1_mat<- get.matrix(period1)
period2_mat<- get.matrix(period2)
period3_mat<- get.matrix(period3)
period4_mat<- get.matrix(period4)
#calculate diversity, ubiquity, kci----
#diversity
div1<- as.data.frame(EconGeo::diversity(period1_mat, RCA=TRUE))%>%
rownames_to_column("lur")%>%
rename(diversity=2)%>%
mutate(time=1)
div2<- as.data.frame(EconGeo::diversity(period2_mat, RCA=TRUE))%>%
rownames_to_column("lur")%>%
rename(diversity=2)%>%
mutate(time=2)
div3<- as.data.frame(EconGeo::diversity(period3_mat, RCA=TRUE))%>%
rownames_to_column("lur")%>%
rename(diversity=2)%>%
mutate(time=3)
div4<- as.data.frame(EconGeo::diversity(period4_mat, RCA=TRUE))%>%
rownames_to_column("lur")%>%
rename(diversity=2)%>%
mutate(time=4)
div_all<- rbind(div1, div2, div3, div4)
#calculate ubiquity
ubi1<- as.data.frame(ubiquity(period1_mat, RCA=TRUE))%>%
rownames_to_column("code")%>%
rename(ubiquity=2)%>%
mutate(time=1)
ubi2<- as.data.frame(ubiquity(period2_mat, RCA=TRUE))%>%
rownames_to_column("code")%>%
rename(ubiquity=2)%>%
mutate(time=2)
ubi3<- as.data.frame(ubiquity(period3_mat, RCA=TRUE))%>%
rownames_to_column("code")%>%
rename(ubiquity=2)%>%
mutate(time=3)
ubi4<- as.data.frame(ubiquity(period4_mat, RCA=TRUE))%>%
rownames_to_column("code")%>%
rename(ubiquity=2)%>%
mutate(time=4)
ubi_all<- rbind(ubi1, ubi2, ubi3, ubi4)
#KCI
kci<- as.data.frame(KCI(period1_mat, RCA=TRUE))
rownames(kci)<- rownames(period1_mat)
kci<- kci%>%
rownames_to_column("lur")%>%
rename(kci=2)%>%
mutate(time=1)
kci2<- as.data.frame(KCI(period2_mat, RCA=TRUE))
rownames(kci2)<- rownames(period2_mat)
kci2<- kci2%>%
rownames_to_column("lur")%>%
rename(kci=2)%>%
mutate(time=2)
kci3<- as.data.frame(KCI(period3_mat, RCA=TRUE))
rownames(kci3)<- rownames(period3_mat)
kci3<- kci3%>%
rownames_to_column("lur")%>%
rename(kci=2)%>%
mutate(time=3)
kci4<- as.data.frame(KCI(period4_mat, RCA=TRUE))
rownames(kci4)<- rownames(period4_mat)
kci4<- kci4%>%
rownames_to_column("lur")%>%
rename(kci=2)%>%
mutate(time=4)
kciall<- rbind(kci, kci2, kci3, kci4)