prob_abxprior, prob_StartBact_R are new parameters instead of prob_StartBact
model_simple.R
change N -> s
prob s = prob_abxprior(1-gen_StartBact_R)
prob S = (1-prob_abxprior)(1-gen_StartBact_R)
prob R = gen_StartBact_R
model_binary.R
prob of having had antibiotics prior to admission = 0.2 (prob_abxprior)
prob_start of S+s = 1 - gen_StartBact_R
prob_start of sR+Sr+sr = gen_Startbact_R
prob S = (1-prob_abxprior)(1-gen_StartBact_R)
prob s = prob_abxprior(1-gen_StartBact_R)
prob sR = (prob_abxprior/2)gen_StartBact_R
prob Sr = (1-prob_abxprior)gen_StartBact_R
prob sr = (prob_abxprior/2)*gen_StartBact_R
prob_abxprior, prob_StartBact_R are new parameters instead of prob_StartBact
model_simple.R
model_binary.R prob of having had antibiotics prior to admission = 0.2 (prob_abxprior) prob_start of S+s = 1 - gen_StartBact_R prob_start of sR+Sr+sr = gen_Startbact_R prob S = (1-prob_abxprior)(1-gen_StartBact_R) prob s = prob_abxprior(1-gen_StartBact_R) prob sR = (prob_abxprior/2)gen_StartBact_R prob Sr = (1-prob_abxprior)gen_StartBact_R prob sr = (prob_abxprior/2)*gen_StartBact_R
model_frequency.R same as in model_binary.R