siman/siman_examples.do at master · UCL/siman · GitHub

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
/*
siman demonstration for paper
IW 24jul2025
*/

// OPEN DATA AND INTRODUCE SOME ERRORS

use simcheck, clear
gen estimand = "logOR"
* changes to insert some errors
* dgms are repeated, for dgm=MNAR
replace  b =  b[_n-300] if rep>100 & dgm=="MNAR"
replace se = se[_n-300] if rep>100 & dgm=="MNAR"
replace df = df[_n-300] if rep>100 & dgm=="MNAR"
* 1 in 3 of SEs are too small, for method=MI
replace se = se/2 if method==3 & mod(rep,3)==0

// RUN SIMAN SETUP
siman setup, rep(rep) dgm(dgm) method(method) target(estimand) est(b) se(se) df(df) true(0)

// GRAPHS OF ESTIMATES

* siman swarm
siman swarm, row(1) xsize(5) ysize(3) meangr(mcol(black))

* siman scatter
siman scatter, xsize(5) ysize(3)

* siman cms
siman comparemethodsscatter if dgm == 3, xsize(4) ysize(3)

* siman blandaltman
siman blandaltman if dgm == 3, xsize(5) ysize(3)

* siman zipplot
siman zipplot, xsize(5) ysize(5)

// COMPUTE PERFORMANCE ESTIMATES

use simcheck, clear
qui siman setup, rep(rep) dgm(dgm) method(method) est(b) se(se) df(df) true(0)
siman analyse
siman table bias empse, row(_perf dgm) nformat(%6.3f)

// GRAPHS OF PERFORMANCE

* siman lollyplot
siman lollyplot, legend(row(1)) labformat(%6.3f %6.0f) xsize(5) ysize(3)

* siman nestloop
use extendedtestdata, clear
siman setup, rep(rep) dgm(beta pmiss mech) method(method) target(estimand) est(b) se(se) true(true)
siman analyse
siman nestloop empse if estimand=="effect", stagger(.05) lcol(black = =) lwidth(1 = =) xsize(5) ysize(3)