-
Notifications
You must be signed in to change notification settings - Fork 0
/
Copy pathEVPI_Poster.Rmd
101 lines (76 loc) · 5.69 KB
/
EVPI_Poster.Rmd
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
---
title: Expected Value of Perfect Information
csl: nature.csl
author:
- name: Katja Schiffers
affil: 1
orcid: '0000-0003-2559-6579'
- name: Cory Whitney
affil: 1
orcid: '0000-0003-4988-4583'
- name: Eike Luedeling
affil: 1
orcid: '0000-0002-7316-3631'
affiliation:
- num: 1
address: INRES Horticultural Sciences, University of Bonn
column_numbers: 2
logoright_name: images/Uni_Bonn_Picture.png
logoleft_name: images/INRES_Logo.png
primary_colour: "white"
title_textcol: "black"
author_textcol: "black"
affiliation_textcol: "black"
sectitle_textcol: "black"
poster_height: "53in"
poster_width: "33in"
output:
posterdown::posterdown_html:
self_contained: false
bibliography: EVPI.bib
---
```{r setup, include=FALSE}
knitr::opts_chunk$set(echo = FALSE)
```
# What is the EVPI?
The expected value of perfect information (EVPI) is a concept developed in decision analysis. It quantifies the potential gain in units of an output variable of a projection model given accurate knowledge on the value a tested input variable will take. EVPI can help identify those variables for which more information would maximize the model output [@lanzanova2019]. When the output variable is monetary, the EVPI can be interpreted as the highest price that decision makers should be willing to pay for perfect information [@hubbard2014].
EVPI is calculated as the difference between the expected value of the outcome variable $EV$ given accurate knowledge on the value the tested input variable will take (perfect information - $PI$) and the expected maximum value ($EMV$) of the outcome variable given only knowledge about the probability distribution of the input variable (imperfect knowledge $II$) [@Strong2014].
$$
EVPI = EV|PI - EMV
$$
Note that the imperfect information is defined as exact knowledge on the probability distribution of the input variable. It is not meant to capture knowledge uncertainty about the true value of the variable.
# Examples
## Discrete probability distribution
Calculating $EVPI$ with discrete probability distributions is relatively straightforward. We use the example of investment in the stock market vs. deposits (bonds):
<img src="images/stock_market_example.jpg" alt="stock market example" align="center" width = "75%" height="75%"/>
The expected value of deposit investment $Expdeposit$ is calculated as the probability of different states of the economy (x-axis) times the expected loss or gain in each condition (y-axis):
$$
Expdeposit : 0.2 \cdot 500 + 0.3 \cdot 500
+ 0.5 \cdot 500 = 500
$$
Likewise, expected value of stock investment $Expstock$ is the probability of each possible state of the economy (x-axis) times the expected losses or gains (y-axis):
$$ Expstock: 0.2 \cdot -800 + 0.3 \cdot 600
+ 0.5 \cdot 1500 = 770 $$
Because $Expstock$ is the more likely decision to result in a gain it is referred to as the Expected Maximum Value $EMV$.
The expected value of the decision *given* perfect information $EV|PI$ prior to decision is calculated as the sum of all the best options (always making the more gainful choice) multipled by the respective probabilities:
$$
EV|PI : 0.2 \cdot 500 + 0.3 \cdot 600
+ 0.5 \cdot 1500 = 1030
$$
The Expected Value of Perfect Information $EVPI$ is then calculated as the difference between the decision given perfect information $EV|PI$ and the the Expected Maximum Value $EMV$.
$$
EVPI : 1030 – 770 = 260
$$
Knowing the direction the market will go (having perfect information) before making our decision would help us take the best deision here. We should be willing to pay up to 350€ for perfect information on the future state of the economy.
## Continuous probability distribution
In the following example of an agricultural decision to apply fertilizer, we assume that the nutrient content of the soil will follow a normal distribution and that the additional gain in yield under the decision not to apply fertilizer will remain constant (<span style="color:green">green line</span>). Therefore the expected gain for the decision to not apply fertilizer is simply zero.
Further we assume when applying fertilizer the effect on additional gain is negatively correlated with the nutrient content of the soil and becomes even negative at a certain level of nutrients and causes losses in yield due to toxic effects (<span style="color:red">dashed red line</span>).
To calculate the expected gain in yield when applying fertilizer, we multiply the dashed red line with the normal density distribution (analog to multiplying the gains and losses with the probability of the marked going in the respective directions in the example above). For the resulting weighted response (<span style="color:red">red line</span>) we add up the areas under the curve: $A-B$. As the positive part $A$ is greater than the negative part $B$, the expected value of applying fertilizer has a small positive value. Therefore, applying fertilizer is the decision with the expected maximum value $EMV$ under imperfect knowledge.
<img src="images/EVPI_figure2.png" alt="agriculture example" align="center" width = "75%" height="75%"/>
The expected value given perfect information $EV|PI$ is equal to $A$, because up to the nutrient content, where the line intersects zero, the best decision is to apply fertilizer, but for higher nutrient contents no fertilizer should be added. The $EVPI$ is then again caculated as the difference between $EV|PI$ and $EMV$.
$$
EVPI : A – (A - B) = B
$$
# R Tools
The *empirical_EVPI()* function in R's decisionSupport library [@decisionSupport2017] calculates EVPI for a simple model with continuous data like the one above. The *multi_EVPI()* function does the same with more complex models with multiple variables.
# References