The present invention discloses a Mendelian
randomization causal inference method based on the summary statistical data of whole
genome association study, comprising: S1, integrating the
inference model, the polygenic effect model and the
estimation error model, adjusting the LD effect, and obtaining the MR-MU model; S2, collecting the GWAS data of the studied
phenotype and its
risk factor candidate set; S3, screening the
risk factor candidate set; S4, excluding the characteristics showing weak pleiotropy evidence; S5, selecting appropriate SNPs from the GWAS data as instrumental variables; S6, performing parameter
estimation on the conditional likelihood function of the MR-MU model, respectively obtaining two sets of parameter
estimation values under the conditions of fixed and unfixed causal effects, and the estimation value of the causal effect under unfixed conditions, and calculating the function values of the conditional likelihood function corresponding to the two sets of parameter estimation values; S7, calculating the
likelihood ratio test statistic, identifying the risk factors of the studied
phenotype, and statistically inferring whether there is a causal effect, so as to improve the accuracy and reliability of the
causal inference.