Decision theory is a formal theory of decision making under uncertainty.
A decision problem consists of:
We work to know What is the best action?
\[ a^*=arg_a max\{E[u(c) | a]\} \]
1.- Decision problem: Should patient be treated for disease?
2.- Decision model:
Consequences:
\( c_{DN} \) (disease for long time, no side effects)
| Disease | Well now | |
|---|---|---|
| Treat | \( c_{WS} \) | \( c_{WS} \) |
| No Treat | \( c_{DN} \) | \( c_{WN} \) |
Probabilities and Utilities:
u(\( c_{DN} \)) = 0 (no treatment - no side effects, patient had sickness, which he suffers long-time)
| Disease | Well now | |
|---|---|---|
| Treat | 90 | 90 |
| No Treat | 0 | 100 |
3.- Expected utility:
Best action is \( a_T \) (treat patient)
| Disease | Well now | Expected Utility | |
|---|---|---|---|
| Treat | 0.3 x 90 | 0.7 x 90 | 90 |
| No Treat | 0.3 X 0 | 0.7 x 100 | 70 |
4.- Sensitivity Analysis: Optimal Decision as Function of Sickness Probability
Expected utility of the two actions:
Typically we are uncertain about the value of p
ExpU_T = 90
FUN_ExpU_T <- function(p) (ExpU_T)
FUN_ExpU_N <- function(p) 100 * (1 - p)
DIFF_FUN <- function(x) FUN_ExpU_T(x) - FUN_ExpU_N(x)
p_vector = seq(from = 0, to = 1, length.out = 100)
plot(p_vector, FUN_ExpU_N(p_vector), type = "l", col = "red", axes = F)
axis(1, at = seq(0, 1, length.out = 11), labels = T)
axis(2)
abline(h = FUN_ExpU_T(p_vector), col = "blue")
p_intersection = uniroot(DIFF_FUN, c(0, 1))$root
abline(v = p_intersection, col = "darkgreen", lty = 2)
box()
The chart shows we should treat if p > 0.1,and don’t treat if p < 0.1
We may be able to perform a test before deciding whether to treat the patient:
We will assume:
How does the model change if test results are available?
Objective: use evidence E to update beliefs about a hypothesis H
Bayes as Odds: \begin{equation} \frac{P(H_1|E)}{P(H_2|E)}=\frac{\frac{P(E|H_1)P(H_1)}{P(E)}}{\frac{P(E|H_2)P(H_2)}{P(E)}}=\frac{P(E|H_1)P(H_1)}{P(E|H_2)P(H_2)}, ~~if ~P(E)>0, P(H2)>0 \end{equation}
Disease Example with Test
– If test is negative P(sD | tN) = (0.3 x 0.05)/(0.3 x 0.05 + 0.7 x 0.85) = 0.025 – If test is positive P(sD | tP) = (0.3 x 0.95)/(0.3 x 0.95 + 0.7 x 0.15) = 0.731 • Expected utility of treating patient is independent of test – EU(aT) = EU(aT | tN) = EU(aT | tP) = 90 • Expected utility of not treating patient depends on test – If test is negative EU(aN | tN) = 0.975 × 100 = 97.5 – If test is positive EU(aN | tP) = 0.269 × 100 = 26.9 – (Recall if no test EU(aN) = 0.7 × 100 = 70 ) • Optimal choice is strategy aF (FollowTest): – Treat if test is positive; don’t treat if test is negative