* Example for the SJ paper 
* Deterministic and probabilistic sensitivity analysis of epidemiological results

episensi 45 94 257 945 , st(cc) dseca(c(.9)) dspca(c(.9)) dsenc(c(.9)) dspnc(c(.9)) 

episensi 45 94 257 945, st(cc) dseca(c(.9)) dspca(c(.8)) dsenc(c(.8)) dspnc(c(.8))

episensi  45 94 257 945 , st(cc)  dpscex(c(.7)) dpscun(c(.7))  dpsnex(c(.6)) dpsnun(c(.6))

episensi  45 94 257 945 , st(cc)  dpscex(c(.9)) dpscun(c(.5))  dpsnex(c(.5)) dpsnun(c(.7))

episensi  45 94 257 945 , st(cc)  dpscex(c(.5)) dpscun(c(.9))  dpsnex(c(.7)) dpsnun(c(.5))

episensi 45 94 257 945, dpexp(c(.7)) dpunexp(c(.5)) drrcd(c(5))

episensi  45 94 257 945 , st(cc) reps(20000)  nodots    ///
dseca(trap(.75  .85 .95  1) )   dspca(trap(.75  .85 .95  1) )  ///
dsenc(trap(.75  .85 .95  1) )   dspnc(trap(.75  .85 .95  1))   ///
seed(123)

episensi  45 94 257 945 , st(cc) reps(20000)  nodots    ///
dseca(trap(.75  .85 .95  1) )   dspca(trap(.75  .85 .95  1) )  ///
dsenc(trap(.75  .85 .95  1) )   dspnc(trap(.75  .85 .95  1))   ///
seed(123)

episensi  45 94 257 945  , st(cc) reps(20000) nodots ///
 dseca(trap(.75  .85 .95  1) )  dspca(trap(.75  .85 .95  1) )  ///
 dsenc(trap(.7  .8 .9  .95) )   dspnc(trap(.7  .8 .9  .95) )   ///  
 corrsens(.8)  corrspec(.8)   seed(123)


episensi  45 94 257 945 , st(cc) reps(20000) nodots ///
  dsbfactor(log-n(0 0.21))  grprior  seed(123)

episensi  45 94 257 945 , st(cc) reps(20000) nodots ///
dpexp(uni(.4 .7)) dpunexp(uni(.4 .7)) drrcd(log-n(2.159 .280))  seed(123) grarrtot grprior 

episensi  45 94 257 945 , st(cc) reps(20000) nodots ///
dpunexp(uni(.4 .7)) dorce(log-n(0 0.639)) drrcd(log-n(2.159 .280))  seed(123) grarrtot grprior 


episensi  45 94 257 945 , st(cc) reps(20000) nodots    ///
dseca(trap(.75  .85 .95  1) )  dspca(trap(.75  .85 .95  1) )  ///
dsenc(trap(.7  .8 .9  .95) )  dspnc(trap(.7  .8 .9  .95))  ///
corrsens(.8) corrspec(.8) ///
dsbfactor(log-n(0 0.21)) ///
dpunexp(uni(.4 .7)) dorce(log-n(0 0.639)) drrcd(log-n(2.159 .280)) ///
seed(123) combined