  Mendelian randomization (MR) is an instrumental variable (IV) approach to
infer causal relationships between exposures and outcomes with genome-wide
association studies (GWAS) summary data. However, the multivariable
inverse-variance weighting (IVW) approach, which serves as the foundation for
most MR approaches, cannot yield unbiased causal effect estimates in the
presence of many weak IVs. To address this problem, we proposed the MR using
Bias-corrected Estimating Equation (MRBEE) that can infer unbiased causal
relationships with many weak IVs and account for horizontal pleiotropy
simultaneously. While the practical significance of MRBEE was demonstrated in
our parallel work (Lorincz-Comi (2023)), this paper established the statistical
theories of multivariable IVW and MRBEE with many weak IVs. First, we showed
that the bias of the multivariable IVW estimate is caused by the
error-in-variable bias, whose scale and direction are inflated and influenced
by weak instrument bias and sample overlaps of exposures and outcome GWAS
cohorts, respectively. Second, we investigated the asymptotic properties of
multivariable IVW and MRBEE, showing that MRBEE outperforms multivariable IVW
regarding unbiasedness of causal effect estimation and asymptotic validity of
causal inference. Finally, we applied MRBEE to examine myopia and revealed that
education and outdoor activity are causal to myopia whereas indoor activity is
not.
