Expected value of sample information

Formulation

Let

d ∈ D the decision being made, chosen from space  D x ∈ X an uncertain state, with true value in space  X z ∈ Z an observed sample composed of  n  observations  ⟨ z 1 , z 2 , . . , z n ⟩ U ( d , x ) the utility of selecting decision  d  from  x p ( x ) the prior subjective probability distribution (density function) on  x p ( z | x ) the conditional prior probability of observing the sample  z

It is common (but not essential) in EVSI scenarios for Z i = X , p ( z | x ) = ∏ p ( z i | x ) and ∫ z p ( z | x ) d z = x , which is to say that each observation is an unbiased sensor reading of the underlying state x , with each sensor reading being independent and identically distributed.

The utility from the optimal decision based only on the prior, without making any further observations, is given by

E [ U ] = max d ∈ D   ∫ X U ( d , x ) p ( x )   d x .

If the decision-maker could gain access to a single sample, z , the optimal posterior utility would be

E [ U | z ] = max d ∈ D   ∫ X U ( d , x ) p ( x | z )   d x

where p ( x | z ) is obtained from Bayes' rule:

p ( x | z ) = p ( z | x ) p ( x ) p ( z ) ; p ( z ) = ∫ p ( z | x ) p ( x )   d x .

Since they don't know what sample would actually be obtained if one were obtained, they must average over all possible samples to obtain the expected utility given a sample:

E [ U | S I ] = ∫ Z E [ U | z ] p ( z ) d z = ∫ Z max d ∈ D   ∫ X U ( d , x ) p ( z | x ) p ( x )   d x   d z .

The expected value of sample information is then defined as

E V S I = E [ U | S I ] − E [ U ] = ( ∫ Z max d ∈ D   ∫ X U ( d , x ) p ( z | x ) p ( x )   d x   d z ) − ( max d ∈ D   ∫ X U ( d , x ) p ( x )   d x ) .

Computation

It is seldom feasible to carry out the integration over the space of possible observations in E[U|SI] analytically, so the computation of EVSI usually requires a Monte Carlo simulation. The method involves randomly simulating a sample, z i = ⟨ z 1 i , z 2 i , . . , z n i ⟩ , then using it to compute the posterior p ( x | z i ) and maximizing utility based on p ( x | z i ) . This whole process is then repeated many times, for i = 1 , . . , M to obtain a Monte Carlo sample if optimal utilities. These are averaged to obtain the expected utility given a hypothetical sample.

Example

A regulatory agency is to decide whether to approve a new treatment. Before making the final approve/reject decision, they ask what the value would be of conducting a further trial study on n subjects. This question is answered by the EVSI.

The diagram shows an influence diagram depiction of an Analytica model for computing the EVSI in this example. For the reader who wishes to study this computation in greater detail, the model can be viewed and evaluated from Analytica Web Player.

The model classifies the outcome for any given subject into one of five categories:

Z i = {"Cure", "Improvement", "Ineffective", "Mild side-effect", "Serious side-effect"}

And for each of these outcomes, assigns a utility equal to an estimated patient-equivalent monetary value of the outcome.

A decision state, x in this example is a vector of five numbers between 0 and 1 that sum to 1, giving the proportion of future patients that will experience each of the five possible outcomes. For example, a state x = [ 5 % , 60 % , 20 % , 10 % , 5 % ] denotes the case where 5% of patients are cured, 60% improve, 20% find the treatment ineffective, 10% experience mild side-effects and 5% experience dangerous side-effects.

The prior, p ( x ) is encoded using a Dirichlet distribution, requiring five numbers (that don't sum to 1) whose relative values capture the expected relative proportion of each outcome, and whose sum encodes the strength of this prior belief. In the diagram, the parameters of the Dirichlet distribution are contained in the variable dirichlet alpha prior, while the prior distribution itself is in the chance variable Prior. The probability density graph of the marginals is shown here:

In the chance variable Trial data, trial data is simulated as a Monte Carlo sample from a Multinomial distribution. For example, when Trial_size=100, each Monte Carlo sample of Trial_data contains a vector that sums to 100 showing the number of subjects in the simulated study that experienced each of the five possible outcomes. The following result table depicts the first 8 simulated trial outcomes:

Combining this trial data with a Dirichlet prior requires only adding the outcome frequencies to the Dirichlet prior alpha values, resulting in a Diriclet posterior distribution for each simulated trial. For each of these, the decision to approve is made based on whether the mean utility is positive, and using a utility of zero when the treatment is not approved, the Pre-posterior utility is obtained. Repeating the computation for a range of possible trial sizes, an EVSI is obtained at each possible candidate trial size as depicted in this graph:

Historical background

The 1961 book Applied Statistical Decision Theory by Schlaifer and Raiffa was among the earliest to utilize EVSI extensively.

Comparison to related measures

Expected value of sample information (EVSI) is a relaxation of the expected value of perfect information (EVPI) metric, which encodes the increase of utility that would be obtained if one were to learn the true underlying state, x . Essentially EVPI indicates the value of perfect information, while EVSI indicates the value of some limited and incomplete information.

The expected value of including uncertainty (EVIU) compares the value of modeling uncertain information as compared to modeling a situation without taking uncertainty into account. Since the impact of uncertainty on computed results is often analysed using Monte Carlo methods, EVIU appears to be very similar to the value of carrying out an analysis using a Monte Carlo sample, which closely resembles in statement the notion captured with EVSI. However, EVSI and EVIU are quite distinct—a notable difference between the manner in which EVSI uses Bayesian updating to incorporate the simulated sample.