-
Notifications
You must be signed in to change notification settings - Fork 0
/
Copy pathSA_parameters.sas
119 lines (69 loc) · 4 KB
/
SA_parameters.sas
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
* These parameters are up to date at 17th Oct 2023;
* caldate1, caldate_never_dot and startyr specified here to overwrite core value;
* country added to switch on SA-spcific abort statements;
* Note that rate_anc_inc has been updated from rate_testanc_inc (used in some other SA analyses) to be consistent with the current core program;
* JAS Oct23;
****% on ART of those diag too high***
/***************************************************************************/
****************** SOUTH AFRICA SPECIFIC PARAMETERS ******************
/***************************************************************************/
* country; country='South Africa';
* caldate1; caldate1 = 1984; *core: 1989;
* caldate_never_dot; caldate_never_dot=1984; *core: 1989;
* POPULATION GROWTH AND DEMOGRAPHY;
* inc_cat; inc_cat = 3;
* population aged >= 15 in 1990: 22 million ;
* fold_change_ac_death_rate;fold_change_ac_death_rate_w = 0.65; fold_change_ac_death_rate_m = 0.65;
* INTRODUCTION OF HIV;
* startyr; startyr = 1993;
* newp_seed; newp_seed = 10 ;
* SEXUAL BEHAVIOUR;
* ych_risk_beh_newp; ych_risk_beh_newp = 1;
* ych_risk_beh_ep; ch_risk_beh_ep = 1;
* ych2_risk_beh_newp; %sample(ych2_risk_beh_newp,
0.95 0.99 1 1/0.99 1/0.95,
0.15 0.15 0.7 0 0 );
* p_rred_p; %sample(p_rred_p, 0.1 0.2 0.3, 0.3 0.4 0.3);
* p_hsb_p; %sample(p_hsb_p, 0.05 0.10 0.15 0.20, 0.1 0.3 0.5 0.1 );
* sex_age_mixing_matrix_m; %sample(sex_age_mixing_matrix_m, 1 2 3 4 5 6 , 0.1 0.1 0.1 0.1 0.1 0.5);
* sex_age_mixing_matrix_w; %sample(sex_age_mixing_matrix_w, 1 2 3 4 5 6 , 0.3 0.3 0.1 0.1 0.1 0.1);
* TRANSMISSION;
* fold_tr; fold_tr = 1;
* HIV TESTING;
* an_lin_incr_test; an_lin_incr_test = 0.017 ;
* rate_anc_inc; rate_anc_inc =0.05;
* date_test_rate_plateau; date_test_rate_plateau= 2020.5;
* NATURAL PROGRESSION AND RISK OF HIV RELATED CONDITIONS;
* fold_change_in_risk_base_rate;
* fold_change_in_risk_base_rate = 0.85;
* fold_decr_hivdeath; fold_decr_hivdeath = 0.15 ;
* fx; fx = 0.85;
* gx; * gx = 1;
* tb_base_prob_diag_l; * %sample_uniform(tb_base_prob_diag_l, 0.50 0.75); * base probability that tb is diagnosed late ;
* crypm_base_prob_diag_l; * %sample_uniform(crypm_base_prob_diag_l, 0.50 0.75); * base probability that crypm is diagnosed late ;
* sbi_base_prob_diag_l; * %sample_uniform(sbi_base_prob_diag_l, 0.50 0.75); * base probability that sbi is diagnosed late ;
* oth_adc_base_prob_diag_l; * %sample_uniform(oth_adc_base_prob_diag_l, 0.50 0.75); * base probability that other adc is diagnosed late;
* incr_death_rate_oth_adc; incr_death_rate_oth_adc = 1.5 ; * note values linked to adjustments to fold_decr_hivdeath ;
* incr_death_rate_crypm; incr_death_rate_crypm = 5 ;
* incr_death_rate_sbi; incr_death_rate_sbi = 1.5 ;
* incr_death_rate_tb; incr_death_rate_tb = 1.5 ;
* LINKAGE, RETENTION, MONITORING, LOSS, RETURN, INTERRUPTION OF ART AND RESTARTING, ART;
* prob_loss_at_diag; prob_loss_at_diag = 0.4;
* rate_return; rate_return = 0.2;
* prob_return_adc; prob_return_adc = 0.95 ;
* art_intro_date; art_intro_date = 2005.5;
* pr_art_init; pr_art_init = 0.95;
* rate_int_choice; rate_int_choice = 0.03 ;
* prob_vl_meas_done; prob_vl_meas_done = 0.9;
* rate_ch_art_init_str;
rate_ch_art_init_str_4 = 1.0; rate_ch_art_init_str_9 = 0.05; rate_ch_art_init_str_10 = 0.1 ; rate_ch_art_init_str_3 = 0.1;
* adh_pattern; adh_pattern = 3;
* ind_effect_art_hiv_disease_death; ind_effect_art_hiv_disease_death =0.3 ;
* rr_non_aids_death_hiv_off_art; rr_non_aids_death_hiv_off_art = 2;
* rr_non_aids_death_hiv_on_art; rr_non_aids_death_hiv_on_art = 1.3;
* SEX WORKERS;
* base_rate_sw; %sample(base_rate_sw, 0.0007 0.0010, 0.45 0.55);
* CIRCUMCISION;
* prob_birth_circ; prob_birth_circ=0.39;
* rate_test_startprep_any; %sample_uniform(rate_test_startprep_any, 0.10 0.15 0.30);
* pref_prep_oral_beta_s1; %sample_uniform(pref_prep_oral_beta_s1, 0.6 0.7 0.8 0.9 1.0 1.1) ;