ohi.model.pressures: Ocean Health Index: Pressures Model
In bbest/ohigui: Ocean Health Index - Graphical User Interface
Description
Usage
Arguments
Details
Value
See Also
Examples
Description
The pressures model function computes a pressures score
for each region given a weighting matrix for a goal and
the individual pressures values.
Usage
1
ohi.model.pressures(p, w, GAMMA = 0.5, browse = F)
Arguments
p
the pressures value matrix [region_id x
pressure]. Each score must be a real number between 0 and
1 inclusive, or NA. The pressure names must be of the
form category_pressure where
category is one of the categories listed in
ohi.pressure.category. Use ss to denote the
social category.
pressure region_id cc_acid cc_sst cc_uv
fp_art_hb 1 0.879 0.360 0.764 NA 2 0.579 0.396 0.531 NA 3
0.926 0.235 0.769 NA 4 0.914 0.554 0.795 NA 5 0.860 0.609
0.802 0.001 6 0.871 0.325 0.788 0.001 7 0.846 0.410 0.677
0.000 8 0.806 0.671 0.752 NA 9 0.844 0.595 0.678 NA 10
0.860 0.575 0.781 0.109
w
the weighting matrix of the form [region_id x
pressure]. Each rank weight must be a real number between
0 and 3 inclusive, or NA.
pressure
region_id cc_acid cc_sst cc_uv fp_art_hb 1 2 1 0.6 NA 2 2
1 0.5 NA 3 2 1 2.1 NA 4 2 1 3.0 NA 5 2 1 2.8 1 6 2 1 2.2
1 7 2 1 1.3 1 8 2 1 1.7 NA 9 2 1 3.0 NA 10 2 1 1.2 1
GAMMA
Multiplier used to combine environmental and
social pressures.
Details
Each pressure layer p(i,j) is either environmental
or social, belongs to a pressures category K \in
\{cc,fp,hd,po,sp,ss\}, and has a value (0..1) for
each region i and pressures layer j. Each
goal has a weight matrix w that has a rank weight
between 0 and 3 inclusive, or NoData, for each region
i and each pressure layer j on a per goal
g basis.
The pressures scores calculations go through 5 steps,
using a complex weighting scheme that varies across
goals, subgoals, pressures categories, and regions:
-
g is the goal or subgoal (e.g., AO, CW, LIV,
ECO, ...),
-
i is the region (e.g., 1, 2, 3,
...),
-
j is the pressures layer or stressor
(e.g., cc_acid, fp_art_lb, etc.).
Calculations
Apply weights for each goal g, region
i, and pressure layer j: Each weighted
pressure p_w(g,i,j) is the pressure layer value
p(i,j) per region i and pressure layer
j multiplied by the rank weight w(g,i,j) for
that goal g, region i, and pressure layer
j. If the w(g,i,j) is NoData or 0, the
weighted pressure p_w(g,i,j) is NoData.
p_{w}(g,i,j) = w(g,i,j) * p(i,j)
Category-level aggregation: The pressures category
score p_K is the sum of all p_w within each
category, then rescaled to 0..1 using a linear scale
range transformation (from 0..3 to 0..1). Any score
p_K greater than 1 is capped to 1:
p_K(g,i) = \frac{\min(∑_{j \in K} p_w(g,i,j),
3)}{3}
Environmental aggregation: The environmental
pressures score p_E(g,i) is the weighted sum of
p_K(g,i), where each weight is the maximum weight
in the pressure category K, and then divided by the
sum of the maximum weights:
w_{K,max}(g,i) = max(\{\forall_j \in K |
w(g,i,j)\})
p_E(g,i) = \frac{∑_K w_{K,max}(g,i)
p_K(g,i)}{∑_K w_{K,max}(g,i)}
Social aggregation: The social pressures score
p_S(g,i) is the mean of the unweighted
social pressure scores p(i,j):
p_S(g,i) = \frac{∑_{j \in S}^{} p(i,j)}{N}
Gamma combination: The pressures score
p_X(g,i):
p_X(g,i) = γ p_E(g,i) + (1-γ)p_S(g,i)
Value
Returns a named vector with the pressures score for each
named region.
See Also
Examples
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
## Not run:
> ohi.pressure.category
$environmental
[1] "po" "hd" "fp" "sp" "cc"
$social
[1] "ss"
> p
pressure
region_id fp_art_hb fp_art_lb fp_com_hb fp_com_lb hd_intertidal
1 0.122 0.25 0.35 0.395 0.954
2 0.096 0.94 0.85 0.252 0.649
3 0.858 0.46 0.84 0.097 0.425
4 0.814 0.63 0.60 0.672 0.659
5 0.247 0.51 0.58 0.941 0.046
6 0.853 0.34 0.15 0.370 0.385
7 0.601 0.31 0.39 0.873 0.064
8 0.355 0.89 0.74 0.159 0.273
9 0.289 0.94 0.52 0.743 0.094
10 0.887 0.89 0.87 0.660 0.746
pressure
region_id hd_subtidal_hb hd_subtidal_sb po_chemicals po_nutrients
1 0.535 0.651 0.042 0.931
2 0.454 0.069 0.234 0.025
3 0.297 0.428 0.970 0.679
4 0.953 0.485 0.063 0.565
5 0.963 0.045 0.552 0.828
6 0.598 0.213 0.907 0.220
7 0.476 0.641 0.980 0.214
8 0.285 0.858 0.447 0.793
9 0.591 0.702 0.719 0.472
10 0.072 0.431 0.685 0.102
pressure
region_id sp_alien sp_genetic ss_wgi
1 0.979 0.761 0.181
2 0.345 0.091 0.631
3 0.223 0.986 0.646
4 0.035 0.078 0.559
5 0.992 0.643 0.432
6 0.963 0.416 0.221
7 0.752 0.627 0.257
8 0.100 0.245 0.333
9 0.316 0.373 0.347
10 0.283 0.224 0.031
> w
pressure
region_id fp_art_hb fp_art_lb fp_com_hb fp_com_lb hd_intertidal
1 2 1 0.92 1 1
2 2 1 0.48 1 1
3 2 1 2.81 1 1
4 2 1 1.19 1 1
5 2 1 2.82 1 1
6 2 1 1.07 1 1
7 2 1 1.48 1 1
8 2 1 0.46 1 1
9 2 1 0.56 1 1
10 2 1 0.90 1 1
pressure
region_id hd_subtidal_hb hd_subtidal_sb po_chemicals po_nutrients
1 2 2 1.00 1
2 2 2 0.79 1
3 2 2 0.37 1
4 2 2 0.91 1
5 2 2 1.06 1
6 2 2 0.72 1
7 2 2 0.49 1
8 2 2 1.18 1
9 2 2 0.18 1
10 2 2 0.28 1
pressure
region_id sp_alien sp_genetic ss_wgi
1 1 1 1
2 1 1 1
3 1 1 1
4 1 1 1
5 1 1 1
6 1 1 1
7 1 1 1
8 1 1 1
9 1 1 1
10 1 1 1
> p_x <- ohi.model.pressures(p, w)
> p_x
1 2 3 4 5 6 7 8 9 10
0.40 0.53 0.68 0.63 0.60 0.43 0.48 0.47 0.50 0.30
> data.frame(region_id=names(p_x), pressure=p_x)
region_id pressure
1 1 0.40
2 2 0.53
3 3 0.68
4 4 0.63
5 5 0.60
6 6 0.43
7 7 0.48
8 8 0.47
9 9 0.50
10 10 0.30
>
>
## End(Not run)
bbest/ohigui documentation built on May 11, 2019, 9:25 p.m.
R Package Documentation
Browse R Packages
We want your feedback!
Note that we can't provide technical support on individual packages. You should contact the package authors for that.
Description Usage Arguments Details Value See Also Examples
Description
The pressures model function computes a pressures score for each region given a weighting matrix for a goal and the individual pressures values.
Usage
1 | ohi.model.pressures(p, w, GAMMA = 0.5, browse = F)
|
Arguments
p |
the pressures value matrix [region_id x
pressure]. Each score must be a real number between 0 and
1 inclusive, or NA. The pressure names must be of the
form category_pressure where
category is one of the categories listed in
pressure region_id cc_acid cc_sst cc_uv fp_art_hb 1 0.879 0.360 0.764 NA 2 0.579 0.396 0.531 NA 3 0.926 0.235 0.769 NA 4 0.914 0.554 0.795 NA 5 0.860 0.609 0.802 0.001 6 0.871 0.325 0.788 0.001 7 0.846 0.410 0.677 0.000 8 0.806 0.671 0.752 NA 9 0.844 0.595 0.678 NA 10 0.860 0.575 0.781 0.109 |
w |
the weighting matrix of the form [region_id x pressure]. Each rank weight must be a real number between 0 and 3 inclusive, or NA. pressure region_id cc_acid cc_sst cc_uv fp_art_hb 1 2 1 0.6 NA 2 2 1 0.5 NA 3 2 1 2.1 NA 4 2 1 3.0 NA 5 2 1 2.8 1 6 2 1 2.2 1 7 2 1 1.3 1 8 2 1 1.7 NA 9 2 1 3.0 NA 10 2 1 1.2 1 |
GAMMA |
Multiplier used to combine environmental and social pressures. |
Details
Each pressure layer p(i,j) is either environmental or social, belongs to a pressures category K \in \{cc,fp,hd,po,sp,ss\}, and has a value (0..1) for each region i and pressures layer j. Each goal has a weight matrix w that has a rank weight between 0 and 3 inclusive, or NoData, for each region i and each pressure layer j on a per goal g basis.
The pressures scores calculations go through 5 steps, using a complex weighting scheme that varies across goals, subgoals, pressures categories, and regions:
-
g is the goal or subgoal (e.g., AO, CW, LIV, ECO, ...),
-
i is the region (e.g., 1, 2, 3, ...),
-
j is the pressures layer or stressor (e.g.,
cc_acid,fp_art_lb, etc.).
Calculations
Apply weights for each goal g, region i, and pressure layer j: Each weighted pressure p_w(g,i,j) is the pressure layer value p(i,j) per region i and pressure layer j multiplied by the rank weight w(g,i,j) for that goal g, region i, and pressure layer j. If the w(g,i,j) is NoData or 0, the weighted pressure p_w(g,i,j) is NoData.
p_{w}(g,i,j) = w(g,i,j) * p(i,j)
Category-level aggregation: The pressures category score p_K is the sum of all p_w within each category, then rescaled to 0..1 using a linear scale range transformation (from 0..3 to 0..1). Any score p_K greater than 1 is capped to 1:
p_K(g,i) = \frac{\min(∑_{j \in K} p_w(g,i,j), 3)}{3}
Environmental aggregation: The environmental pressures score p_E(g,i) is the weighted sum of p_K(g,i), where each weight is the maximum weight in the pressure category K, and then divided by the sum of the maximum weights:
w_{K,max}(g,i) = max(\{\forall_j \in K | w(g,i,j)\})
p_E(g,i) = \frac{∑_K w_{K,max}(g,i) p_K(g,i)}{∑_K w_{K,max}(g,i)}
Social aggregation: The social pressures score p_S(g,i) is the mean of the unweighted social pressure scores p(i,j):
p_S(g,i) = \frac{∑_{j \in S}^{} p(i,j)}{N}
Gamma combination: The pressures score p_X(g,i):
p_X(g,i) = γ p_E(g,i) + (1-γ)p_S(g,i)
Value
Returns a named vector with the pressures score for each named region.
See Also
Examples
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 | ## Not run:
> ohi.pressure.category
$environmental
[1] "po" "hd" "fp" "sp" "cc"
$social
[1] "ss"
> p
pressure
region_id fp_art_hb fp_art_lb fp_com_hb fp_com_lb hd_intertidal
1 0.122 0.25 0.35 0.395 0.954
2 0.096 0.94 0.85 0.252 0.649
3 0.858 0.46 0.84 0.097 0.425
4 0.814 0.63 0.60 0.672 0.659
5 0.247 0.51 0.58 0.941 0.046
6 0.853 0.34 0.15 0.370 0.385
7 0.601 0.31 0.39 0.873 0.064
8 0.355 0.89 0.74 0.159 0.273
9 0.289 0.94 0.52 0.743 0.094
10 0.887 0.89 0.87 0.660 0.746
pressure
region_id hd_subtidal_hb hd_subtidal_sb po_chemicals po_nutrients
1 0.535 0.651 0.042 0.931
2 0.454 0.069 0.234 0.025
3 0.297 0.428 0.970 0.679
4 0.953 0.485 0.063 0.565
5 0.963 0.045 0.552 0.828
6 0.598 0.213 0.907 0.220
7 0.476 0.641 0.980 0.214
8 0.285 0.858 0.447 0.793
9 0.591 0.702 0.719 0.472
10 0.072 0.431 0.685 0.102
pressure
region_id sp_alien sp_genetic ss_wgi
1 0.979 0.761 0.181
2 0.345 0.091 0.631
3 0.223 0.986 0.646
4 0.035 0.078 0.559
5 0.992 0.643 0.432
6 0.963 0.416 0.221
7 0.752 0.627 0.257
8 0.100 0.245 0.333
9 0.316 0.373 0.347
10 0.283 0.224 0.031
> w
pressure
region_id fp_art_hb fp_art_lb fp_com_hb fp_com_lb hd_intertidal
1 2 1 0.92 1 1
2 2 1 0.48 1 1
3 2 1 2.81 1 1
4 2 1 1.19 1 1
5 2 1 2.82 1 1
6 2 1 1.07 1 1
7 2 1 1.48 1 1
8 2 1 0.46 1 1
9 2 1 0.56 1 1
10 2 1 0.90 1 1
pressure
region_id hd_subtidal_hb hd_subtidal_sb po_chemicals po_nutrients
1 2 2 1.00 1
2 2 2 0.79 1
3 2 2 0.37 1
4 2 2 0.91 1
5 2 2 1.06 1
6 2 2 0.72 1
7 2 2 0.49 1
8 2 2 1.18 1
9 2 2 0.18 1
10 2 2 0.28 1
pressure
region_id sp_alien sp_genetic ss_wgi
1 1 1 1
2 1 1 1
3 1 1 1
4 1 1 1
5 1 1 1
6 1 1 1
7 1 1 1
8 1 1 1
9 1 1 1
10 1 1 1
> p_x <- ohi.model.pressures(p, w)
> p_x
1 2 3 4 5 6 7 8 9 10
0.40 0.53 0.68 0.63 0.60 0.43 0.48 0.47 0.50 0.30
> data.frame(region_id=names(p_x), pressure=p_x)
region_id pressure
1 1 0.40
2 2 0.53
3 3 0.68
4 4 0.63
5 5 0.60
6 6 0.43
7 7 0.48
8 8 0.47
9 9 0.50
10 10 0.30
>
>
## End(Not run)
|
R Package Documentation
Browse R Packages
We want your feedback!
Note that we can't provide technical support on individual packages. You should contact the package authors for that.
Embedding an R snippet on your website
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.