

# Two sample MR report

---

## order.Actinomycetales.id.420 against BK polyomavirus seropositivity

Date: **22 January, 2026**

---

### Results from two sample MR:


|method                    | nsnp|         b|        se|      pval|
|:-------------------------|----:|---------:|---------:|---------:|
|MR Egger                  |    5| 1.0585749| 1.0425854| 0.3847092|
|Weighted median           |    5| 0.7629380| 0.5196908| 0.1420876|
|Inverse variance weighted |    5| 0.9386928| 0.4015371| 0.0194003|
|Simple mode               |    5| 0.7983199| 0.6996623| 0.3175400|
|Weighted mode             |    5| 0.7559855| 0.6390038| 0.3022937|

---

### Heterogeneity tests


|method                    |        Q| Q_df|    Q_pval|
|:-------------------------|--------:|----:|---------:|
|MR Egger                  | 1.963028|    3| 0.5801157|
|Inverse variance weighted | 1.978552|    4| 0.7397039|

--- 

### Test for directional horizontal pleiotropy


| egger_intercept|        se|      pval|
|---------------:|---------:|---------:|
|      -0.0131961| 0.1059103| 0.9087228|

--- 

### Test that the exposure is upstream of the outcome


| snp_r2.exposure| snp_r2.outcome|correct_causal_direction | steiger_pval|
|---------------:|--------------:|:------------------------|------------:|
|        3.79e-05|          3e-07|TRUE                     |    0.7298839|

Note - R^2 values are approximate

---

### Forest plot of single SNP MR

![plot of chunk chunk5](figure/orderActinomycetalesid420_against_BK_polyomavirus_seropositivitychunk5-1.png)

---

### Comparison of results using different MR methods

![plot of chunk chunk6](figure/orderActinomycetalesid420_against_BK_polyomavirus_seropositivitychunk6-1.png)

---

### Funnel plot

![plot of chunk chunk7](figure/orderActinomycetalesid420_against_BK_polyomavirus_seropositivitychunk7-1.png)

---

### Leave-one-out sensitivity analysis

![plot of chunk chunk8](figure/orderActinomycetalesid420_against_BK_polyomavirus_seropositivitychunk8-1.png)
