Hi,
It would be great if you can help me with the following:
mr_results <- mr(dat = data, method_list = c("mr_two_sample_ml","mr_egger_regression", "mr_egger_regression_bootstrap","mr_simple_median","mr_weighted_median","mr_ivw","mr_ivw_radial","mr_ivw_mre","mr_ivw_fe","mr_simple_mode","mr_weighted_mode","mr_weighted_mode_nome"))
When I select the following methods on the website/shiny app MRBase , everything works properly
But, when i try to do this in a r script I'm getting the following error:
Analysing 'UKB-a:230' on '7'
Radial IVW
Estimate Std.Error t value Pr(>|t|)
Effect (Mod.2nd) 0.05760315 0.02928647 1.966886 0.049196349
Iterative 0.05760315 0.02928647 1.966886 0.049196349
Exact (FE) 0.06030600 0.02211508 2.726917 0.006392912
Exact (RE) 0.05910238 0.02987084 1.978598 0.052237094
Residual standard error: 1.324 on 63 degrees of freedom
F-statistic: 3.87 on 1 and 63 DF, p-value: 0.0536
Q-Statistic for heterogeneity: 110.5133 on 63 DF , p-value: 0.0002031823
Outliers detected
Number of iterations = 2
Error in pnorm(abs(b/se), low = FALSE) :
Non-numeric argument to mathematical function
Can you suggest me a way on how to rectify the same.
Best,
Shantanu