Skip to content

Commit ab59419

Browse files
add V&V plan for model 2; edits to V&V plan and summary for model 1
1 parent a0e25a5 commit ab59419

1 file changed

Lines changed: 44 additions & 7 deletions

File tree

docs/source/models/concept_models/vivarium_alzheimers/concept_model.rst

Lines changed: 44 additions & 7 deletions
Original file line numberDiff line numberDiff line change
@@ -383,7 +383,11 @@ scenario, and input draw.
383383
with 100K simulants each)
384384
- None
385385
* - 1.0
386-
- * Verify crude birth rate (CBR) against GBD
386+
- **Note:** All these checks can be done separately for each age
387+
group and sex, but due to the large number of age groups, it may
388+
be more prudent to start by looking at aggregated results.
389+
390+
* Verify crude birth rate (CBR) against GBD
387391
* Verify ACMR against GBD
388392
* Validate Alzheimer's CSMR against GBD
389393
* Verify Alzheimer's incidence rate against GBD
@@ -393,12 +397,45 @@ scenario, and input draw.
393397
* Check whether overall population remains stable over time
394398
* Check whether Alzheimer's prevalence remains stable over time
395399
* For comparison with model 2, calculate total "real world"
396-
Alzheimer's population over time as :math:`p_\text{AD} \cdot X_t
397-
/ S`, where :math:`p_\text{AD}` is prevalence of AD, :math:`X_t`
398-
is the simulated population at time :math:`t`, and :math:`S =
399-
X_{2025}` / (real population in 2025) is the model scale.
400+
Alzheimer's population over time as :math:`p_\text{AD}(t) \cdot
401+
X_t / S`, where :math:`p_\text{AD}(t)` is prevalence of AD at
402+
time :math:`t`, :math:`X_t` is the simulated population at time
403+
:math:`t`, and :math:`S = X_{2025}` / (real total population in
404+
2025) is the model scale
400405
- * Birth observer was missing, so we couldn't verify CBR
406+
* Total population per draw was 200k instead of 100k, and there
407+
were 10 draws instead of 25
401408
* Timestep was 182 days, resulting in 3 timesteps in 2025, making
402-
population counts 1.5 times what they should be; we'll change
403-
timestep to 183 days for future models
409+
population counts 1.5 times what they should be in 2025; we'll
410+
change the timestep to 183 days for future models
411+
* Total population decreased monotonically during the 76 years of
412+
the sim from 200k to about 170k in USA and about 125k in China
413+
-
414+
* - 2.0
415+
- **Note:** All these checks can be done separately for each age
416+
group and sex, but it may be more prudent to start by looking at
417+
aggregated results.
418+
419+
* Verify the number of new simulants per year against the :ref:`AD
420+
population model <other_models_alzheimers_population>`
421+
* Use interactive sim to verify initial population structure
422+
against the :ref:`AD population model
423+
<other_models_alzheimers_population>`
424+
* Verify that all simulants in the model have AD (i.e., all
425+
recorded person-time is in the "AD" state, not the "susceptible"
426+
state)
427+
* Verify that there are no transitions between AD states during
428+
the simulation (since it's an SI model and all simulants should
429+
be in the I state the whole time)
430+
* Verify ACMR against GBD
431+
* Validate Alzheimer's CSMR against GBD
432+
* Validate Alzheimer's EMR against GBD
433+
* Validate Alzheimer's YLLs and YLDs against GBD
434+
* For comparison with model 1, calculate total "real world"
435+
Alzheimer's population over time as :math:`X_t / S`, where
436+
:math:`X_t` is the simulated population at time :math:`t`, and
437+
:math:`S = X_{2025}` / (real population with AD in 2025) is the
438+
model scale (I'm not sure how closely we expect this to match
439+
model 1)
440+
-
404441
-

0 commit comments

Comments
 (0)