Watanabe2018_State Transition Model with Treatment and Costs

Model with functions depending on Age, Male, BP (Blood Pressure).
There are 3 disease states: Healthy, Sick, and Dead, where the Dead state is terminal. The yearly transition probabilities are: Healthy to Dead: Age/1000; Healthy to Sick: According to function F1 depending on Age and Male and BP; Sick to Healthy: 0.1; Sick to Dead: according to function F2 depending on Age and Male.
Pre-Transition Rules: Age increased by 1 and BP by Age/10 each simulation cycle. Post-Transition Rules: Treatment = BP>140 , becomes 1 when BP crosses 140 threshold; BP =BP-Treatment*10 , meaning a drop of 10 once treatment is applied; CostThisYear = Age + \Treatment*10 , cost depends on age and if treatment was taken; Cost= Cost + CostThisYear , it accumulates cost over time.
Initial conditions: Healthy = (50 Male, 50 Female with Age =1,2,...,50 for each individual), BP =120, Sick = (0,0) and Dead = (0,0).
Output: Number of men and women in each disease state for years 1-10 and their ages and costs in each state. A stratified report by male and female and young – up to age 30 and old above age 30 is produced.