Resistance mutations

[dimitrasal]

Resistance mutations need to be accounted for when calculating probability of infection from infected long term partner
Variables on total proportion of people with resistance mutations t_prop_rm and for viral load groups t_prop_vlg1...6_rm are used as dummy variables and need to be estimated.

if epi=1 then do;  * dependent_on_time_step_length ;  
	if epvls=1    then do; risk_eip = max(0,tr_rate_undetec_vl+(0.000025*rand('normal')));          vl_source=1;	t_prop_rm=t_prop_vlg1_rm; end;
	if epvls ne 1 then do;
		if epi_tm1 =0 then do; risk_eip = max(0,tr_rate_primary+(0.075*rand('normal'))); ep_primary=1; vl_source=6;	t_prop_rm=t_prop_vlg6_rm; end;* infection in primary;
		if epi_tm1 =1 then do; risk_eip = max(0,(0.05*fold_tr)+(0.0125*rand('normal')));               vl_source=4;	t_prop_rm=t_prop_vlg4_rm; end;