-
Notifications
You must be signed in to change notification settings - Fork 9
/
mr.sthlp
85 lines (63 loc) · 2.57 KB
/
mr.sthlp
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
{smcl}
{* *! version 0.1.0 19jan2020 Tom Palmer}{...}
{vieweralsosee "mrrobust" "help mrrobust"}{...}
{viewerjumpto "Syntax" "mr##syntax"}{...}
{viewerjumpto "Description" "mr##description"}{...}
{viewerjumpto "Options" "mr##options"}{...}
{viewerjumpto "Examples" "mr##examples"}{...}
{viewerjumpto "Stored results" "mr##results"}{...}
{viewerjumpto "Author" "mr##author"}{...}
{title:Title}
{p 5}
{bf:mr} {hline 2} Mendelian randomization programs
{p_end}
{marker syntax}{...}
{title:Syntax}
{p 8 16 2}
{opt mr} {it:subcommand} ... [{it:aweight}] {ifin} [{cmd:,} {it:options}]
{marker description}{...}
{title:Description}
{pstd}
{cmd:mr} is a simple wrapper to the commands in the {help mrrobust} package.
{pstd}
The {it:subcommand} is specified as the mrrobust program name without its mr prefix, i.e. {cmd:mregger ...} can alternatively be run using the syntax {cmd:mr egger ...}.
{marker options}{...}
{title:Options}
{phang}
{cmd:mr} takes the options for the specified subcommand.
{marker examples}{...}
{title:Examples}
{pstd}Using the data provided by {help mr##do:Do et al. (2013)} recreate
{help mr##bowden:Bowden et al. (2016)}, Table 4, LDL-c "All genetic variants" estimates.{p_end}
{pstd}Setup{p_end}
{phang2}{cmd:.} {stata "use https://raw.github.com/remlapmot/mrrobust/master/dodata, clear"}{p_end}
{pstd}Select observations ({it:p}-value with exposure < 10^-8){p_end}
{phang2}{cmd:.} {stata "gen byte sel1 = (ldlcp2 < 1e-8)"}{p_end}
{pstd}IVW (with fixed effect standard errors){p_end}
{phang2}{cmd:.} {stata "mr egger chdbeta ldlcbeta [aw=1/(chdse^2)] if sel1==1, ivw fe"}{p_end}
{pstd}Scatter plot of MR-Egger model{p_end}
{phang2}{cmd:.} {stata "mr eggerplot chdbeta chdse ldlcbeta ldlcse if sel1==1"}{p_end}
{pstd}Multivariable IVW{p_end}
{phang2}{cmd:.} {stata "mr mvivw chdbeta ldlcbeta hdlcbeta tgbeta [aw=1/(chdse^2)] if sel1==1"}{p_end}
{marker results}{...}
{title:Stored results}
{pstd}
{cmd:mr} returns the results from the specified subcommand.
{marker references}
{title:References}
{marker bowden}{...}
{phang}
Bowden J, Davey Smith G, Haycock PC, Burgess S. Consistent estimation in
Mendelian randomization with some invalid instruments using a weighted median
estimator. Genetic Epidemiology, 2016, 40, 4, 304-314.
{browse "https://dx.doi.org/10.1002/gepi.21965":DOI}
{p_end}
{marker do}{...}
{phang}
Do et al., 2013. Common variants associated with plasma triglycerides and risk
for coronary artery disease. Nature Genetics. 45, 1345-1352.
{browse "https://dx.doi.org/10.1038/ng.2795":DOI}
{p_end}
{marker author}{...}
{title:Author}
INCLUDE help mrrobust-author