-
Notifications
You must be signed in to change notification settings - Fork 0
/
Copy pathCI_fluxes_hgII_diffexchange_boxes.R
100 lines (80 loc) · 3.27 KB
/
CI_fluxes_hgII_diffexchange_boxes.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
setwd("C:/Users/gi/Desktop/ULTIME_hgII___/In1");
d6<-read.csv(file='exchange_hgII.csv', header=T)
setwd("C:/Users/gi/Desktop/ULTIME_hgII___/In2")
d2<-read.csv(file='exchange_hgII.csv', header=T)
setwd("C:/Users/gi/Desktop/ULTIME_hgII___/IIe_IIIe_IVe/IIIe/In1c")
dn1<-read.csv(file='exchange_hgII.csv', header=T)
setwd("C:/Users/gi/Desktop/ULTIME_hgII___/IIe_IIIe_IVe/IIIe/In2c")
dn2<-read.csv(file='exchange_hgII.csv', header=T)
setwd("C:/Users/gi/Desktop/ULTIME_hgII___/IIe_IIIe_IVe/IIIe/In3c")
dn3<-read.csv(file='exchange_hgII.csv', header=T)
setwd("C:/Users/gi/Desktop/ULTIME_hgII___/IIe_IIIe_IVe/IIIe/In4c")
dn4<-read.csv(file='exchange_hgII.csv', header=T)
setwd("C:/Users/gi/Desktop/ULTIME_hgII___/IIe_IIIe_IVe/IIe/In2b")
dn5<-read.csv(file='exchange_hgII.csv', header=T)
setwd("C:/Users/gi/Desktop/ULTIME_hgII___/IIe_IIIe_IVe/IIe/In1b")
dn6<-read.csv(file='exchange_hgII.csv', header=T)
setwd("C:/Users/gi/Desktop/nuoveMEt3/Fosfati/F6")
nm5<-read.csv(file='exchange_hgII.csv', header=T)
setwd("C:/Users/gi/Desktop/nuoveMEt3/Fosfati/F2")
nm4<-read.csv(file='exchange_hgII.csv', header=T)
setwd("C:/Users/gi/Desktop/nuoveMEt3/add_sim2")
nm2<-read.csv(file='exchange_hgII.csv', header=T)
setwd("C:/Users/gi/Desktop/nuoveMEt3/add_sim3")
nm3<-read.csv(file='exchange_hgII.csv', header=T)
diff_OL<-c(d2$diff_OL_kmol_y,d6$diff_OL_kmol_y,
dn1$diff_OL_kmol_y,dn2$diff_OL_kmol_y,
dn3$diff_OL_kmol_y,
dn4$diff_OL_kmol_y,
dn5$diff_OL_kmol_y,dn6$diff_OL_kmol_y,
nm2$diff_OL_kmol_y,
nm3$diff_OL_kmol_y,nm4$diff_OL_kmol_y,
nm5$diff_OL_kmol_y)
diff_SOL<-c(d2$diff_SOL_kmol_y,d6$diff_SOL_kmol_y,
dn1$diff_SOL_kmol_y,dn2$diff_SOL_kmol_y,
dn3$diff_SOL_kmol_y,
dn4$diff_SOL_kmol_y,dn5$diff_SOL_kmol_y,
dn6$diff_SOL_kmol_y,
nm2$diff_SOL_kmol_y,
nm3$diff_SOL_kmol_y,nm4$diff_SOL_kmol_y,
nm5$diff_SOL_kmol_y)
diff_Sed<-c(d2$diff_Sed_kmol_y,d6$diff_Sed_kmol_y,
dn1$diff_Sed_kmol_y,dn2$diff_Sed_kmol_y,
dn3$diff_Sed_kmol_y,
dn4$diff_Sed_kmol_y,
dn5$diff_Sed_kmol_y,dn6$diff_Sed_kmol_y,
nm2$diff_Sed_kmol_y,
nm3$diff_Sed_kmol_y,nm4$diff_Sed_kmol_y,
nm5$diff_Sed_kmol_y)
adv_to_SOL_hgII<-c(d2$adv_to_SOL,d6$adv_to_SOL,
dn1$adv_to_SOL,dn2$adv_to_SOL,
dn3$adv_to_SOL,
dn4$adv_to_SOL,
dn5$adv_to_SOL,dn6$adv_to_SOL,
nm2$adv_to_SOL,
nm3$adv_to_SOL,nm4$adv_to_SOL,nm5$adv_to_SOL)
SE_advhgII_to_SOL<-sd(adv_to_SOL_hgII)/sqrt(13)
CI_advhgII_to_SOL_min<-mean(adv_to_SOL_hgII)-1.96*SE_advhgII_to_SOL
CI_advhgII_to_SOL_max<-mean(adv_to_SOL_hgII)+1.96*SE_advhgII_to_SOL
mean(diff_OL)
SE_diff_OL<-sd(diff_OL)/sqrt(13)
CI_diff_OL_min<-mean(diff_OL)-1.96*SE_diff_OL
CI_diff_OL_max<-mean(diff_OL)+1.96*SE_diff_OL
SE_diff_SOL<-sd(diff_SOL)/sqrt(13)
SE_diff_SOL
mean(diff_SOL)
SE_diff_Sed<-sd(diff_Sed)/sqrt(length(diff_Sed))
mean(diff_Sed)
SE_diff_Sed
mean(diff_OL);SE_diff_OL
mean(diff_SOL);SE_diff_SOL
mean(adv_to_SOL_hgII);SE_advhgII_to_SOL
mean(diff_Sed);SE_diff_Sed
CI_diff_OL_min
CI_diff_OL_max
CI_diff_SOL_min
CI_diff_SOL_max
CI_diff_Sed_min
CI_diff_Sed_max
CI_advhgII_to_SOL_min
CI_advhgII_to_SOL_max