(EC.1) — the net benefit of testing as a quadratic in the risk
ProvedFamilyGenTest.ec_1genetic-testingmarkov-decision-processoptimal-stoppingp2o-gran-per-chapterp2o-plan-paperp2o-strengthenedp2o-v1
Let be feasible and untested in . Then
This is the one-person net benefit of testing, the source of the two-sided form of the rule.
Differences from the source:
- reading: a set of genes in the model of §2; Assumption 4, Assumption 5 and Theorem 1's proof are single-gene → one gene: ; reason: Theorem 1 compares with two real thresholds, which needs a scalar; its proof, Assumption 4 and §4.3 are single-gene (pp.15, ec2), and no multi-gene threshold statement is printed
- generalization: the Mendelian prior (5) → an arbitrary mass function on ; reason: the argument does not use the pedigree and p.11 says priors incorporate all prior information; (5) is the special case given by the theorem eq_5_mendelian_prior
- encoding: a state as a partition of the tested individuals, state space →
State I := I → Option (Fin 3), withFeasiblemeaning ; reason: same objects when is the support of the prior; the family is aFintype, possibly empty - encoding: with Definition 3 and Assumption 1 →
testedReward: at result , at results ; reason: after a test the belief is the true carrier status, which is what Assumption 1 and Definition 3 give - generalization: the reward assumptions of §4 (Assumptions 1–4) → arbitrary real rewards and costs; reason: the identity is algebraic in the data and does not use them; every instance of the source is covered
Preamble
import Mathlib import Definitions.Def_FamilyGenTest_Rewards
Formal statement
namespace FamilyGenTest
variable {I : Type*} [Fintype I] [DecidableEq I]
/-- (EC.1) (Adelman & Wang, Optimal Genetic Testing of Families, SSRN 4624694 (Jan 24 2023), p.ec2):
`Eᵢ[rᵢ(ℙᵢ(S′)) | s] − cᵢ(s) − r̄ᵢ(ℙᵢ(s)) = −rᵢ^δ ℙᵢ(s)² + (rᵢ^δ + dᵢ − cᵢᶠ) ℙᵢ(s) − cᵢˢ`. -/
theorem ec_1 (R : Rewards I) (μ : Prior I) (s : State I) (hs : Feasible μ s) (i : I)
(hi : s i = none) :
(∑ γ, marginal μ s i γ * testedReward R i γ) - cost R μ s i - rbar R i (risk μ s i)
= -rdelta R i * risk μ s i ^ 2 + (rdelta R i + censoring R i - R.cf i) * risk μ s i
- R.cs i := by sorry
end FamilyGenTest
Source
Adelman & Wang, Optimal Genetic Testing of Families, SSRN 4624694 (Jan 24 2023), (EC.1), p.ec2
Human review
Confirmed by the mission captain (proposal self-audit).
Confirmed by the moderator at approval.